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

AAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGG  >  NC_000913/2228075‑2228349
                                                                                                                                                   |                                                                                                                                   
aagTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCCGCCAGCAATTGGCGGGGTTg                                                                                                                                  >  2:15065/1‑151 (MQ=255)
             cTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGTTccgc                                                                                                                     >  1:57178/1‑150 (MQ=255)
               gATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCCACGCCGCCAGCCCGCAATTGGCGGGGTTgtccccggttccccc                                                                                                                    >  1:214493/1‑136 (MQ=255)
                                gTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCCCGCGTTTCGCCGCCCCTGCTTTaaaccta                                                                                                  >  1:377748/1‑145 (MQ=255)
                                 tCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTAcc                                                                                                <  2:296135/151‑1 (MQ=255)
                                      aTAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAgg                                                                                            <  2:29629/150‑1 (MQ=255)
                                        aaTGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCGTTCCCCCCCCCGTGCTTTTACCCCACCCCGgggg                                                                                         >  2:99642/1‑151 (MQ=255)
                                                 cTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAAcc                                                                                <  2:89010/151‑1 (MQ=255)
                                                   gCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGGGTTGCTCCCGTTCCGCCGCCCCTGCCTTTACCCGACCCCGGCGGGCAGTAACCt                                                                               >  1:372726/1‑150 (MQ=255)
                                                                               gtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGGTGCCAGCCCGCAATTGGCGGGGTTTACCGGGGTTCGCAGCCCGGGGTTTAACCGGACCCAGGCGGGGAGCAACCCCCCCCCCCCTTTTTCCCCAGATATagagg                                                 >  2:205336/1‑150 (MQ=255)
                                                                               gtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCGCGGTCCGCCGCCCCGGCTTTTAACGTACCCCAGCGGGCAGGAACCCCCCCCCCCTTCTTTTCCCCGagattcagg                                                  >  1:149572/1‑146 (MQ=255)
                                                                                                   gCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTGCTCCCCGTTCCCCGCCCCCGCTTTAAACCGACCCAGGCGGGGCGTAACCTCCGCCCGATTCTTTTCCCCGATATCCAGATGCCAGcccccccgttagc                            >  1:66813/1‑143 (MQ=255)
                                                                                                       cTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTTCTCGCGGTCCGCCGCCCCGGCTTTTACCCCACCCCGGCGGGCAGTAAACTCCGCCCCGTTCTTTTCCCCGATATCCAGGATCCCACCCGCCCTTAAACCcaaa                       >  1:287798/1‑151 (MQ=255)
                                                                                                             gTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTTTCCGGTTCCGCCGCCCGTGCTTTAACCGTACCCCGGGGGGGGGGAAACCCCGCCCCGTTCTTTTCCCCGATATCCGGATGCCAACCCCCCCGTTAAGCCAAAAGAAcg                  >  2:47413/1‑150 (MQ=255)
                                                                                                             gTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGGTTTCTCCCGGTCCCCCGCCCCGGCCTTTAACGCAACCCGGGGGGGAGGAACCCCCCCCCCATTCTTTCCCCCGATATCGAGATGGCAACCCCCCCGTTAACCCaaaaaaacg                  >  1:308969/1‑146 (MQ=255)
                                                                                                                             cTGCCAGCCAGCAATTGGCGGGTTTGCGCGCGTTCCGCCGCCCGTGCTTTAACCCTACCCCGGGGGGCAGGAAACTCCGCCCCACTCTTTCCCCAACAATCGACATGGCATCCCGCCGTTCACCCCAAAGTACCGTTCGTGTACCCCggg  >  1:255625/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                   
AAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGGTTTGCTCGCGTTCCGCCGCCCGTGCTTTAACCGTACCCAGGCGGGCAGTAACCTCCGCCCGATTCTTTCCCCAGATATCGAGATGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTACCACGGG  >  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: