Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 31.6% intergenic (‑258/+50) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,706,0480GT31.6% 382.0 / 173.8 201intergenic (‑258/+50)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (62/71);  new base T (55/3);  total (127/74)
Fisher's exact test for biased strand distribution p-value = 1.72e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

Only 100 of 202 total aligned reads displayed.

GCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705950‑3706146
                                                                                                            |                                                                                                  
gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                                                                                   >  1:1134815/1‑100 (MQ=255)
gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                                                                                   >  2:1992020/1‑100 (MQ=255)
 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                                                  <  1:2282282/100‑1 (MQ=255)
 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                                                                                  <  1:2875506/100‑1 (MQ=255)
  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGACTT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 >  1:2112048/1‑100 (MQ=255)
  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 >  1:1899801/1‑100 (MQ=255)
  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 >  1:458707/1‑100 (MQ=255)
   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                <  1:349317/100‑1 (MQ=255)
    attttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAgcgc                                                                                               >  2:796438/10‑100 (MQ=18)
    attttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                                                                                               >  2:1839977/10‑100 (MQ=18)
     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                                              <  2:2088136/100‑1 (MQ=255)
     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                                              >  2:582712/1‑100 (MQ=255)
      ttcgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCGT‑‑ACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                                             >  2:2311783/8‑100 (MQ=18)
       ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATc                                                                                            >  1:2351461/7‑100 (MQ=18)
       ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                           >  2:2138085/7‑100 (MQ=21)
       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                            <  1:2447633/100‑1 (MQ=255)
        tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATc                                                                                            >  1:1920233/6‑100 (MQ=21)
          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                                                                                         >  2:2291973/4‑100 (MQ=21)
          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                         >  2:2519674/4‑100 (MQ=21)
          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                         >  2:1396946/4‑100 (MQ=21)
          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                         <  1:1396944/97‑1 (MQ=21)
          taggCCGGATAAGGCGTTAACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                         >  1:1123923/4‑100 (MQ=18)
           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  2:889574/3‑100 (MQ=17)
           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  2:1880716/3‑100 (MQ=17)
           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  1:2129657/3‑100 (MQ=17)
           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  1:1492369/3‑100 (MQ=17)
           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAACCTGCAAATTCAATATATTGCAGGATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  1:150656/3‑100 (MQ=17)
           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAACCGCGCAAATTCAATATATTGCGGAGAA‑‑‑AAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                                                                                       >  1:2602641/3‑99 (MQ=21)
            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                        >  2:2212556/2‑100 (MQ=21)
            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:2741719/2‑100 (MQ=18)
            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAACTTTGCAAATTCAATATATTGCAATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAgatc                                                                                       >  1:1402919/2‑97 (MQ=18)
            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:297573/2‑100 (MQ=17)
             gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  2:60616/1‑100 (MQ=255)
             gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  1:2211422/1‑100 (MQ=255)
             gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  2:1317900/1‑100 (MQ=255)
             gCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATGCAATACATTGCAGAATT‑‑‑TTC‑‑CTAGGCCTGATAAGCGGTGCGCATCAGGCAt                                                                                      >  2:1143506/1‑100 (MQ=255)
             gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                      >  2:308360/1‑99 (MQ=17)
             gCCGGATAAGGCGTTCGCGCCGCATCCGGCCT‑‑ACAAAAGTTTGCAAATTCAATAAATTGCAGAATT‑‑‑CAT‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAt                                                                                      >  2:2223264/1‑100 (MQ=17)
             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                      >  2:715891/1‑99 (MQ=18)
             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                                      >  2:1764218/1‑99 (MQ=25)
             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                                      >  1:1629349/1‑99 (MQ=21)
             gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                      >  2:2513525/1‑99 (MQ=21)
              ccggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  2:191299/9‑98 (MQ=17)
              ccggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  1:561257/9‑98 (MQ=17)
              ccggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                                                     >  1:1658508/9‑98 (MQ=21)
              ccggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                      >  1:1681024/9‑98 (MQ=17)
              ccggataaggcgTTAACGCCGCATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                                     >  2:278676/9‑100 (MQ=18)
              ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                                     >  1:849248/1‑100 (MQ=255)
              ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                                                     >  1:1448217/1‑100 (MQ=255)
              ccGGATAAGGCGTTAACGCCGCATCCGGCATAAACAAAATCGTGCAAAATCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  1:2065620/1‑100 (MQ=17)
               cggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttat                                                                                   >  2:185518/8‑98 (MQ=18)
               cggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                                                                                    >  2:2898489/8‑100 (MQ=17)
                ggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  2:1379629/7‑100 (MQ=18)
                ggataaggcgTTTACGCCGCATCCGGCCT‑‑ACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:1188413/7‑100 (MQ=17)
                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  2:2682618/7‑100 (MQ=18)
                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:1979177/7‑100 (MQ=17)
                ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttta                                                                                   >  1:754814/7‑98 (MQ=21)
                ggataaggcgTTCACGCCGCATCCGACCC‑‑ACAAAGTCTTGCAAAATCAATAAATTGCAGGAACG‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                                                                                   >  2:88122/7‑96 (MQ=25)
                 gATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                                                  >  2:1616832/1‑100 (MQ=255)
                   tGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttttt                                                                                >  1:1390426/1‑98 (MQ=255)
                    gcggcgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                               >  1:1199982/1‑100 (MQ=255)
                    gcggcgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                               <  2:1448220/100‑1 (MQ=255)
                     cggcgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                              >  2:390814/1‑100 (MQ=255)
                        cgTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                           <  2:2413054/100‑1 (MQ=255)
                           tAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                                        <  2:463066/100‑1 (MQ=255)
                             aCGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                                      <  2:2587543/100‑1 (MQ=255)
                               gCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                                                    <  1:2565758/100‑1 (MQ=255)
                                  tcttccgatcT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                                                 <  2:2529070/91‑1 (MQ=255)
                                    aTCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                               <  1:595569/100‑1 (MQ=255)
                                      ccGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                             <  2:2403941/100‑1 (MQ=255)
                                      ccGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                                             >  1:427177/1‑100 (MQ=255)
                                       cGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGAGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                                                 >  1:109915/1‑95 (MQ=255)
                                        ggCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                           >  1:184531/1‑100 (MQ=255)
                                            t‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       >  1:2150628/1‑100 (MQ=255)
                                            t‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       >  1:2408914/1‑100 (MQ=255)
                                                 aaaaaTCGTGCAAATTCGATACATTGCAGAATT‑‑‑TTC‑‑GGAGGCCTGAGAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                                      <  1:2866558/98‑1 (MQ=255)
                                                  aaaaTCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                                                    <  1:882415/100‑1 (MQ=255)
                                                    aaTCGTGCAAATTCGATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                  <  1:2845744/100‑1 (MQ=255)
                                                    aaTCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                  <  2:2359571/100‑1 (MQ=255)
                                                       cGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                                               >  2:745275/1‑100 (MQ=255)
                                                         tGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             >  1:1800363/1‑100 (MQ=255)
                                                          gCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            <  1:289808/100‑1 (MQ=255)
                                                           cAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                           >  2:1646894/1‑100 (MQ=255)
                                                            aaaTTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                          <  1:904804/100‑1 (MQ=255)
                                                            aaaTTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                          >  2:398980/1‑100 (MQ=255)
                                                              aTTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                                        >  1:2015610/1‑100 (MQ=255)
                                                                       aTTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTGTGCCATCAGTCTTGTATGGCTGTTAATTATTCGATAGAGACGttt                                >  1:2467596/1‑99 (MQ=255)
                                                                        ttGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTTCTTTTTTATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                              >  1:268705/1‑100 (MQ=38)
                                                                        ttGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                              <  2:1714930/100‑1 (MQ=255)
                                                                          gCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            >  2:2011583/1‑100 (MQ=255)
                                                                            aGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  1:1521932/100‑1 (MQ=255)
                                                                            aGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  2:830496/100‑1 (MQ=255)
                                                                              gAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          >  1:615295/1‑99 (MQ=255)
                                                                                aTT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                        <  2:1199984/99‑1 (MQ=255)
                                                                                 tt‑‑‑ttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGACACGTTTTCGAAGTGAt                      >  1:1702270/1‑100 (MQ=255)
                                                                                       ttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGt                    <  1:1290452/100‑1 (MQ=255)
                                                                                            gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 >  2:207370/1‑100 (MQ=255)
                                                                                              aGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  1:1036274/100‑1 (MQ=255)
                                                                                                gCCTGATAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTaa            >  1:37437/1‑100 (MQ=255)
                                                                                                    gATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         >  2:1223202/1‑100 (MQ=255)
                                                                                                            |                                                                                                  
GCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑ACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705950‑3706146

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: