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.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACG  >  NC_000913/3705918‑3706113
                                                                                                    |                                                                                                       
gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGACTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                                                         >  2:180282/1‑100 (MQ=255)
 cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACa                                                                                                        >  1:479657/1‑100 (MQ=255)
   catccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                                                      <  2:3993677/97‑1 (MQ=255)
       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                                                                >  1:2885079/1‑100 (MQ=255)
         ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                                                             >  1:988076/1‑100 (MQ=255)
         ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                                                             <  2:1346090/100‑1 (MQ=255)
         ggTATACAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                                                             >  1:2759421/1‑100 (MQ=255)
          gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                                                   <  2:3732767/95‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCg                                                                                   <  2:4197409/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCg                                                                                   <  1:923343/100‑1 (MQ=255)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTg                                                                           >  2:4065163/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGa                                                                          >  1:3346544/1‑100 (MQ=255)
                         aCAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGAt                                                                         >  2:836076/1‑100 (MQ=255)
                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATa                                                                        >  2:308829/1‑100 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGc                                                                     <  2:219682/100‑1 (MQ=255)
                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAg                                                                 <  1:1347262/100‑1 (MQ=255)
                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAgc                                                                <  2:3966626/100‑1 (MQ=255)
                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAg                                                                 <  1:1042068/99‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGACTT‑‑‑TTCGTAGGCCTGATAAGCGGAgc                                                                >  2:1645393/1‑99 (MQ=255)
                                      tttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CACGTAGGCCTGATAAGCGTAGCGCAt                                                            >  1:2389468/8‑100 (MQ=17)
                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCAt                                                            >  2:1606354/1‑100 (MQ=255)
                                       ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCa                                                          >  1:3675378/7‑100 (MQ=18)
                                        tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCa                                                          >  2:2963406/6‑100 (MQ=17)
                                        tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑‑GACGTAGGCCTGATAAGCGTAGCGCATCAg                                                         >  1:1102970/6‑100 (MQ=18)
                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAg                                                         <  1:3412513/100‑1 (MQ=255)
                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                         >  1:252764/5‑100 (MQ=18)
                                         gaaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAATCTT‑‑‑CTGGTAGGCCTGATAAGCATCGCGCATCAg                                                         >  1:2830811/5‑100 (MQ=18)
                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAgg                                                        <  2:2627754/100‑1 (MQ=255)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                       >  2:1979436/3‑100 (MQ=18)
                                           aTGACGGATAAGGCGTTCACGCCGCATCCGTCCTACAAGTTCATGCAATTTTAATAAGTTGCCAGATAT‑‑‑CCTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                        >  2:1527434/1‑100 (MQ=25)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                      <  2:3006328/100‑1 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                      >  2:2714127/1‑100 (MQ=255)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                      >  1:3237597/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGct                                                      >  1:3112745/2‑99 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                        >  2:776051/2‑98 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                      >  1:1502194/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                      >  1:4050688/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTC‑‑‑CGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                      >  2:1888515/2‑100 (MQ=18)
                                            gcccggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATC‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                      >  1:3893438/11‑100 (MQ=17)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAggaa                                                      >  2:2838325/1‑97 (MQ=255)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                     >  1:2650065/1‑100 (MQ=255)
                                             gCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAAT‑‑‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                     >  1:1816695/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAAT‑‑‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                     >  2:4136408/1‑99 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                     >  2:326633/1‑99 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                     >  2:3868412/1‑98 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                     >  2:1423190/1‑98 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                      >  1:2267381/1‑100 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑‑GACGTAGGCCTGATAAGCGTAGCGCATCAGGcgat                                                    >  1:1712766/1‑97 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                     >  2:1309585/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                     >  1:2177069/1‑99 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAA‑‑‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGCag                                                     >  2:3626007/1‑99 (MQ=18)
                                              ccggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                    >  1:2446499/9‑98 (MQ=18)
                                              ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                    >  1:1072803/9‑98 (MQ=17)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                    >  2:254173/9‑98 (MQ=18)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                   >  2:2163276/9‑100 (MQ=21)
                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                    >  2:1190945/1‑100 (MQ=255)
                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                    >  1:4144737/1‑100 (MQ=255)
                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                    <  2:2502378/100‑1 (MQ=255)
                                               cggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATC‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                   >  2:2109531/8‑97 (MQ=17)
                                               cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                   >  1:396157/1‑100 (MQ=255)
                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  2:2346684/7‑100 (MQ=18)
                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  1:3278712/7‑100 (MQ=17)
                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                                                   >  1:3842269/7‑97 (MQ=21)
                                                ggataaggcgTTCACGCCTTATCCTGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGAT‑‑‑GGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttt                                                  >  1:3540081/7‑100 (MQ=14)
                                                ggataaggcgTTCACGCCTTATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  1:621797/7‑100 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  2:3599403/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  1:860604/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  2:1358885/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTCGTAGGACTGATAAGCGTAGCGCATCAGGCAattc                                                  >  1:771647/7‑99 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                  >  1:1468528/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                  >  1:1899437/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                 >  2:935024/7‑98 (MQ=21)
                                                 gataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg                                                >  1:203771/6‑97 (MQ=21)
                                                 gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                 <  2:3679299/100‑1 (MQ=255)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                <  1:1589448/100‑1 (MQ=255)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTTGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                               >  1:3771128/1‑100 (MQ=255)
                                                         gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                         <  2:3476437/100‑1 (MQ=255)
                                                         gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                         <  1:3591913/100‑1 (MQ=255)
                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                        <  1:1620320/100‑1 (MQ=255)
                                                           tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                       >  1:1182755/1‑100 (MQ=255)
                                                             aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                     <  2:1706227/100‑1 (MQ=255)
                                                               gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                   <  1:3713940/100‑1 (MQ=255)
                                                               gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                   >  2:2655945/1‑100 (MQ=255)
                                                                ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                  >  1:2167757/1‑100 (MQ=255)
                                                                 cTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                 >  1:131155/1‑100 (MQ=255)
                                                                  ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                <  1:738607/100‑1 (MQ=255)
                                                                    aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                              <  1:2629059/100‑1 (MQ=255)
                                                                      ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                            <  1:672771/100‑1 (MQ=255)
                                                                         gCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                         >  1:3936350/1‑100 (MQ=255)
                                                                          ccTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                        >  2:807629/1‑100 (MQ=255)
                                                                           cTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                       >  2:3146875/1‑100 (MQ=255)
                                                                            tACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                      >  2:1068362/1‑100 (MQ=255)
                                                                             aCAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCttt                     <  2:3782099/100‑1 (MQ=255)
                                                                              cAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTTAGTCTTGTATGGCtttt                    >  2:3354978/1‑100 (MQ=18)
                                                                                  aaTCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                 <  1:2605217/100‑1 (MQ=255)
                                                                                      gTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt             <  2:549819/100‑1 (MQ=255)
                                                                                      gTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt             >  1:496008/1‑100 (MQ=255)
                                                                                        ccgatctCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg           <  1:4070815/94‑1 (MQ=255)
                                                                                             ttCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga      <  1:2967100/100‑1 (MQ=255)
                                                                                                aaTACATTGCAGAATT‑‑‑TCCGTAGGCCTGATAAGCGGAGCGCATAAGGCATTTTTGACTTTGCCATCAGTCACGTATGGCTTTTACTTAGTCGTTAGAGAc   >  2:1139590/1‑100 (MQ=25)
                                                                                                    |                                                                                                       
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑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: