Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,028 T→C 23.3% intergenic (‑238/+70) 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,0280TC23.3% 609.4 / 192.9 324intergenic (‑238/+70)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (139/103);  new base C (71/3);  total (217/107)
Fisher's exact test for biased strand distribution p-value = 1.53e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 326 total aligned reads displayed.

TACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTG  >  NC_000913/3705931‑3706124
                                                                                                    |                                                                                                      
taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAt‑t‑‑‑t‑t                                                                                                  <  1:2567488/100‑1 (MQ=255)
  cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑g                                                                                              <  2:799482/100‑1 (MQ=255)
   aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑Gt                                                                                             >  1:1320482/1‑100 (MQ=255)
    gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTa                                                                                            >  1:3090983/1‑100 (MQ=255)
     aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAg                                                                                           <  1:2286545/100‑1 (MQ=255)
       tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGc                                                                                         <  1:860543/100‑1 (MQ=255)
        cttccgatcttGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGcc                                                                                        <  2:1856683/90‑1 (MQ=255)
           gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGa                                                                                     <  1:1805404/100‑1 (MQ=255)
                aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGc                                                                                >  2:3349954/1‑100 (MQ=255)
                  tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCgg                                                                              <  1:2174010/100‑1 (MQ=255)
                     gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAg                                                                            <  1:1042068/99‑1 (MQ=255)
                      cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GACT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAgc                                                                           >  2:1645393/1‑99 (MQ=255)
                        ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                        <  2:36542/100‑1 (MQ=255)
                         tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                       >  1:2618801/8‑100 (MQ=17)
                           tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                     <  2:1601408/100‑1 (MQ=255)
                           tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑ATA‑T‑‑‑G‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                    >  1:1102970/6‑100 (MQ=18)
                            gaaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑ATCT‑T‑‑‑C‑TG‑‑GTAGGCCTGATAAGCATCGCGCATCAg                                                                    >  1:2830811/5‑100 (MQ=18)
                             taggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                   >  2:2482918/4‑100 (MQ=17)
                             tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                    <  2:1484001/99‑1 (MQ=255)
                              aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                  >  2:1979436/3‑100 (MQ=18)
                              aggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAAGCGTGCAAATTCAATACATTGCATGGGC‑C‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                  >  1:3366364/3‑100 (MQ=17)
                              aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑C‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                  >  1:1359021/3‑100 (MQ=21)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCC‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                >  1:3628725/2‑100 (MQ=18)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAAT‑C‑‑‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                >  1:864509/2‑100 (MQ=18)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 >  2:1747157/2‑100 (MQ=18)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 >  2:3765236/2‑100 (MQ=17)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 >  1:75441/2‑100 (MQ=17)
                               ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑AAT‑TCAAC‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                  >  1:3585531/2‑100 (MQ=21)
                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CGAT‑C‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 >  2:3096200/2‑100 (MQ=21)
                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑C‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 >  1:472568/2‑100 (MQ=21)
                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑C‑‑‑T‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                 >  2:3855021/2‑100 (MQ=18)
                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑ATA‑T‑‑‑G‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga                                                                >  1:3519790/2‑98 (MQ=21)
                                gCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCA‑GAAA‑T‑‑‑C‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                >  2:4136408/1‑99 (MQ=18)
                                gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGCA‑‑ATTCT‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                                >  2:1423190/1‑98 (MQ=18)
                                gCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                                  >  2:3690750/1‑98 (MQ=17)
                                 cctgataaggcgTTTACGCCTTATCAGGCCTAC‑AAGATCGTGCAAATTCAACATATTGC‑‑CACT‑C‑‑‑A‑CCCAGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                                                                >  2:2651227/9‑99 (MQ=21)
                                 ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGGGC‑C‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                               >  1:1072803/9‑98 (MQ=17)
                                 ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TGAG‑C‑‑‑C‑AC‑‑GTAGGCCTGATAAGCATAACGCATCAGGCa                                                                 >  1:2668686/9‑98 (MQ=17)
                                 ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑C‑‑‑C‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                               >  2:914270/9‑98 (MQ=21)
                                 ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AGAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                                                               >  1:4028903/9‑98 (MQ=17)
                                 ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                               >  1:4144737/1‑100 (MQ=255)
                                 ccGGATAAGGCGGTCACGCCGCATCCGGCCTAC‑GAAACCGCAAGAATTCAATATATTGCA‑GGAG‑C‑‑‑G‑GT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                               >  2:3085/1‑98 (MQ=18)
                                  cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                              <  2:3259963/100‑1 (MQ=255)
                                   ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCC‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                            >  1:2464264/7‑98 (MQ=21)
                                   ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑C‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  2:1209954/7‑100 (MQ=18)
                                   ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑AAT‑TCAAC‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                                                              >  1:3842269/7‑97 (MQ=21)
                                   ggataaggcgTTCACGCCTTATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CGAT‑C‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  1:621797/7‑100 (MQ=21)
                                   ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑C‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  1:3784957/7‑100 (MQ=18)
                                   ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑C‑‑‑A‑T‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  1:2400924/7‑100 (MQ=18)
                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑C‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  1:2149184/7‑100 (MQ=21)
                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑C‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  2:1358885/7‑100 (MQ=18)
                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAG‑C‑‑‑T‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                             >  1:1468528/7‑100 (MQ=17)
                                   ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGA‑T‑‑‑T‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                             >  1:1899437/7‑100 (MQ=17)
                                    gataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCC‑‑‑C‑TT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg                                                           >  1:203771/6‑97 (MQ=21)
                                    gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑T‑‑‑C‑AC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                            <  2:2100866/95‑1 (MQ=17)
                                    gataaggcgTTCACGCCGCATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTGCA‑GAGA‑T‑‑‑G‑GC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttttt                                                            >  1:3250564/6‑100 (MQ=14)
                                     aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                           <  2:396156/100‑1 (MQ=255)
                                            gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCATATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                    >  2:2517834/1‑100 (MQ=255)
                                            gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                                    <  1:1497072/100‑1 (MQ=255)
                                             ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATCCATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                   <  2:2698226/100‑1 (MQ=255)
                                              tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                                  >  1:1182755/1‑100 (MQ=255)
                                               aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt                                                 >  2:3610401/1‑100 (MQ=255)
                                                aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                <  2:1332259/100‑1 (MQ=255)
                                                aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                                <  2:1706227/100‑1 (MQ=255)
                                                       aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                         <  2:3005166/100‑1 (MQ=255)
                                                        tCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                        >  1:710336/1‑100 (MQ=255)
                                                         ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                                       <  1:1698568/100‑1 (MQ=255)
                                                           ggCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                     <  2:2185394/100‑1 (MQ=255)
                                                           ggCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                     <  2:257182/100‑1 (MQ=255)
                                                            gCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                                    >  1:3936350/1‑100 (MQ=255)
                                                                aC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                <  1:2935653/100‑1 (MQ=255)
                                                                aC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                                <  1:3149966/100‑1 (MQ=255)
                                                                   aaaaTCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                              <  2:4046120/100‑1 (MQ=255)
                                                                        cGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                         <  1:3955061/100‑1 (MQ=255)
                                                                         gTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                        <  2:549819/100‑1 (MQ=255)
                                                                         gTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                        >  1:496008/1‑100 (MQ=255)
                                                                         gTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                        >  2:1069108/1‑100 (MQ=255)
                                                                          tGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                       <  1:940906/100‑1 (MQ=255)
                                                                          tGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                       >  2:3494064/1‑100 (MQ=255)
                                                                          tGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                       <  1:2629258/100‑1 (MQ=255)
                                                                          tGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                       <  1:1606358/100‑1 (MQ=255)
                                                                           gCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                      <  2:3975158/100‑1 (MQ=255)
                                                                           ccgatctCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                      <  1:4070815/94‑1 (MQ=255)
                                                                            cAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATGTTTGCCTTTGCAATCAGTCTTGTATGGCTTTTAattatt                        >  2:2951599/1‑97 (MQ=255)
                                                                              aaTTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa                   <  1:660167/100‑1 (MQ=255)
                                                                                ttCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                 >  2:829143/1‑100 (MQ=255)
                                                                                ttCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga                 <  1:2967100/100‑1 (MQ=255)
                                                                                   aaTACATTGCA‑GAAT‑T‑‑‑T‑CC‑‑GTAGGCCTGATAAGCGGAGCGCATAAGGCATTTTTGACTTTGCCATCAGTCACGTATGGCTTTTACTTAGTCGTTAGAGAc              >  2:1139590/1‑100 (MQ=25)
                                                                                    aTACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg             >  2:215955/1‑100 (MQ=255)
                                                                                         ttGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc        <  1:3683144/100‑1 (MQ=255)
                                                                                         ttGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc        >  1:3895367/1‑100 (MQ=255)
                                                                                         ttGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc        >  2:223491/1‑100 (MQ=255)
                                                                                          tGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGtttt         <  2:3441592/98‑1 (MQ=255)
                                                                                           gCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa      >  1:3726310/1‑100 (MQ=255)
                                                                                           gCA‑GAAT‑G‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTGGCCATCAGTCTTGTATGGCGTTTAATTATTCGATAGAGACGTTTTCGa      <  2:2771372/100‑1 (MQ=255)
                                                                                           gCA‑GAAT‑C‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGa      <  1:1580369/100‑1 (MQ=255)
                                                                                               gAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt   <  2:4184273/100‑1 (MQ=255)
                                                                                               gAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt   >  1:3326233/1‑100 (MQ=255)
                                                                                                aaT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg  >  1:2241978/1‑100 (MQ=255)
                                                                                                 aT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCCGAGCACATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTg  <  2:505188/99‑1 (MQ=255)
                                                                                                    |                                                                                                      
TACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑T‑‑‑T‑TC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTG  >  NC_000913/3705931‑3706124

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