Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,7470TG53.3% 3.1 / 11.2 16T255P (ACC→CCC) citGtriphosphoribosyl‑dephospho‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base G (8/0);  total (11/5)
Fisher's exact test for biased strand distribution p-value = 6.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAA‑CGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAAT‑GCCAGTTCAGG‑ATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGT  >  NC_000913/646499‑646994
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           
cGGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCCATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGACAATTAATATTTCCAAATTTATCACGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGtgc                                                                                                                                                                                                                                                           >  2:11643/1‑249 (MQ=255)
                                                                 tCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGAGGATCGCAAGGCCCCCTTTTT‑CAATAATAGTTGGGCACCCCGCTGTAGCCCGCGCCGGCCCCCCTCGccc                                                                                                                                                                                         >  1:334382/1‑250 (MQ=255)
                                                                                aTACGGGTTCTCCGAAAATTAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTTG‑GATGCCCCCCTTTTGTACAACTGTTTTCGTCTCGCGCTGTAGCAACCGCAGGgcccccccccccccccgcgcagcaacg                                                                                                                                                                          >  1:6320/1‑230 (MQ=255)
                                                                                                    aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑GATGCCCCCCTTTTTCAATAATGTTTTCGCCTCCGGCTGTAGCCACCGCAGGGCCCCCCCGCCACCCCCCGCGGGAA‑CGTTGGTATCGCCGTTGAttgc                                                                                                                                                     >  1:262202/1‑248 (MQ=255)
                                                                                                    aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑AAAGCCCCCCTTTTGCAATAATGTTTTCGCCTCCCGCTGTAGCCACCGCCGGCCCCCCCCGCCCCCCCGCGCTGCCA‑CGTTgggtttgccgctgttccc                                                                                                                                                     >  1:218678/1‑235 (MQ=255)
                                                                                                             aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAA‑CGTTGGTATCGCCGTTGATCGCCATCAGTAg                                                                                                                                            <  1:245280/251‑1 (MQ=255)
                                                                                                                  tATCAAGTTCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCCCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAA‑CGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGc                                                                                                                                        <  2:175157/250‑1 (MQ=255)
                                                                                                                       aGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGG‑GAAGCCCCCCTTTTTCAATAAAGTTTTCGCCTCCCGCTGTAGGCCGCGCAGGCCCCCCCCGGCACCGCCCGAGGCAAAGGTTGG‑ATAGCCGTTGATCggcctccgatagggccaggg                                                                                                                                   >  1:234595/1‑231 (MQ=255)
                                                                                                                               aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑AATGCCCCCTTTTTTCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCCCGCCCCCCCGCGATGCCA‑CGGTGGTATCGCCGCTGGTCGCCATCCGTAGGAGGAAGGGATCGAGCa                                                                                                                           >  1:93377/1‑250 (MQ=255)
                                                                                                                                                       cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑AATGCCCCCCTTTTTCAATAATGTTTTCGCCTCGCGCTTTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGGTGCAA‑CGTTGGGATCGCCGTTGTTCGCCCTCAGTAGGAGCAAGGGATCGAGCCAT‑GCCCGCTCAGG‑GTCTAAcccccc                                                                                                >  1:69882/1‑249 (MQ=255)
                                                                                                                                                               gATCAGTAGGTCAGCACTGCCGACTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑AATGCCCCCCTTTTGCAATAATGTTTGCGCCTCGCGCGGTAGCCAGCGCAGGCCCCCCTCCCCACCCCGCGCTGCAA‑CGTTGGGATAGCCGGTGGTTGCCCCCAGTAGGGGGAAGGGATTGAGCAAA‑GCCCGTTAAGG‑GTCTAAcccccgcaccagcc                                                                                        >  1:85743/1‑241 (MQ=255)
                                                                                                                                                                      aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAA‑CGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAAT‑GCCAGTTCAGG‑ATCTAACCCCTGATCCAGCAGAGTGAg                                                                                 <  1:118753/251‑1 (MQ=255)
                                                                                                                                                                              aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑GATGCCCCCTTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGCCAGCGCAGGCCCCCCTCCCCACCCCCCCATGCAAACGT‑GGGATCTCCGGTGGTCGCCCTCAGTAGGAGCAAAGTATTGAGGAATTGCC‑GTTCAGG‑ATTTTAACCCTGCTCCCGCAGGGGGGGGTAATGCg                                                                         >  1:109767/1‑251 (MQ=255)
                                                                                                                                                                                                   tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAA‑CGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAAT‑GCCAGTTCAGG‑ATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCAcc                                                     <  2:218678/250‑1 (MQ=255)
                                                                                                                                                                                                    cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAA‑CGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAAT‑GCCAGTTCAGG‑ATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCa                                                    <  1:180838/250‑1 (MQ=255)
                                                                                                                                                                                                         gATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCG‑ACTGCCCCCCTTTTTCCATAATGTTTGCGCCCCGCGCTGTAGCCCGCGCAGGGCCCCCTCCCCACCCCGCGATGCAA‑CGTTTGTATCCCCGGTGGTCGGCCTCCGTCAGAGCAAGGGATCGAGCAAA‑GCCAGTGCAGGGACATAACCC‑TGATCCCACAGAGTGAGGTAATGGGGCAAGGCGGGGTTGATCACCatcgta                                               >  2:12833/1‑245 (MQ=255)
                                                                                                                                                                                                                                                     gggTTCG‑AATGCCCCCTTTTTGCCATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAA‑CGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAAT‑GCCAGTTCAGG‑ATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGTTTGATCACCAGTGGATAACCCGCTTCTGCTTCCCCGCGTGCGCCGGTAAGGCCAAGCTGt  >  2:298775/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           
CGGGGACAGGGATAAGAAACATGACACCCATCACCACCAGTGGGGCCAATAGTTTCCATATATTATCTTTTGCTAAAGACATACGGGTTCTCCGAAAATTAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCG‑AATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAA‑CGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAAT‑GCCAGTTCAGG‑ATCTAACCCCTGATCCAGCAGAGTGAGGTAATGCGGCAAGGCGTGATTGATCACCAGTGGATAACCCGCTTCGGCTTCACCGCGTGCGCCGGTAAGGCCAAGCTGT  >  NC_000913/646499‑646994

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