Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,486,502 C→T T128I (ACT→ATT)  crp → DNA‑binding transcriptional dual regulator CRP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,486,5020CT89.8% 115.9 / 6.4 49T128I (ACT→ATT) crpDNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (23/21);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01

TGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGG  >  NC_000913/3486271‑3486742
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
tGATCAAAGACGCAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTACACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTACAAGTCATTTCAGAGAAAGTGGGCAAc                                                                                                                                                                                                                               >  2:336534/1‑251 (MQ=255)
  aTCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGTCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGTCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCt                                                                                                                                                                                                                             >  2:140012/1‑251 (MQ=255)
        gACGAAGGGGGTAAAGAAATGATCCTCCCCTAGCGTAATCAGGGGGTTTTTATTGGCGAAGTGGGCCTGTGTGAAGAGGCCCAGGAAGGTAGCGCATGGGTTCGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGtt                                                                                                                                                                                                                       <  2:187669/251‑1 (MQ=255)
          cGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCCTTTCAGAGAACGTGGGCAACCTGGCGTTcc                                                                                                                                                                                                                     >  2:480736/1‑251 (MQ=255)
           gAAGAGGGCAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGCGGGCCAGGAACGTAGCGCTTGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCAAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCCCAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGGAAGTGGGCAACCTGGCGTTCCt                                                                                                                                                                                                                    >  2:255009/1‑251 (MQ=255)
            aagggggTAAAGAAATACCCCTCTCCTATCTAAATGAGGGTGTTTTTATGGGCGAACGGGGCTTGTTTAAAGAGGGCAAGGAAGGTAGCCCTTGGGTTCGCGAAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTc                                                                                                                                                                                                                   <  1:475521/247‑1 (MQ=255)
               agGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAATGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGAc                                                                                                                                                                                                                >  2:31673/1‑251 (MQ=255)
                                           aaattggggtttttttATTGGCAAAAGGGGCCTGTTTGAAGGGGCTCGGGAACGTAGCACATGGGGAGTGGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAAGTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTctgc                                                                                                                                                                                    <  2:104569/239‑1 (MQ=255)
                                            aaTCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTctgct                                                                                                                                                                                   >  1:248326/1‑251 (MQ=255)
                                             aTCAGGGGGATTTTATTGGCGAACTGGGCCTGTTTGAAGATGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCCCAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGACTTGCACAGACTctgctg                                                                                                                                                                                  >  2:15845/1‑251 (MQ=255)
                                                     gATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGCGAAAGTGGGCAACCTGGCGTTCCt                                                                                                                                                                                                                    >  2:257338/1‑209 (MQ=255)
                                                       ttttATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTAAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCa                                                                                                                                                                         >  1:150880/1‑250 (MQ=255)
                                                                     tGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAACACAACCAGACGcc                                                                                                                                                          >  2:220333/1‑250 (MQ=255)
                                                                            gTTTAAAGAGGGCCGTGACCGTAGCGCTTGGGTACGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACt                                                                                                                                                    <  1:228117/250‑1 (MQ=255)
                                                                                  agagGGCCAGGAACGCAGCGCATGGGTAAGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAACACCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACATGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGATGTATCTGGCAAAACAACCCGACGCTATGACTCACCCgg                                                                                                                                             >  2:500701/1‑251 (MQ=255)
                                                                                            ggtccctgccaatctggTACGTGAGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATg                                                                                                                                     <  1:222796/235‑1 (MQ=255)
                                                                                             gAACGTAGCGCATGGGTACGTGAGAAAAACGCCTGTGAAGTGGCTGAACTTTCGTACAAAAGATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGGTTGTATGCACCGATGGCGCGTCGTCTGCAAGTCATTTCAGAGACAGTGGGCAAGCTGGCGTGCCTCAACGTGACGGGCCGCCTTGCACAGACTATTGTGGATATGGCAATACAACCAGACGCTATGACACACCCGGACGGTATGc                                                                                                                                    >  2:299199/1‑250 (MQ=255)
                                                                                                       cATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAACATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAAtt                                                                                                                         >  2:490969/1‑251 (MQ=255)
                                                                                                           gggaCGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTAccc                                                                                                                     <  1:267736/248‑1 (MQ=255)
                                                                                                           ggaaCGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTAccc                                                                                                                     >  1:393052/4‑251 (MQ=255)
                                                                                                           ggTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGCTTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCAAAGACTCTGCTGAATCTGGCAAAACAACCAGAGGCTATGACTCACCCGGACGGTATg                                                                                                                                     >  2:222796/1‑235 (MQ=255)
                                                                                                            gaaatgggaaaaaaCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGGGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTGCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAAGTAccc                                                                                                                     <  2:465658/241‑1 (MQ=255)
                                                                                                             tACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGt                                                                                                                   >  2:28939/1‑251 (MQ=255)
                                                                                                                 ggCGAAATCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGTCATTCTGATGCGTTTGTCTTCACAGATGGCGCGTCGTCTTCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAgg                                                                                                               <  1:288299/250‑1 (MQ=255)
                                                                                                                     gAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAACAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTTAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCCCAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGGATGCACATCCAAATTACCCGTCAGGaa                                                                                                             >  2:451851/1‑250 (MQ=255)
                                                                                                                     gAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATc                                                                                                                               >  2:361553/1‑232 (MQ=255)
                                                                                                                     aaaaaCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATc                                                                                                                               <  1:361553/231‑1 (MQ=255)
                                                                                                                      aaaaCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGAATGCAAATCAAAATTACCCGTCAGGCAAt                                                                                                           >  1:383487/1‑251 (MQ=255)
                                                                                                                           cGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTc                                                                                                      >  2:114715/1‑251 (MQ=255)
                                                                                                                                     gTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTAAATCTGGCAAAACAACCAGACGCTATGACACACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCgg                                                                                            <  1:388586/251‑1 (MQ=255)
                                                                                                                                       ggCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGc                                                                                           <  1:114715/250‑1 (MQ=255)
                                                                                                                                           gAAATTGCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATGCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGtt                                                                                       <  1:500701/250‑1 (MQ=255)
                                                                                                                                            aaaTTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTtct                                                                                     <  1:255009/251‑1 (MQ=255)
                                                                                                                                             aaTTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGAGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTtctc                                                                                    >  1:191270/1‑251 (MQ=255)
                                                                                                                                                      cAAAAAAGTTCGGCAAGTGATTCATGTAAACCCGGACATTCTTACGCGTTTGTCGGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGTGAAAGTGGGCAACCTTGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACATCCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTTTTCGCGTGAAACCg                                                                           <  1:369537/251‑1 (MQ=255)
                                                                                                                                                        aaaaaTTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTg                                                                         <  1:336534/251‑1 (MQ=255)
                                                                                                                                                                  ccAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGActct                                                                                                                                                                                      <  1:85762/132‑1 (MQ=255)
                                                                                                                                                                  ccAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGActct                                                                                                                                                                                      >  2:85762/1‑132 (MQ=255)
                                                                                                                                                                    aaTTGATTCATGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCCCCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTg                                                             <  2:101411/251‑1 (MQ=255)
                                                                                                                                                                         aTTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTgc                                                            >  2:1968/1‑246 (MQ=255)
                                                                                                                                                                           tcaagTAAACCCGGAACTGCTGAGGTGTTTGTCGGCCCAGATGGCGCAGCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGTTAAGATTGTCGGCTGTTCTCGTGAAACCGTGGGCCGCATTCTGAAGATg                                                       <  2:142317/246‑1 (MQ=255)
                                                                                                                                                                                      cGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGGCGTCTGCAAGTCATTTCAGAGAAAGGGTGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATTACTCACCCGGATGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGGCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTTGAAGATCAg                                           <  1:50168/251‑1 (MQ=255)
                                                                                                                                                                                        gACATTCTTATGCGTTGGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACCCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGGTCTCGTGAAACCGTGGGTCGCATTCTGAAGATGCTGGAAGATCAGaa                                         <  1:145590/251‑1 (MQ=255)
                                                                                                                                                                                        gACATTCTGATGCGTTTGTCTGCACAGATGGGGCGTCGGCTGCAAGACATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGaa                                         <  2:157869/251‑1 (MQ=255)
                                                                                                                                                                                             tCTGATGCGTTTGTCTGCACAGCTGGCGCGTCGTCTGCAAGTCATTTCAGCGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCCAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGACATTGGTCAGATTGTCGGCTGTTCTCGTGAAAACGTGGGACGAATTCTGAAGATGCTGGACGATCAGAACCTg                                     >  2:60338/1‑250 (MQ=255)
                                                                                                                                                                                                                 aGATGGCGCGTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAAcc                >  1:234983/1‑251 (MQ=255)
                                                                                                                                                                                                                       ggcgTCGTCTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCAtcgtc          <  2:86553/250‑1 (MQ=255)
                                                                                                                                                                                                                              tcTGCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACg   >  1:170397/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gCAAGTCATTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACgg  <  1:429214/249‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
TGATCAAAGACGAAGAGGGTAAAGAAATGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTCTGCAAGTCACTTCAGAGAAAGTGGGCAACCTGGCGTTCCTCGACGTGACGGGCCGCATTGCACAGACTCTGCTGAATCTGGCAAAACAACCAGACGCTATGACTCACCCGGACGGTATGCAAATCAAAATTACCCGTCAGGAAATTGGTCAGATTGTCGGCTGTTCTCGTGAAACCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGG  >  NC_000913/3486271‑3486742

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: