Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,640,9810CT44.4% 33.4 / 17.2 36F40F (TTC→TTT) dtpBdipeptide and tripeptide permease B
Reads supporting (aligned to +/- strand):  ref base C (3/17);  new base T (16/0);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.35e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTT  >  NC_000913/3640836‑3641119
                                                                                                                                                  |                                                                                                                                          
aGGATGGGATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTttatttct                                                                                                                                       >  1:238633/1‑148 (MQ=255)
      gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCTTTtaa                                                                                                                                   >  2:431345/1‑149 (MQ=255)
          tGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCCGTTaacccct                                                                                                                               >  1:391832/1‑145 (MQ=255)
               ggTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTGCGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGAta                                                                                                                        <  1:127998/151‑2 (MQ=255)
               ggTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCAGTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGAtt                                                                                                                        <  1:164102/151‑1 (MQ=255)
                   aataaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAAAAACTTTGATtct                                                                                                                      >  2:373666/1‑149 (MQ=255)
                      aaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCa                                                                                                                  <  2:699114/150‑1 (MQ=255)
                       aataTGACTACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTaaaaaaactgggttt                                                                                                                       >  1:296220/1‑132 (MQ=255)
                         tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGGTTTTTTGTTAAACactctgggttct                                                                                                                      >  1:200692/1‑132 (MQ=255)
                                           cATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACtttt                                                                                            <  2:604009/151‑1 (MQ=255)
                                               gggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTg                                                                                        <  2:260408/151‑1 (MQ=255)
                                                 gATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGGTAAACAACTTGGGTTCTCCCAAGAGGAGGCTTTTTTCACtttttg                                                                                          >  1:175440/1‑145 (MQ=255)
                                                  aTGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCtt                                                                                     <  1:286354/151‑1 (MQ=255)
                                                     ctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCTTTTAAAAGCTTTGGATCTCGCAAGAGCAGGCTTTTTTCACTTTTTGTGGttttt                                                                                  >  2:289101/1‑150 (MQ=255)
                                                      tgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTTAAAAGCCTGGATTCTTGCCAGAGGCGGCTTTTTTCACTTTTTGGGGtttt                                                                                   >  1:288142/1‑149 (MQ=255)
                                                       gcagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTtgc                                                                                 <  1:447222/150‑1 (MQ=255)
                                                          gcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTTAACAACCTTGGTTCTCGCGAGAGCAGGCTTTTTTTACTTTTTGTGCttttttt                                                                                >  1:47358/1‑145 (MQ=255)
                                                           caACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTgcgc                                                                            <  2:472457/151‑1 (MQ=255)
                                                            aaCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTTTTTTAAAAGCTTTGGTTTTTGCCAAAGCAGGGTTTTTTTTCTTTTTGGGGttttttt                                                                                >  1:240606/1‑143 (MQ=255)
                                                                   cccATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTAt                                                                    <  2:290799/151‑1 (MQ=255)
                                                                         tttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAAAAGCTTGGGTTTTCGCCAGAGCAGGGTTTTTTCACTTTTTGTGGTTTTTTTTCGCTGGGCTATgggccc                                                              >  2:638141/1‑147 (MQ=255)
                                                                                     ttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAACAACTTTGATTTTCGCAAGAGCAGGCTTTTTTCACTTTTTGTGCTTTTTCTGCGCTGGGCTATGGCCTCCttttcctt                                                     >  2:695997/1‑143 (MQ=255)
                                                                                     ttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTAAACCGCTTGGGTTTTCGCAAGAGCAGGCTTTTTTTCCTTTTGGTGGTTTTTTTGTGCCGGTCTTTGGCCTCCttttccct                                                     >  1:378141/1‑143 (MQ=255)
                                                                                            gAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTAAAAAGCCTGGGTTTTCCCAAGAGCAGGCCTTTTTCTCTTTTTGGGGTTTTGCTGCGCTGGGTTATGGCCTCATTTCCAtttgggg                                                >  2:411833/1‑141 (MQ=255)
                                                                                             aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTAAAAAGCTTGGATTTTCGCAAGAGGAGGGTTTTTTCACTTTTGGTGCTTTTTCTGCGCTGGGCTATGGCCCCATTTTCATTTGGGGGTATGt                                          >  2:258245/1‑151 (MQ=255)
                                                                                             aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTCTTAAAAAGCTTGGATTATCGCCAGAGCAGGCTTTTTTTCCCTTTTGTGCTTTTTCTGCGGTTGGCTTTGGCCTCCTTTTCCTTtgcgg                                                >  1:420656/1‑145 (MQ=255)
                                                                                             aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGt                                          <  1:143992/151‑1 (MQ=255)
                                                                                                tATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTCTTTAACCCCATGGATTCTACCCAAAGGCGGCTTTTTGGCCTTTTTGGGCTTTTTTTGCGGCGGCTTATGGGCCCCTTTCCCTTggcgg                                                >  2:579616/1‑142 (MQ=255)
                                                                                                   gggAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCg                                     <  1:431887/150‑1 (MQ=255)
                                                                                                                  actacGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTgggg                       <  1:713379/149‑1 (MQ=255)
                                                                                                                    tacGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAACAGCCTTGATTTTCGCAAGAGCAGGCTTTTTTTACTTTTTGTGGTTTTTCTGCGCTGGTCTTTTGCCTCATTTTCCTTTGGGGGTATGTTGGCGgccccccccttggggac                    >  2:535607/1‑139 (MQ=255)
                                                                                                                     acGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCaa                  <  2:593432/151‑1 (MQ=255)
                                                                                                                     acGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCa                   <  2:378360/150‑1 (MQ=255)
                                                                                                                       ggCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTTACTTTTTGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTTGCGGCTATGTCGGCGACCCCCCGCTGGGGACCaaa                 >  2:79106/1‑150 (MQ=255)
                                                                                                                            gCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCAcc           <  2:438998/151‑1 (MQ=255)
                                                                                                                                 gCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGt      <  2:93425/151‑1 (MQ=255)
                                                                                                                                     aCTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCATGAGCAGGCTTGTTTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCtt  <  1:535825/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
AGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTT  >  NC_000913/3640836‑3641119

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑