Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,017 C→T 24.3% intergenic (‑227/+81) 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,0170CT24.3% 554.6 / 161.6 289intergenic (‑227/+81)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (94/108);  new base T (63/2);  total (179/110)
Fisher's exact test for biased strand distribution p-value = 4.61e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 292 total aligned reads displayed.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116
                                                                                                     |                                                                                                        
  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACAt                                                                                                        >  1:224486/1‑100 (MQ=255)
   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACAtt                                                                                                       <  2:2602896/100‑1 (MQ=255)
    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTg                                                                                                      <  2:2369967/100‑1 (MQ=255)
     aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGc                                                                                                     >  2:3288547/1‑100 (MQ=255)
          gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAACGTGCAAA‑TTCAATACATTGCA‑GA‑At                                                                                              <  1:3861547/100‑1 (MQ=255)
           caatcAGATTGCTGACAACGTGCGCGTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAC‑AAAGTCGTATAAACTGCAATACATTgaattta                                                                                                <  1:2977455/96‑7 (MQ=17)
             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑At‑ttt                                                                                          >  1:1164220/1‑100 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑g                                                                                      <  2:2109352/100‑1 (MQ=255)
                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTa                                                                                    <  2:1379094/100‑1 (MQ=255)
                   atcttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAgg                                                                                  <  1:1686241/97‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGcc                                                                                >  1:727115/1‑100 (MQ=255)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAggagat                                                                              >  2:1686246/1‑96 (MQ=255)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTg                                                                              <  1:2007268/100‑1 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGa                                                                             <  1:3536186/100‑1 (MQ=255)
                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATa                                                                           <  1:1218830/100‑1 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATaa                                                                          >  1:3012156/1‑98 (MQ=255)
                              gTGCGCTTTGTTTATGCCGGATGTGGCGTTAACTCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCg                                                                       <  2:1842329/100‑1 (MQ=255)
                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCg                                                                       >  2:2175300/1‑100 (MQ=255)
                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                                                                      >  1:2823462/1‑100 (MQ=255)
                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                   <  1:1911322/100‑1 (MQ=255)
                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                   >  2:4162666/1‑100 (MQ=255)
                                   cTTTGTTTCTGCCGGATGCTGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GC‑AT‑TTTC‑‑GTAGGCCTGATAAGTGGAgc                                                                   <  1:2049830/99‑1 (MQ=17)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                  >  2:289564/1‑100 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                  <  2:1782461/100‑1 (MQ=255)
                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                >  2:3497768/1‑100 (MQ=255)
                                      tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAACATATTGCA‑AT‑TC‑TCTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                               >  2:27686/8‑100 (MQ=18)
                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                               >  1:2553322/1‑100 (MQ=255)
                                        tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAA‑TTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                             >  2:2419993/6‑100 (MQ=18)
                                        tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAA‑TTCAATAGATTGCA‑GA‑GA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                             >  1:639892/6‑100 (MQ=17)
                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                            >  1:1605099/1‑100 (MQ=255)
                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAATTC‑TGCAAA‑TTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                           >  2:679164/5‑100 (MQ=18)
                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAA‑TTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                           >  2:251099/4‑100 (MQ=18)
                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAA‑TCCAAAATATTGCA‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:3457958/4‑100 (MQ=21)
                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAGTCGTGCAAA‑TTCACTACATTGAG‑GA‑AT‑TTAC‑‑GTAGGCCTGATAAGGGGAGCGCATCAgg                                                           >  1:861178/1‑100 (MQ=17)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAACATATTGC‑‑‑A‑AT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:2839328/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAA‑CTCAATATATTGCA‑GAGAT‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:1582414/3‑100 (MQ=17)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAA‑TTCAATATATTGCA‑GG‑GA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                                                          >  1:2434933/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                          >  2:3035489/3‑100 (MQ=18)
                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                          <  1:1432903/100‑1 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  2:1910432/1‑100 (MQ=255)
                                            ggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAA‑TTCAATATATTGCA‑GA‑AA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:1547835/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAA‑TTCAATATATTGCA‑GA‑AA‑CCTT‑‑CTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:2564641/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAACATATTGC‑‑‑A‑AT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGct                                                         >  2:598624/2‑99 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAA‑TTCAATATATTGTG‑GG‑CT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:1959008/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAA‑TTCAATATATTGCA‑GG‑AA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa                                                         >  2:2077125/2‑100 (MQ=18)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                        >  2:4083402/1‑100 (MQ=255)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                        >  2:3743114/1‑100 (MQ=255)
                                             gCCGGATAAGGCGTTTACGCCTTATCCGGCCTAC‑AAAATCGTGCTAA‑TTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                        >  2:4055131/1‑99 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAA‑TTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                        >  2:2704278/1‑99 (MQ=18)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAA‑TTCAATGTATTGCA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                        >  2:3164096/1‑99 (MQ=25)
                                               cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAA‑TTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                                                      >  1:1786424/8‑100 (MQ=17)
                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAA‑TTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                    >  2:2957040/7‑98 (MQ=21)
                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAA‑TTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:3491711/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAA‑TTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:1501097/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAA‑TTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:2472738/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAA‑TTCAATGTATTGCA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                     >  2:1280973/7‑96 (MQ=25)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAA‑TTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                                                     >  1:3155307/7‑96 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAA‑TTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattc                                                     >  2:1327207/7‑99 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAA‑TTCAATAGATTGCA‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                     >  1:101432/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAA‑TCCAAAATATTGCA‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                    >  2:1149377/7‑98 (MQ=25)
                                                 gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAA‑TTCAATAAATTGAT‑CG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                    >  1:1486924/6‑100 (MQ=18)
                                                 gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAA‑TTCAATATATTGCAGGG‑AC‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctta                                                    >  2:3466834/6‑95 (MQ=21)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                   <  2:2629454/100‑1 (MQ=255)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                   >  1:1841852/1‑100 (MQ=255)
                                                   taaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAA‑TTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc                                                  >  1:4078602/4‑100 (MQ=17)
                                                   taaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAA‑TTCAATATATTGCA‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAAGCAGTtttgc                                                  >  2:1971434/4‑100 (MQ=17)
                                                    gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                 <  2:3906001/100‑1 (MQ=255)
                                                    gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACACTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                 <  1:3745252/100‑1 (MQ=255)
                                                       gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt                                              >  2:335722/1‑100 (MQ=255)
                                                        cgTTCACGCCTTATCCGGCCTAC‑AAAATCATGCAAA‑CTCAATACATTGCC‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAACGCATCAGGCATTtttgcct                                                >  1:228780/1‑97 (MQ=17)
                                                        cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTttttccttt                                              >  2:429108/1‑99 (MQ=255)
                                                        cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                             <  1:2889977/100‑1 (MQ=255)
                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          <  1:1347015/100‑1 (MQ=255)
                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATGGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          <  1:2116188/100‑1 (MQ=255)
                                                            accGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                         >  2:2888554/3‑100 (MQ=255)
                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                        >  1:3835576/1‑100 (MQ=255)
                                                                 cTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                    <  2:1605103/100‑1 (MQ=255)
                                                                      ccGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                               >  1:3990390/1‑100 (MQ=255)
                                                                      acGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGACCTGCTAAGCGAAACGCATCAGGCATTTTTGCCTTTGCTATCAGTCTTGt                                >  1:3940670/2‑99 (MQ=21)
                                                                       cGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              <  2:2564645/100‑1 (MQ=255)
                                                                        ggCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGACCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                             >  1:3059117/1‑100 (MQ=255)
                                                                         gacTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                             <  1:3591753/97‑1 (MQ=255)
                                                                         gCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTCAGTCTTGTATgg                            >  1:692630/1‑100 (MQ=18)
                                                                         gCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                            >  2:3700401/1‑100 (MQ=255)
                                                                          ccTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                           <  2:3862681/100‑1 (MQ=255)
                                                                          ccTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                           <  1:335720/100‑1 (MQ=255)
                                                                          ccTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                           <  2:2811919/100‑1 (MQ=255)
                                                                             aC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                        >  2:532106/1‑100 (MQ=255)
                                                                             aC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                        <  2:908367/100‑1 (MQ=255)
                                                                              c‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCGGGCATTTTTGCCTTTGCCATCAGTATTGTATGGCtttt                       <  1:2630123/100‑1 (MQ=255)
                                                                                    tCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                  <  1:1669262/100‑1 (MQ=255)
                                                                                       tGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               >  2:1997054/1‑100 (MQ=255)
                                                                                        gCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              <  2:790905/100‑1 (MQ=255)
                                                                                          aaa‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCTGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  1:2739882/97‑1 (MQ=255)
                                                                                                cAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga       >  2:1524824/1‑100 (MQ=255)
                                                                                                 aaTACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga       >  1:415643/1‑99 (MQ=255)
                                                                                                 aaTACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCCATagag        <  1:3725661/98‑1 (MQ=255)
                                                                                                  aTACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAtttcgatc     >  2:3787735/1‑93 (MQ=255)
                                                                                                     cATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt  <  2:3247635/100‑1 (MQ=255)
                                                                                                     cATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt  <  2:755463/100‑1 (MQ=255)
                                                                                                     |                                                                                                        
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAA‑TTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116

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