Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,2010TG52.9% 5.3 / 19.3 17T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (2/6);  new base G (9/0);  total (11/6)
Fisher's exact test for biased strand distribution p-value = 2.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-02

ATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCG  >  NC_000913/1356061‑1356349
                                                                                                                                            |                                                                                                                                                      
aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTG‑TTAAATGcc                                                                                                                                             <  1:191402/151‑1 (MQ=255)
   ggCCCACGAGGCGCGCGGGAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTATGCGGGG‑TTTAATACCCGa                                                                                                                                          >  2:532040/1‑151 (MQ=255)
                  cgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGG‑GTTAATACCAGAAAAgcggcggttat                                                                                                                            >  1:252058/1‑145 (MQ=255)
                    gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGG‑GTTAAAGCCAGATAGGGGGCGGTCATcc                                                                                                                          >  2:13173/1‑150 (MQ=255)
                       aaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGTG‑TTTAATGCCCCATAGGCGGCGTTTCTcccccg                                                                                                                      >  1:170873/1‑148 (MQ=255)
                         aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGAc                                                                                                                    <  1:424623/151‑1 (MQ=255)
                                   ccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTG‑GTAAATTCCCGATAGGCGGCGTTCATCCCCCGACGCAGCGttta                                                                                                          >  1:518034/1‑149 (MQ=255)
                                               tACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGG‑TTTAAAGCCAGAACAGCGGCGTTTATCCCCGGACGCAGCCGTTAACGCaaaagccg                                                                                              >  1:592747/1‑145 (MQ=255)
                                                 ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGG‑GTTAAAGCCAGAAAGGGGGCGTTCATCCCAGGACGCCGCGGTTAAAGCaaaggggggg                                                                                            >  2:339045/1‑143 (MQ=255)
                                                                             ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGG‑TTAAATGCCAGATAGGGGGCGGTCCTCCCCGGGCGCAGGGTTTAACGCCAGCGCGGGGGGGCCCG‑CAAAAAGGGTAgatgggctg                                                                >  2:440487/1‑146 (MQ=255)
                                                                             ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGG‑TTAAATGCCAGATAGGCGACGGTTATCCCCGGAAGCAGCGTTTAAAGCAAAGGCGGGGGGGCCCC‑CAAAAAGGAAAGctggctggg                                                               >  2:370761/1‑150 (MQ=255)
                                                                                                      aTGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGG‑TTAAAAGGCAGAAAAGCGAGGTTCATCCCAGGACGCAGCGGTAAAAACAAACGCGGGGGGCCACG‑CAAAAAGGATAGCTGGGGGGGAGCAGGGCCGGGCAGAACGTCGCa                                       >  1:172303/1‑150 (MQ=255)
                                                                                                                gCGGGATTATTTAGCGGCGGGTTGGGGGGGTTTAAATCCCAAAAAGGGGCGTTTCTCCCCAGGAGCCGGGGTTTAAAGC‑AGCGGGGGGGGTCCCGCCAAAATGGAAAGCGGGCGGGCAGCAGGCCAGGCACGAACGACCCATTCCCCggc                              >  1:318080/1‑149 (MQ=255)
                                                                                                                     attattTAGCGGCGGTTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAg                        <  1:511697/150‑1 (MQ=255)
                                                                                                                                      tGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGt      <  1:531876/151‑1 (MQ=255)
                                                                                                                                        gCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTg     <  1:81163/150‑1 (MQ=255)
                                                                                                                                         cGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgc    <  2:151091/150‑1 (MQ=255)
                                                                                                                                           gtg‑tTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccg  >  1:395443/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                      
ATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTG‑TTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCG  >  NC_000913/1356061‑1356349

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