Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,951,0180CA40.0% 11.0 / 12.1 20intergenic (‑123/+34)fliL/fliKflagellar protein FliL/flagellar hook‑length control protein FliK
Reads supporting (aligned to +/- strand):  ref base C (5/7);  new base A (5/3);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCC  >  NC_002947/4950888‑4951154
                                                                                                                                  |                                                                                                                                        
ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaa                                                                                                                                 <  2:86945/140‑1 (MQ=255)
ttCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaa                                                                                                                                 >  1:86945/1‑140 (MQ=255)
             taatCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGtt                                                                                                          <  1:777078/150‑1 (MQ=255)
                         cATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                <  1:406009/148‑1 (MQ=255)
                         cATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgta                                                                                                >  2:406009/1‑148 (MQ=255)
                               gggTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGTATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGCGATCTCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGa                                                                                        <  1:396609/150‑1 (MQ=255)
                                            aaTGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGcc                                                                           >  2:133086/1‑150 (MQ=255)
                                                 cttggcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACaa                                                                      <  1:133086/150‑1 (MQ=255)
                                                  ttggcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAac                                                                     <  2:185135/150‑1 (MQ=255)
                                                  ttggcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaagaa                                                                                                                              >  1:416707/1‑93 (MQ=255)
                                                  ttggcttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCGAGCCCGCgaagaa                                                                                                                              <  2:416707/93‑1 (MQ=255)
                                                       ttgAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGACGACAACCACct                                                                >  1:840726/1‑150 (MQ=255)
                                                                                       atggatTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgt                                                                                                 >  1:192672/1‑85 (MQ=255)
                                                                                       atggatTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCgt                                                                                                 <  2:192672/85‑1 (MQ=255)
                                                                                            ttCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGc                           <  2:840726/150‑1 (MQ=255)
                                                                                                                  ggCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGgcgc     >  1:905731/1‑150 (MQ=255)
                                                                                                                     ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc  >  1:709043/1‑150 (MQ=255)
                                                                                                                     ttcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCAGGGAAACATTCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGcgcc  >  2:711345/1‑150 (MQ=255)
                                                                                                                      tcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGc            <  2:458716/139‑1 (MQ=255)
                                                                                                                      tcttcGCGGGCGAGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGc            >  1:458716/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                        
TTCGCCATGCCAATAATCCGTCGTCCATCGGGGGTTCAAAGGAGAATGGCTTGGCTTGAGCAAGTGTTATGCCAGATTTTGCAGGGGATGGATTCCAGCAGGTTTTGCGGTGCTGGCTTCTTCGCGGGCTCGCCCGCGAAGAAGCCAGCCCAGGCCAACAGTTATCAGGCGTAGTAATCGACCGCACTGTCGCCGACAACCACCTGCTGCTCCAAAGGCCGGGCGCTATCGACCGCTTCGGCCCCACCCTGCTCGGCCCGGCGCGCC  >  NC_002947/4950888‑4951154

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: