Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,2390GC53.8% 1.3 / 13.5 13pseudogene (371/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base C (4/3);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCG‑TTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCG‑TTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGA‑GCAAACGGCCCCATTAAGCCTTCTGCCGACA‑TTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCA  >  NC_000913/2228020‑2228466
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                       
aGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCG‑TTCCGCCGCCCGTGCtt                                                                                                                                                                                                                                <  1:788738/229‑1 (MQ=255)
          gggTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCCGCCAGCAATTGGCGGGTTTTTTCCCG‑TTTCCCCCCCCCCGCTTTTAccccacc                                                                                                                                                                                                                      >  2:788738/1‑224 (MQ=255)
                                  tttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCG‑TTCCGCCGCCCg                                                                                                                                                                                                                                     <  1:665693/190‑1 (MQ=255)
                                   ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCG‑TTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTccc                                                                                                                                                                                    <  1:604828/238‑1 (MQ=255)
                                   ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGAGGGCTCCGCTGACAGCCAGCAATTGGCGGGTTTGCTCGCG‑TTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTccc                                                                                                                                                                                    <  1:801121/238‑1 (MQ=255)
                                                                   tCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCG‑TTCCGCCCCCCGGGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCCATTCTTTCCCCCGATATCGAGAAGCCCGCCCCCCG‑CTTAGCCa                                                                                                                                                  >  1:282063/1‑239 (MQ=255)
                                                                                                           ctctTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCG‑TTTCGCCGCCCGTGCTTTTACCGTACCCCGGCGGGCAGTAACCTCCCCCCGCTTCTTTCCCCAGATATCGAGATTCCAGCCCGCCG‑CTAAGCCAAAAGTACCCTTGGTGGA‑CCACCGGCCGGGCGTACCTTCg                                                                                                          >  1:756888/1‑238 (MQ=255)
                                                                                                                                                                cccAGTGCGGGCGGGTTACGCTGCCAGCCAGCAATGGGCGGTTTTGCTCGGG‑TCCCGCCCCCTGTGCTTAACCCGCCCCCGGGGGGGCGGTACCCCCCGCCCTTTTCTTCCCCAAGATAGCGAGCTGCCAGCCCGCCG‑TTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGA‑GCAAACGGCCCCATTAAGCCTTCTGCCGACA‑TTTTTTGc                                                  <  2:162956/239‑1 (MQ=255)
                                                                                                                                                                   aGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGTTTCCCCCCCCCGGCTTTTAACCG‑ACCCCCGGGGGGAGCAAACCCCGCCCCATTCTTTTCCCCGATATTGAGATGGCAGCCCCCCGTTTAACCCAAAAAAACC‑CTGGGGTACCCAGGGGCCGG‑CGGACCCTCGGGCGGGGCGGTTACAGGAAAAAGGCCCCCTTAAGCCTTCTCCGGCCA‑tttttttcg                                                 >  2:370611/1‑235 (MQ=255)
                                                                                                                                                                                       ccgcccagcAGTTGGGGGTTTTCCCCGCG‑CCCCCCCCCCCGTGCTGTAACGCGCCCGAGGCGGCCGGTACCCCCCCCCCGTTTCTTCCCCCGGATACCGGGATCCCCCCCCGCCG‑TTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGA‑GCAAACGGCCCCATTAAGCCTTCTGCCGACA‑TTTTTTGCCGCTCCATATCCGCCGAAAAGTc                           <  2:756888/235‑1 (MQ=255)
                                                                                                                                                                                            cagcAATTGGCGGGTTTTCTCGCG‑GTCCGCCCCCCCTGCTTTAACCGGACCCAGGCGGGCAGTAACCTCCCCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCG‑TTTAACCAAAAGTACCCCCGGTGTA‑TCCCGGGCCCGGCGGACCCGCCGGCGGGTCGTTCCGA‑GCCAAAGGGCCCCTTCAGGCTTCTTTCGCCACTTTTTTGCCGC‑CCCTATCCGCCGCAAAGACgc                         >  2:271002/1‑235 (MQ=255)
                                                                                                                                                                                                 ttgggggggTTTCCCCGCT‑TCCCCCCCCCGGTGTTTTACCCGTCCCCGGGCGGGCATTACCCCCCCCCCGATTCTTCCCCCAGATACCGAGATGCCCGCCCGCCG‑TTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGA‑GCAAACGGCCCCATTAAGCCTTCTGCCGACA‑TTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGa                 <  1:370611/233‑1 (MQ=255)
                                                                                                                                                                                                 aTTGGCGGGTTTGCTCGCG‑TTCCGCCGCCCGTGCTTTAACCGTACCCCGGCGGGCAGTAAACTCCGCCCGATTCTTTCCCCAGATATCGCGATGCCAGCCCGCCG‑TTAAGCCAAAAGTACCGTTGGGGTA‑CCACGGGCCGGTCGTACCTGCGGCCGCGCCGTTCCGA‑GCAAACCGCCCCCTTAAGCCTTCCTCCGGCA‑TTTTTTTCCGCCCCACTTTCGCCCCAAAGTAGACATGCgg                  >  2:822897/1‑238 (MQ=255)
                                                                                                                                                                                                                       ccgccgCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCG‑TTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGA‑GCAAACGGCCCCATTAAGCCTTCTGCCGACA‑TTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCa  >  1:934120/1‑233 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                       
AGCTGCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCG‑TTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCG‑TTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGA‑GCAAACGGCCCCATTAAGCCTTCTGCCGACA‑TTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCA  >  NC_000913/2228020‑2228466

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: ‑