Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA1,250,9611.C26.7% 27.9 / 15.4 15intergenic (‑123/‑105)dhaK/dhaRdihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit/dhaKLM operon transcription activator
Reads supporting (aligned to +/- strand):  ref base . (6/5);  new base C (2/2);  total (8/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 = 8.24e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGATGGTCGGGGCATATCTCCCTT  >  NC_000913_3_pae_tpiA/1250830‑1251109
                                                                                                                                    |                                                                                                                                                    
ttttttCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACataatttgg                                                                                                                                    <  2:19917/150‑7 (MQ=255)
                               tAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATAGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACgtg                                                                                                     >  1:87282/1‑150 (MQ=255)
                                   gCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCGGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTAtt                                                                                                 >  3:60972/1‑150 (MQ=255)
                                                   aCGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCACTATGAAACataatt                                                                                                                                       >  1:72650/1‑93 (MQ=255)
                                                   aCGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCACTATGAAACataatt                                                                                                                                       <  2:72650/96‑4 (MQ=255)
                                                    cGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGa                                                                                <  2:58935/150‑1 (MQ=255)
                                                          ggCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGa                                                                          >  1:30952/1‑150 (MQ=255)
                                                          ggCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGa                                                                          >  1:46676/1‑150 (MQ=255)
                                                                            aTGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAgg                                                        >  1:102642/1‑150 (MQ=255)
                                                                                        tgCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGatat                                            <  2:102642/150‑1 (MQ=255)
                                                                                                gCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGgcg                                    <  2:30952/150‑1 (MQ=255)
                                                                                                                    taaagcatATGTTTCACTATGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGATGGTCgg                >  4:64304/8‑150 (MQ=255)
                                                                                                                      aagcatATGTTTCACTATGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGATGGTCgggg              <  3:64304/145‑1 (MQ=255)
                                                                                                                        aaTTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGATGGTCGGGGCa            >  4:20024/1‑150 (MQ=255)
                                                                                                                                  cATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGATGGTCGGGGCATATCTCCCtt  <  4:60972/150‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                                    
TTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGATGGTCGGGGCATATCTCCCTT  >  NC_000913_3_pae_tpiA/1250830‑1251109

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: