Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,017 C→T 32.8% 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,0170CT32.8% 454.6 / 249.0 314intergenic (‑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 (92/101);  new base T (90/3);  total (210/104)
Fisher's exact test for biased strand distribution p-value = 1.01e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

Only 100 of 315 total aligned reads displayed.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACG  >  NC_000913/3705918‑3706113
                                                                                                    |                                                                                                       
gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGACTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                                                         >  2:180282/1‑100 (MQ=255)
 cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACa                                                                                                        >  1:479657/1‑100 (MQ=255)
   catccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                                                      <  2:3993677/97‑1 (MQ=255)
       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑g                                                                                                >  1:2885079/1‑100 (MQ=255)
         ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑Ga‑a                                                                                             >  1:988076/1‑100 (MQ=255)
         ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑Ga‑a                                                                                             <  2:1346090/100‑1 (MQ=255)
         ggTATACAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑Ga‑a                                                                                             >  1:2759421/1‑100 (MQ=255)
          gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑c                                                                                                   <  2:3732767/95‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCg                                                                                   <  2:4197409/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCg                                                                                   <  1:923343/100‑1 (MQ=255)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTg                                                                           >  2:4065163/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGa                                                                          >  1:3346544/1‑100 (MQ=255)
                         aCAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGAt                                                                         >  2:836076/1‑100 (MQ=255)
                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATa                                                                        >  2:308829/1‑100 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGc                                                                     <  2:219682/100‑1 (MQ=255)
                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAg                                                                 <  1:1347262/100‑1 (MQ=255)
                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAgc                                                                <  2:3966626/100‑1 (MQ=255)
                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAg                                                                 <  1:1042068/99‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑CT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAgc                                                                >  2:1645393/1‑99 (MQ=255)
                                      tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GA‑GT‑T‑‑‑CACGTAGGCCTGATAAGCGTAGCGCAt                                                            >  1:2389468/8‑100 (MQ=17)
                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCAt                                                            >  2:1606354/1‑100 (MQ=255)
                                       ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GA‑AT‑CC‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCa                                                          >  1:3675378/7‑100 (MQ=18)
                                        tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GA‑GT‑T‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCa                                                          >  2:2963406/6‑100 (MQ=17)
                                        tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGTTCGTGCAAATTCAATAAATTG‑CA‑‑A‑TA‑T‑‑‑GACGTAGGCCTGATAAGCGTAGCGCATCAg                                                         >  1:1102970/6‑100 (MQ=18)
                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAg                                                         <  1:3412513/100‑1 (MQ=255)
                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTG‑CA‑‑A‑TTCT‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                         >  1:252764/5‑100 (MQ=18)
                                         gaaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTG‑CA‑AT‑CT‑T‑‑‑CTGGTAGGCCTGATAAGCATCGCGCATCAg                                                         >  1:2830811/5‑100 (MQ=18)
                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAgg                                                        <  2:2627754/100‑1 (MQ=255)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                       >  2:1979436/3‑100 (MQ=18)
                                           aTGACGGATAAGGCGTTCACGCCGCATCCGTCCTAC‑AAGTTCATGCAATTTTAATAAGTTGCCA‑GA‑TA‑T‑‑‑CCTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                        >  2:1527434/1‑100 (MQ=25)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                      <  2:3006328/100‑1 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                      >  2:2714127/1‑100 (MQ=255)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                      >  1:3237597/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTG‑CA‑AT‑TC‑T‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGct                                                      >  1:3112745/2‑99 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCAAATTCAACATATTG‑CA‑AT‑TC‑T‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                        >  2:776051/2‑98 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTG‑TG‑GG‑CT‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                      >  1:1502194/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTG‑TG‑GG‑CT‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                      >  1:4050688/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTG‑CA‑AT‑CT‑C‑‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                      >  2:1888515/2‑100 (MQ=18)
                                            gcccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTG‑CA‑GAGAT‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                      >  1:3893438/11‑100 (MQ=17)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAggaa                                                      >  2:2838325/1‑97 (MQ=255)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                     >  1:2650065/1‑100 (MQ=255)
                                             gCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTG‑CA‑GA‑AA‑T‑‑‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                     >  1:1816695/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTG‑CA‑GA‑AA‑T‑‑‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                     >  2:4136408/1‑99 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                     >  2:326633/1‑99 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTG‑CA‑AT‑TC‑T‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                     >  2:3868412/1‑98 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTG‑CA‑‑A‑TTCT‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                     >  2:1423190/1‑98 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTG‑C‑‑‑A‑AT‑TCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                      >  1:2267381/1‑100 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTG‑CA‑‑A‑TA‑T‑‑‑GACGTAGGCCTGATAAGCGTAGCGCATCAGGcgat                                                    >  1:1712766/1‑97 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTG‑CA‑TG‑AG‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                     >  2:1309585/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CC‑AT‑C‑‑‑CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                     >  1:2177069/1‑99 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTG‑CA‑GG‑AA‑A‑‑‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGCag                                                     >  2:3626007/1‑99 (MQ=18)
                                              ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                    >  1:2446499/9‑98 (MQ=18)
                                              ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTG‑CATGG‑GC‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                    >  1:1072803/9‑98 (MQ=17)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTG‑TG‑GG‑CT‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                    >  2:254173/9‑98 (MQ=18)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTG‑CA‑GG‑AG‑T‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                   >  2:2163276/9‑100 (MQ=21)
                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                    >  2:1190945/1‑100 (MQ=255)
                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                    >  1:4144737/1‑100 (MQ=255)
                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                    <  2:2502378/100‑1 (MQ=255)
                                               cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTG‑CA‑GAGAT‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                   >  2:2109531/8‑97 (MQ=17)
                                               cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                   >  1:396157/1‑100 (MQ=255)
                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  2:2346684/7‑100 (MQ=18)
                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GA‑GT‑T‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  1:3278712/7‑100 (MQ=17)
                                                ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTG‑C‑‑‑A‑AT‑TCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                                                   >  1:3842269/7‑97 (MQ=21)
                                                ggataaggcgTTCACGCCTTATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTG‑CA‑GA‑GA‑T‑‑‑GGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttt                                                  >  1:3540081/7‑100 (MQ=14)
                                                ggataaggcgTTCACGCCTTATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTG‑AT‑CG‑AT‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  1:621797/7‑100 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑CA‑GG‑AC‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  2:3599403/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑CA‑GG‑AC‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  1:860604/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTG‑TG‑GG‑CT‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  2:1358885/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CC‑AT‑C‑‑‑CTCGTAGGACTGATAAGCGTAGCGCATCAGGCAattc                                                  >  1:771647/7‑99 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTG‑CA‑GG‑AG‑C‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  1:1468528/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTG‑CA‑GA‑GA‑T‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                  >  1:1899437/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTG‑CA‑GG‑AG‑T‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                 >  2:935024/7‑98 (MQ=21)
                                                 gataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GA‑ATCC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg                                                >  1:203771/6‑97 (MQ=21)
                                                 gATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                 <  2:3679299/100‑1 (MQ=255)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                <  1:1589448/100‑1 (MQ=255)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTTGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                               >  1:3771128/1‑100 (MQ=255)
                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                         <  2:3476437/100‑1 (MQ=255)
                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                         <  1:3591913/100‑1 (MQ=255)
                                                          ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                        <  1:1620320/100‑1 (MQ=255)
                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                       >  1:1182755/1‑100 (MQ=255)
                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                     <  2:1706227/100‑1 (MQ=255)
                                                               gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                   <  1:3713940/100‑1 (MQ=255)
                                                               gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                   >  2:2655945/1‑100 (MQ=255)
                                                                ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                  >  1:2167757/1‑100 (MQ=255)
                                                                 cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                 >  1:131155/1‑100 (MQ=255)
                                                                  ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                <  1:738607/100‑1 (MQ=255)
                                                                    aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                              <  1:2629059/100‑1 (MQ=255)
                                                                      ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                            <  1:672771/100‑1 (MQ=255)
                                                                         gCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                         >  1:3936350/1‑100 (MQ=255)
                                                                          ccTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                        >  2:807629/1‑100 (MQ=255)
                                                                           cTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                       >  2:3146875/1‑100 (MQ=255)
                                                                            tAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                      >  2:1068362/1‑100 (MQ=255)
                                                                             aC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                     <  2:3782099/100‑1 (MQ=255)
                                                                              c‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTTAGTCTTGTATGGCtttt                    >  2:3354978/1‑100 (MQ=18)
                                                                                  aaTCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                 <  1:2605217/100‑1 (MQ=255)
                                                                                      gTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt             <  2:549819/100‑1 (MQ=255)
                                                                                      gTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt             >  1:496008/1‑100 (MQ=255)
                                                                                        ccgatctCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg           <  1:4070815/94‑1 (MQ=255)
                                                                                             ttCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga      <  1:2967100/100‑1 (MQ=255)
                                                                                                aaTACATTG‑CA‑GA‑AT‑T‑‑‑TCCGTAGGCCTGATAAGCGGAGCGCATAAGGCATTTTTGACTTTGCCATCAGTCACGTATGGCTTTTACTTAGTCGTTAGAGAc   >  2:1139590/1‑100 (MQ=25)
                                                                                                    |                                                                                                       
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACG  >  NC_000913/3705918‑3706113

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