Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 24.8% 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,0480GT24.8% 611.8 / 182.0 284intergenic (‑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 (101/107);  new base T (67/2);  total (175/109)
Fisher's exact test for biased strand distribution p-value = 2.61e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146
                                                                                                            |                                                                                                  
tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCgg                                                                                                    <  1:3151016/100‑1 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 <  1:2370531/100‑1 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                 >  2:2706405/1‑100 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                >  1:1099675/1‑100 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                >  1:1172899/1‑100 (MQ=255)
      atttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCAGA‑AA‑T‑C‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCa                                                                                              >  1:2628907/9‑100 (MQ=18)
       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                                                             >  2:1263790/1‑100 (MQ=255)
        gTTTATGCCGGATGCGGCGTTAACGCCTTCTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                                                            >  2:3658973/1‑100 (MQ=255)
         tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGA‑GT‑T‑C‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                           >  1:691795/6‑100 (MQ=14)
          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATAGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                          <  2:2414482/100‑1 (MQ=255)
          gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCAAT‑CT‑T‑C‑‑‑TG‑‑GTAGGCCTGATAAGCATCGCGCATCAg                                                                                          >  1:19611/5‑100 (MQ=18)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑AC‑C‑A‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                          >  1:167847/5‑100 (MQ=18)
           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑CAAACTTTGCAAATTCAATATATTGC‑‑A‑AT‑C‑T‑‑‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                         >  2:2095046/4‑100 (MQ=18)
            nggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCAGA‑AA‑T‑C‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                        >  2:3164560/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCAGG‑AA‑A‑C‑‑‑AC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc                                                                                        >  1:225326/3‑100 (MQ=18)
             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  2:3561352/1‑100 (MQ=255)
             tGCCGGATGCGGCGTTAACGCCTTATACGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAAGCCTGATAAGCGGAGCGCATCAGGCa                                                                                       >  2:290824/1‑100 (MQ=255)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATGTTGCAGA‑AA‑C‑C‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:2504201/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGATCG‑AT‑C‑C‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                       >  2:2267474/2‑100 (MQ=21)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  2:1052906/1‑100 (MQ=255)
              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  1:1622948/1‑100 (MQ=255)
              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑AT‑CCC‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                                                     >  1:3468936/1‑100 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑AC‑C‑A‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                      >  2:2806682/1‑99 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACC‑AT‑C‑C‑‑‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                      >  1:2550234/1‑99 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCAGAGAT‑T‑A‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                                      >  2:1964773/1‑99 (MQ=17)
               caggataaggcgTTCACGCCGCATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTGCAGA‑GA‑T‑G‑‑‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                                                     >  1:2072162/9‑100 (MQ=14)
                cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑ATCC‑C‑‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                                                                                   >  1:2950735/8‑99 (MQ=18)
                cggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑AC‑C‑A‑‑‑T‑‑‑GTAAACCTGATAAGCGTAGCGCATCAGGCAgtt                                                                                    >  1:2114005/8‑100 (MQ=18)
                cggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTGGG‑CT‑C‑C‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                                                                                    >  2:1017215/8‑100 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGA‑GT‑T‑C‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  2:3331247/7‑100 (MQ=17)
                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGCAAT‑TC‑AAC‑‑‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCAatt                                                                                    >  1:1546523/7‑100 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑A‑AT‑T‑CAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                                                                                    >  1:440205/7‑97 (MQ=21)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑AC‑C‑A‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  2:393529/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑AC‑C‑A‑‑‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                   >  1:1605903/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTGGG‑CT‑C‑C‑‑‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                                                     >  1:2016312/7‑96 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGT‑TT‑A‑T‑‑‑TG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                                                   >  2:59395/7‑96 (MQ=25)
                 ggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                                                                                   >  2:534952/7‑100 (MQ=18)
                   aTGCGGCGTTAACTCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 <  1:2408202/100‑1 (MQ=255)
                    ttcggcgTTAACGACTTATCCGGCCTAC‑AAAATCGTGCAAATACAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAAGCCTGATAAGCGGAGCGCATCAGTCATTtttgc                                                                                >  1:2650609/3‑100 (MQ=18)
                    taaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑AT‑C‑CC‑‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg                                                                                 >  1:1830812/4‑95 (MQ=18)
                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                              <  1:2861154/100‑1 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTACAGA‑AT‑T‑A‑‑‑TC‑‑ATAAACATAATAAGCGGAGCACATCAGGCATTtttgcct                                                                              >  1:1370461/1‑100 (MQ=18)
                       ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                                                             <  2:1382847/100‑1 (MQ=255)
                       ggcaaacatcacttcgaaaacgtctctatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  1:3102815‑M1/74‑100 (MQ=255)
                       ggcaaacatcacttcgaaaacgtctctatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                      >  1:1333199‑M1/74‑100 (MQ=255)
                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAGTTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                                          <  1:1265379/100‑1 (MQ=255)
                           ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                         <  1:1200613/100‑1 (MQ=255)
                            acatcacttcgaaaacgtctctatcgaataattaaaagccatacaagactgatggcaaaggcaaaaatGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                 >  1:1786730‑M1/69‑100 (MQ=255)
                             aaCGCCTTATCCGACCTAC‑AAAATCGTGCAAATTCAATTCATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                                       >  2:2408127/1‑100 (MQ=255)
                               cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATAAGGCATTTTTGCCTTTGCCATCa                                                                     >  1:1015696/1‑100 (MQ=255)
                                gCCTTATCCGGTCTAC‑AAAATCGTGCAAATTCAATACATTACAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGTGCGTATCAGGCATTTTTGCCTTTGCCATCAg                                                                    <  2:915478/100‑1 (MQ=21)
                                 ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                                                   <  2:1455428/100‑1 (MQ=255)
                                     aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                               >  1:2584491/1‑100 (MQ=255)
                                     aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                               <  2:3762462/100‑1 (MQ=255)
                                     aTCCGGCCTAA‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                               <  1:802008/100‑1 (MQ=255)
                                      tCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                                              >  1:622140/1‑100 (MQ=255)
                                        cGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                                                            >  1:1044028/1‑100 (MQ=255)
                                         ggCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                           <  1:564194/100‑1 (MQ=255)
                                           ccTAC‑AAAATCGTGCAAATTCTATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                         >  2:3190842/1‑100 (MQ=255)
                                           ccTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                         >  1:2723415/1‑100 (MQ=255)
                                              aC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                                                       <  2:1205512/99‑1 (MQ=255)
                                                  aaaTCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                                                   <  1:3313972/100‑1 (MQ=255)
                                                   aaTCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                  <  1:887696/100‑1 (MQ=255)
                                                   aaTCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                                                  <  2:2324906/100‑1 (MQ=255)
                                                       gTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                              >  2:2149316/1‑100 (MQ=255)
                                                        tGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:3714728/100‑1 (MQ=255)
                                                        tGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:3215207/100‑1 (MQ=255)
                                                        tGCAAATTCAATACACTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                             <  1:3868206/100‑1 (MQ=255)
                                                         gCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                            <  1:3402754/100‑1 (MQ=255)
                                                          cAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatagatcgg                                           >  2:2603597/1‑93 (MQ=255)
                                                          cAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                           <  2:1989175/100‑1 (MQ=255)
                                                           aaaTTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATAAGTCTTGTATGGCTTTTAATTATTCGAt                                          >  1:484995/1‑100 (MQ=255)
                                                                 aaTACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc                                    >  1:775810/1‑100 (MQ=255)
                                                                       ttGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                              >  2:1691914/1‑100 (MQ=255)
                                                                         gCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            >  2:1556844/1‑100 (MQ=255)
                                                                         gCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            <  2:775801/100‑1 (MQ=255)
                                                                         gCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            >  1:1472838/1‑100 (MQ=255)
                                                                         gCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            >  1:2394447/1‑100 (MQ=255)
                                                                         gCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa                            <  2:3025539/100‑1 (MQ=255)
                                                                          cAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCATGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGaa                           >  1:1457587/1‑100 (MQ=255)
                                                                           aGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  1:668154/100‑1 (MQ=255)
                                                                           aGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                          <  1:3616515/100‑1 (MQ=255)
                                                                             a‑aT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg                        <  2:3462058/100‑1 (MQ=255)
                                                                                t‑t‑t‑‑‑tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt                      >  1:591142/1‑100 (MQ=255)
                                                                                t‑t‑t‑‑‑tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt                      >  1:279712/1‑100 (MQ=255)
                                                                                  t‑t‑‑‑tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATg                     <  2:3312742/100‑1 (MQ=255)
                                                                                    t‑‑‑tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGt                    <  2:2194126/100‑1 (MQ=255)
                                                                                            gTAGGCTTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 <  1:2662858/100‑1 (MQ=255)
                                                                                            gTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTg                 <  1:2459902/100‑1 (MQ=255)
                                                                                              aGGCTTGATAAGCGGAGCGGATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGcc               <  2:3110103/100‑1 (MQ=255)
                                                                                               ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGAGAtcggaagag              >  2:799123/1‑92 (MQ=38)
                                                                                                  CTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           >  1:287892‑M1/1‑100 (MQ=255)
                                                                                                  CTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAAt           <  1:2294402‑M1/100‑1 (MQ=255)
                                                                                                   tGATAAGCGGATCGCATCAGGCATTTTTTTCTTTTTCTTCTCTCTTATATGACTTTTAATTATTCGATACAGACGTTTTCTCAGTGATGTTTGCCTAAt           >  2:1805848/1‑99 (MQ=255)
                                                                                                    GATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATgg         <  1:3045446‑M1/100‑1 (MQ=255)
                                                                                                      TAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt       >  1:39681‑M1/1‑100 (MQ=255)
                                                                                                       aaGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATc      <  1:1691929/100‑1 (MQ=255)
                                                                                                        aGCGGAGCGCATCAGGCATTTTTGTCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     <  1:3584991/100‑1 (MQ=255)
                                                                                                        aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATTTTTGCCTAATGGATCa     >  2:2489706/1‑100 (MQ=255)
                                                                                                        aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa     >  2:3583350/1‑100 (MQ=255)
                                                                                                            |                                                                                                  
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑T‑T‑‑‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACC  >  NC_000913/3705949‑3706146

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