Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913989,9430TG57.1% 5.3 / 13.4 14S138S (TCA→TCC) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base G (4/4);  total (7/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAAC‑‑GGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGG‑CGGTTAACGCCGAGAAGTCGACT‑AATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTT‑CACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCA  >  NC_000913/989719‑990160
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                          
gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAAC‑‑GGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGG‑CGGTTAACGCCGAGAAGTCGACT‑AATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTT‑CACt                                                                                                                                                                                                 <  1:186571/251‑1 (MQ=255)
   aGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAATCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAAC‑‑GGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGG‑CGGTTAACGCCGAGAAGTCGACT‑AATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTT‑CACTGAt                                                                                                                                                                                              <  1:243764/251‑1 (MQ=255)
             gACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAAC‑‑GGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGG‑CGGTTTAAGCCGAGAAGTCGGCTTAATTT‑CTTTTCAACGTGTCGGGGGGTTGGGGAACCTCGGCCTGCGGGGAGCGATAGCGATGTAACATTT‑cccttgtaaaccccac                                                                                                                                                                                     >  1:221718/1‑235 (MQ=255)
                                   ggTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTCCTTCGCGAGAAAAAC‑‑GGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGG‑CGGTTAACACCGCGAAGTCGACT‑AATTTTCTTTTCAGCGTGTCGGGGGCTTGCGCAACGTCGGCCCTCGGGGAGGGAAAGCGATGGTCCATTT‑CACTGATAACACACAGCAAAGGAAATTCCCAGAGGATg                                                                                                                                                               >  1:342880/1‑250 (MQ=255)
                                    gggCCCAAAATCCAGGATTCCGGGTGCAGACGCTTAACAAGGTTTTTTTGTACTTGGAAAGAAAAAC‑‑GGCGACGGGTGCTAAAACCCATCAGATGGAGGCGGGAAATATCAAGACCGG‑CGTTTAACGCTGAAAATTGGACT‑ATTTTGCTTCCCAGCGGGCCGGGGGCTGGCGAAACGTCGGCCTGCGTTGAGCGATAGCGATGT‑CCATTTTCCCTTATAACCGCCAGCAAAGGAACTTCCCAGAGGATGAc                                                                                                                                                             <  1:46084/249‑1 (MQ=255)
                                                cAGGATTCCTGTGGCAGACGCTTATCGGTGGTTTCTGGTACTGCGGGAGAAAAGC‑‑GGCGCGGGGTGCCAAAATCCACAAGAGGGAAGCGCAACATAT‑AAGGCCGGGCGTTTAACGCAGAGAAGCCGAAT‑AATTTGTTTCCCGGGGTGGCGAGGGCTTGCGCAACGTCGGCCGGGGGGGAGCAATAGCGATGTCCCTTTC‑CCCTGTTAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACgg                                                                                                                                                 <  1:1683/251‑1 (MQ=255)
                                                                      tAACGATGGTTTCTTGTACTTCGCGAGAAAAAC‑‑GGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGG‑CGGTTAACGACGAGAAGTCGACT‑AAATTTCTTTCCAGCGTGTCGGGGGGTTTCGCAACCTCGGCAGGCGGGGAGCGGTAGCGGAGTACCATTT‑TACTGATAACCGCCCACAACAAAACTTCCCCACGGACGACTTCAAGCCCCGGGCCGCTTAAACCAATATCCAGc                                                                                                                           >  1:60900/1‑251 (MQ=255)
                                                                                     gtactccggaaaaaaaaC‑‑GGGGCCGGGTGCAAAAATCCACAAGAGGGAAGCGCGACATAACAAGACGGG‑GGGGTACCGCCGAGAAGTCGCCA‑ATTTTGCTTTCCAGTGTGTCGGGGGTTGGCGAACCGCCGGCCTGCGGGGAGCAAAAGCGATGTCCCATTT‑CACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTg                                                                                                            <  1:53577/239‑1 (MQ=255)
                                                                                            tcgAGAAAAACAGGGCGAGGGGTGAAAAACTCCATGAGATGGAGGCGCCACATATAAAGCCCGG‑CGTTTAACGCAGAGAGGCGGCCA‑AATTTCCTTTCTGGGGGTTCGAGGTCTTGCGCACGGTCGCCTTGGGGGGAGCAAAAGAGTTGTCCCATTT‑CACGGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGac                                                                                                         <  1:153098/248‑1 (MQ=255)
                                                                                                                                     gAAGCGCGACATATCAAGACCGG‑CGGTTAAAGCCGAGAAGTCGACT‑AATTTTCTTTCCCGCGTGTGGAGGGCTTGCGCCACGTCGGCCTGCGGGGGGCGAGAGCGGTGTACCAATT‑CACTGATAACCCCCAGCAAAGGAACTTTCCCGAGGATGACTTCACCCCACGGGCCGCTTTAAAGAATATCCAGCTTGCCCTTATCGTTGGACACGGTGACTTTTTCCGGGGTAAAGCGGAACCCGCGTAACCAGt                                                              >  1:185873/1‑251 (MQ=255)
                                                                                                                                                                       cGAGAAGTCGACT‑AATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTT‑CACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcagggc                               >  1:118681/1‑249 (MQ=255)
                                                                                                                                                                            aGTCGACT‑AATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTT‑CACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGc                        >  1:98176/1‑251 (MQ=255)
                                                                                                                                                                                gACT‑AATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTT‑CACTGATAACCGCCAGCAAAGGAACTTCCCAGCGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAACGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaa                                           >  1:208327/1‑228 (MQ=255)
                                                                                                                                                                                                   ggtggggggggtttGCGCAACGTCGGCCGGGGGGGAGCGATAGCGATGTCCCATTT‑CACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCa  <  1:119672/239‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                          
GTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAAC‑‑GGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGG‑CGGTTAACGCCGAGAAGTCGACT‑AATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTT‑CACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCA  >  NC_000913/989719‑990160

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