Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,228,2220TG62.5% 0.1 / 11.6 16pseudogene (388/1287 nt)mdtQpseudogene; putative channel/filament proteins
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base G (10/0);  total (13/3)
Fisher's exact test for biased strand distribution p-value = 3.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-02

AAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGA‑GAT‑GCC‑AGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGG  >  NC_000913/2228075‑2228349
                                                                                                                                                   |                                                                                                                                   
aagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCCGCCAGCAATTGGCGGGGTTg                                                                                                                                  >  2:15065/1‑151 (MQ=255)
             cTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCT‑CGCGTTccgc                                                                                                                     >  1:57178/1‑150 (MQ=255)
               gATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCCACGCCGCCAGCCCGCAATTGGCGGGGTTgtccccggttccccc                                                                                                                    >  1:214493/1‑136 (MQ=255)
                                gTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCC‑CGCGTTTCGCCGCCCCTGCTTTaaaccta                                                                                                  >  1:377748/1‑145 (MQ=255)
                                 tCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCCGTGCTTTAACCGTAcc                                                                                                <  2:296135/151‑1 (MQ=255)
                                      aTAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAgg                                                                                            <  2:29629/150‑1 (MQ=255)
                                        aaTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCT‑CCCGTTCCCCCCCCCGTGCTTTTACCCCACCCCGgggg                                                                                         >  2:99642/1‑151 (MQ=255)
                                                 cTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAAcc                                                                                <  2:89010/151‑1 (MQ=255)
                                                   gCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGGGTTGCT‑CCCGTTCCGCCGCCCCTGCCTTTACCCGACCCCGGCGGGCAGTAACCt                                                                               >  1:372726/1‑150 (MQ=255)
                                                                               gtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGGTGCCAGCCCGCAATTGGCGGGGTTTAC‑CGGGGTTCGCAGCCCGGGGTTTAACCGGACCCAGGCGGGGAGCAACCCCCCCCCCCCTTTTTCCCCAGATATaga‑gg                                                 >  2:205336/1‑150 (MQ=255)
                                                                               gtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCT‑CGCGGTCCGCCGCCCCGGCTTTTAACGTACCCCAGCGGGCAGGAACCCCCCCCCCCTTCTTTTCCCCGagattcagg                                                  >  1:149572/1‑146 (MQ=255)
                                                                                                   gCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCCGTTCC‑CCGCCCCCGCTTTAAACCGACCCAGGCGGGGCGTAACCTCCGCCCGATTCTTTTCCCCGATATCCA‑GAT‑GCC‑AGcccccccgttagc                            >  1:66813/1‑143 (MQ=255)
                                                                                                       cTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCT‑CGCGGTCCGCCGCCCCGGCTTTTACCCCACCCCGGCGGGCAGTAAACTCCGCCCCGTTCTTTTCCCCGATATCCAGGAT‑CCC‑A‑CCCGCCCTTAAACCcaaa                       >  1:287798/1‑151 (MQ=255)
                                                                                                             gTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTTT‑CCGGTTCCGCCGCCCGTGCTTTAACCGTACCCCGGGGGGGGGGAAACCCCGCCCCGTTCTTTTCCCCGATATCCG‑GAT‑GCCAACCCCCCCGTTAAGCCAAAAGAAcg                  >  2:47413/1‑150 (MQ=255)
                                                                                                             gTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCT‑CCCGGTCCCCCGCCCCGGCCTTTAACGCAACCCGGGGGGGAGGAACCCCCCCCCCATTCTTTCCCCCGATATCGA‑GATGGCA‑ACCCCCCCGTTAACCCaaaaaaacg                  >  1:308969/1‑146 (MQ=255)
                                                                                                                             cTGCCAGCCAGCAATTGGCGGGTTTGCG‑CGCGTTCCGCCGCCCGTGCTTTAACCCTACCCCGGGGGGCAGGAAACTCCGCCCCACTCTTTCCCCAACAATCGA‑CAT‑GGC‑ATCCCGCCGTTCACCCCAAAGTACCGTTCGTGTACCCCggg  >  1:255625/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                   
AAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCT‑CGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGA‑GAT‑GCC‑AGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGG  >  NC_000913/2228075‑2228349

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