Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,706,0170CT28.3% 423.6 / 169.4 247intergenic (‑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 (80/82);  new base T (63/1);  total (163/84)
Fisher's exact test for biased strand distribution p-value = 2.25e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 251 total aligned reads displayed.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116
                                                                                                   |                                                                                                       
  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                                                                       >  2:312023/1‑100 (MQ=255)
  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                                                                       <  2:1494721/100‑1 (MQ=255)
   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                                                      >  2:181037/1‑100 (MQ=255)
    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                                                                     >  1:2550098/1‑100 (MQ=255)
    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                                                                     >  2:1103448/1‑100 (MQ=255)
    cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                                                                     >  2:1264582/1‑100 (MQ=255)
     aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                                                    >  1:1597698/1‑100 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCg                                                                                      <  1:23573/100‑1 (MQ=255)
                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGt                                                                                     <  2:2372065/100‑1 (MQ=255)
                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTa                                                                                    >  1:1043269/1‑100 (MQ=255)
                    tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGc                                                                                 <  1:2242478/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGcc                                                                                <  2:718562/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGcc                                                                                <  2:2773431/100‑1 (MQ=255)
                      ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCt                                                                               >  2:1637616/1‑100 (MQ=255)
                      taaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCt                                                                               >  1:1291537/4‑100 (MQ=18)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTg                                                                              >  1:3091987/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGa                                                                             >  1:58313/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGa                                                                             >  1:498859/1‑100 (MQ=255)
                         aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGAt                                                                            >  2:3317488/1‑100 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAgcg                                                                  <  1:1544037/100‑1 (MQ=255)
                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCa                                                                <  1:2039716/100‑1 (MQ=255)
                                       ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCa                                                             >  1:2400999/7‑100 (MQ=21)
                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATc                                                              >  2:1169837/1‑100 (MQ=255)
                                        tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCa                                                             >  2:2585220/6‑100 (MQ=17)
                                        tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATT‑‑‑TTCGTAGGCCTGATAAGCGTAGCGCATCa                                                             >  2:529276/6‑100 (MQ=17)
                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATT‑‑‑TTCGTAGGCCTGATAAGCGTAGCGCATCAg                                                            >  2:1654854/5‑100 (MQ=17)
                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAGAA‑‑‑AATGTAGGCCTGATAAGCGTAGCGCATCAgg                                                           >  1:1821335/5‑100 (MQ=18)
                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAgg                                                           >  2:1840102/4‑100 (MQ=17)
                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAgg                                                           >  1:3230938/4‑100 (MQ=21)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:1000104/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  2:602221/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑‑GACGTAGGCCTGATAAGCGTAGCGCATCAGGcg                                                         >  1:1298502/3‑99 (MQ=21)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                          >  1:883459/3‑100 (MQ=21)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  1:1370537/1‑100 (MQ=255)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:3166654/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAATCTT‑‑‑CTGGTAGGCCTGATAAGCATCGCGCATCAGGCa                                                         >  1:2145711/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:2287294/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:1589755/2‑100 (MQ=18)
                                             gctggataaggcgTTCACGCCGCATCCAGCCTACAAAATTGTGCAAAATCAATGGATTGCACAGCT‑‑‑AGCGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                        >  2:1240779/10‑99 (MQ=25)
                                             gCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAAT‑‑‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                        >  2:1937770/1‑99 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                       >  1:3315229/1‑100 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                        >  2:1747760/1‑99 (MQ=18)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                       >  1:2497360/1‑100 (MQ=21)
                                              ccggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                       >  2:159578/9‑98 (MQ=18)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaat                                                       >  2:412421/9‑98 (MQ=25)
                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                                                    >  1:1130178/7‑98 (MQ=21)
                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:3106379/7‑100 (MQ=18)
                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:2591136/7‑100 (MQ=18)
                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:3040246/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  1:411605/7‑100 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                     >  1:1205073/7‑100 (MQ=17)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                                                     >  1:2861830/7‑100 (MQ=17)
                                                 gataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGAT‑‑‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattta                                                    >  2:975239/6‑99 (MQ=17)
                                                 gataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCATTtatg                                                   >  1:1223879/6‑100 (MQ=21)
                                                  ataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCATTtatgc                                                  <  2:2497348/96‑1 (MQ=25)
                                                    ttccgatctaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                 <  1:2207786/91‑1 (MQ=255)
                                                    gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTGGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                 <  2:3345524/100‑1 (MQ=255)
                                                    gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                 <  1:3067158/100‑1 (MQ=255)
                                                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                >  2:2550205/1‑100 (MQ=255)
                                                     aggcgTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑AATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATttggcg                                                 <  1:2588144/99‑2 (MQ=25)
                                                      ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt                                               <  1:2957401/100‑1 (MQ=255)
                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                           <  2:833073/100‑1 (MQ=255)
                                                             aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAtn                                        <  2:3294518/100‑2 (MQ=255)
                                                             aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                        <  1:532806/100‑1 (MQ=255)
                                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                       <  2:2751557/100‑1 (MQ=255)
                                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                       >  1:737464/1‑100 (MQ=255)
                                                                ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt                                     <  2:2944052/100‑1 (MQ=255)
                                                                 cTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc                                    <  1:380573/100‑1 (MQ=255)
                                                                  ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCGTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                   <  1:1938700/100‑1 (MQ=255)
                                                                  ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                   <  2:1043267/100‑1 (MQ=255)
                                                                  ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt                                   >  2:1407773/1‑100 (MQ=255)
                                                                    aTCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                 <  1:143647/100‑1 (MQ=255)
                                                                     tCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGt                                <  2:1511071/100‑1 (MQ=255)
                                                                      ccGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                               >  2:786280/1‑100 (MQ=255)
                                                                       cGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              <  1:2759549/100‑1 (MQ=255)
                                                                       cGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              >  1:1743919/1‑100 (MQ=255)
                                                                           cTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                          >  1:1984517/1‑100 (MQ=255)
                                                                            tACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                         >  1:3325570/1‑100 (MQ=255)
                                                                              cAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                       <  1:640530/100‑1 (MQ=255)
                                                                              cAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                       <  1:282578/100‑1 (MQ=255)
                                                                              cAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                       <  1:124775/100‑1 (MQ=255)
                                                                              cAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGAGAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                       <  1:1009354/100‑1 (MQ=255)
                                                                                aaaTCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTaa                     >  1:1149378/1‑100 (MQ=255)
                                                                                 aaTCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                    <  2:2494877/100‑1 (MQ=255)
                                                                                 aaTCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                    <  1:1577153/100‑1 (MQ=255)
                                                                                  aTCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatt                   <  2:1048064/100‑1 (MQ=255)
                                                                                    cGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                 >  2:192874/1‑100 (MQ=255)
                                                                                    cGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                 <  2:3352237/100‑1 (MQ=255)
                                                                                       gCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg              <  1:1189761/100‑1 (MQ=255)
                                                                                       gCAAATTAAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAGCAGTCTTGTATGGCTTTTAATTATTCg              <  1:1707013/100‑1 (MQ=255)
                                                                                         aaaTTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt            <  1:1524885/100‑1 (MQ=255)
                                                                                         aaaTTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt            >  2:217044/1‑100 (MQ=255)
                                                                                           aTTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag          >  2:2108375/1‑100 (MQ=255)
                                                                                             tCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATagag        >  1:758954/1‑100 (MQ=255)
                                                                                              cAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGGCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga       >  2:2257189/1‑100 (MQ=255)
                                                                                              cAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga       <  2:2654874/100‑1 (MQ=255)
                                                                                               aaTACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc      >  2:1354620/1‑100 (MQ=255)
                                                                                                aTACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg     <  1:801937/100‑1 (MQ=255)
                                                                                                   cATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt  >  2:1560597/1‑100 (MQ=255)
                                                                                                   cATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGttt  <  1:1798122/100‑1 (MQ=255)
                                                                                                   |                                                                                                       
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116

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