Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,017 C→T 23.4% 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,0170CT23.4% 369.2 / 96.7 185intergenic (‑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 (61/68);  new base T (38/2);  total (115/70)
Fisher's exact test for biased strand distribution p-value = 9.99e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 188 total aligned reads displayed.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116
                                                                                                    |                                                                                                            
   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                                                           >  2:1342879/1‑100 (MQ=255)
     aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTgn                                                                                                         <  2:2788830/100‑2 (MQ=255)
      ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                                                                        >  2:1827001/1‑100 (MQ=255)
      ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                                                                        >  2:175748/1‑100 (MQ=255)
       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                                                                      >  2:1575285/1‑100 (MQ=255)
         ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga‑a                                                                                                   >  2:1524677/1‑100 (MQ=255)
         ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga‑a                                                                                                   >  2:1101737/1‑100 (MQ=255)
            ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At                                                                                                  <  1:3301264/98‑1 (MQ=255)
             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑t‑‑‑t‑t                                                                                          >  2:771530/1‑100 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑g                                                                                      <  2:2497041/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑g                                                                                      <  1:2202870/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑g                                                                                      <  1:1891880/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑g                                                                                      <  1:771528/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑g                                                                                      <  1:1736690/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑g                                                                                      <  1:155391/100‑1 (MQ=255)
                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTa                                                                                    >  2:438803/1‑100 (MQ=255)
                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTa                                                                                    >  1:1176200/1‑100 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGcc                                                                                >  1:2951213/1‑100 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGcc                                                                                >  2:2848465/1‑100 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGcc                                                                                >  2:3092866/1‑100 (MQ=255)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTg                                                                              >  2:1685756/1‑100 (MQ=255)
                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCgg                                                                      <  1:3365549/100‑1 (MQ=255)
                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGa                                                                     >  2:1429063/1‑100 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                  <  1:1406199/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                  >  1:2806088/1‑100 (MQ=255)
                                    attttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑C‑TT‑‑GTAGGCCTGATACGCGTAgcgc                                                                 >  1:3193611/10‑100 (MQ=17)
                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                <  1:751435/100‑1 (MQ=255)
                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                <  2:2872182/100‑1 (MQ=255)
                                      tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑T‑‑‑T‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATc                                                              >  2:3131942/8‑100 (MQ=21)
                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                  <  2:2036993/97‑1 (MQ=255)
                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                               >  2:1724467/1‑100 (MQ=255)
                                       ttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCA‑GA‑AA‑T‑‑‑C‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATc                                                              >  2:2922567/7‑100 (MQ=18)
                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                              >  2:2016511/1‑100 (MQ=255)
                                        tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑T‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCATAGCGCATCa                                                             >  2:2954177/6‑100 (MQ=18)
                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                           >  2:910925/4‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑AT‑C‑‑‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:223702/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GG‑GA‑C‑‑‑T‑AT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                                                          >  2:831/3‑100 (MQ=18)
                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                          <  1:1201060/100‑1 (MQ=255)
                                            ngCCGGATAAGGCGTTCGCGCCGCATCCGGCCTAC‑AAAAGTTTGCAAATTCAATAAATTGCA‑GA‑AT‑T‑‑‑C‑AT‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa                                                         >  1:437202/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑A‑‑‑T‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                         >  2:2904754/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑A‑‑‑T‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                         >  1:2454328/2‑100 (MQ=21)
                                             gtcggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCA‑GG‑AA‑C‑‑‑G‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                        >  1:236624/10‑99 (MQ=21)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                        >  2:60731/1‑100 (MQ=255)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGAGAATCAGGCAt                                                        >  2:280631/1‑100 (MQ=255)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                       >  2:721850/1‑100 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAAC‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                         >  1:1845889/1‑100 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑C‑‑‑T‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                        >  2:58121/1‑99 (MQ=18)
                                              ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                      >  1:2737739/9‑100 (MQ=18)
                                               cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagtt                                                      >  2:3009341/8‑97 (MQ=17)
                                               cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                      >  2:3400471/1‑100 (MQ=255)
                                                ggataaggcgTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCA‑GA‑AA‑T‑‑‑C‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                     >  1:115540/7‑100 (MQ=18)
                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                    >  1:1966237/7‑98 (MQ=18)
                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGTAAATTCAATACATTGCA‑TG‑GG‑C‑‑‑C‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                     >  1:41950/7‑100 (MQ=17)
                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAAC‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                                                      >  1:2025235/7‑97 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:235105/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑A‑‑‑T‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                     >  2:1084949/7‑96 (MQ=25)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑T‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                                                     >  2:2506337/7‑96 (MQ=21)
                                                 gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑‑‑C‑TG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaagctt                                                    >  2:982579/6‑95 (MQ=21)
                                                 gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑T‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttt                                                    >  2:147905/6‑100 (MQ=17)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                   <  2:2055126/100‑1 (MQ=255)
                                                   taaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑C‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttagc                                                  >  2:3190606/4‑97 (MQ=17)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                  >  2:2198868/1‑100 (MQ=255)
                                                    aaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcgccc                                                >  1:2443099/3‑94 (MQ=21)
                                                     cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                >  1:532750/1‑100 (MQ=255)
                                                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                               >  2:636271/1‑100 (MQ=255)
                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                            <  1:2766626/100‑1 (MQ=255)
                                                          ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                           <  2:798264/100‑1 (MQ=255)
                                                                ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                     <  1:1766074/100‑1 (MQ=255)
                                                                ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                     >  1:2741244/1‑100 (MQ=255)
                                                                ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                     >  2:2131044/1‑100 (MQ=255)
                                                                 cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                    >  2:2753732/1‑100 (MQ=255)
                                                                  ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TT‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                   >  2:3385853/1‑100 (MQ=255)
                                                                  ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                   <  2:1032330/100‑1 (MQ=255)
                                                                    aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                 <  2:809253/100‑1 (MQ=255)
                                                                      ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGATCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                               <  1:1826993/100‑1 (MQ=255)
                                                                       cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              <  1:615911/100‑1 (MQ=255)
                                                                       cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              <  2:2653136/100‑1 (MQ=255)
                                                                           cTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                          <  2:2203091/100‑1 (MQ=255)
                                                                            tAC‑AAAATCGTGCAAATTCAATACATTGCT‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGTATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                         <  1:2061347/100‑1 (MQ=255)
                                                                            tAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                         <  1:1332952/100‑1 (MQ=255)
                                                                             aC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTCTTTTGTTATCAGTCTTGTATGGCttt                        >  1:1666071/1‑100 (MQ=25)
                                                                              c‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                       >  1:1070762/1‑100 (MQ=255)
                                                                                 aaaTCGTTCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                     <  1:2152895/100‑1 (MQ=255)
                                                                                  aaTCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                    <  2:2365425/100‑1 (MQ=255)
                                                                                    tCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                  <  1:355783/100‑1 (MQ=255)
                                                                                     cGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                 >  1:39255/1‑100 (MQ=255)
                                                                                      gTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                >  1:251217/1‑100 (MQ=255)
                                                                                       tGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  1:1129711/100‑1 (MQ=255)
                                                                                       tGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               <  1:2250687/100‑1 (MQ=255)
                                                                                       tGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               >  1:1823842/1‑100 (MQ=255)
                                                                                        gCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              <  1:1069124/100‑1 (MQ=255)
                                                                                        gCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              <  1:2040945/100‑1 (MQ=255)
                                                                                        gCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              <  2:1150607/100‑1 (MQ=255)
                                                                                         cAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa             <  2:3322459/100‑1 (MQ=255)
                                                                                         cAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa             >  2:2472563/1‑100 (MQ=255)
                                                                                          aaaTTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt            <  1:3178742/100‑1 (MQ=255)
                                                                                           aaTTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa           <  1:1602566/100‑1 (MQ=255)
                                                                                              tCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag        <  1:287265/100‑1 (MQ=255)
                                                                                               cAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga       <  1:522501/100‑1 (MQ=255)
                                                                                                    cATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCTGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGGATGGCTTTTAATTATTCGATAGAGACGtt   <  1:696498/99‑1 (MQ=255)
                                                                                                    |                                                                                                            
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116

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: ‑