Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,135,9750GC41.2% 7.9 / 10.9 18A138P (GCC→CCC) flgHflagellar protein of basal‑body outer‑membrane L ring
Reads supporting (aligned to +/- strand):  ref base G (5/5);  new base C (0/7);  total (5/13)
Fisher's exact test for biased strand distribution p-value = 4.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACG  >  NC_000913/1135922‑1136038
                                                     |                                                               
cTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAAt                                                   >  2:120941/1‑68 (MQ=255)
cTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAAt                                                   >  2:308445/1‑68 (MQ=255)
 ccGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATa                                                  <  2:64145/67‑1 (MQ=255)
  cGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATAc                                                 <  1:100063/68‑1 (MQ=255)
                aGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGt                                   >  1:354317/1‑68 (MQ=255)
                      cGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTtgacgg                             <  2:687561/68‑1 (MQ=255)
                      cGGTGGTAACACGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTATCGGCACGTtgacgg                             <  2:98424/68‑1 (MQ=255)
                      cGGTGGTAACACGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTtgacgg                             <  2:367625/68‑1 (MQ=255)
                       ggtggtAACACGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTtgacggt                            <  1:230583/68‑1 (MQ=255)
                                 cggtCAACGGAAAGGGCGGGTCCAATGCCAGCAATACCTTTAGCGGCACGTTGACTGTGACGGTTGAc                  <  2:50310/65‑1 (MQ=255)
                                 cGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGAc                  <  1:430572/68‑1 (MQ=255)
                                    tCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAg               <  2:289358/68‑1 (MQ=255)
                                     cAACGGAAAGGGCGGGGCCATGGCCAGCAATACCTTTAGCGGCACGTTGACGTTGACGGTTGACCAgg              <  2:616246/68‑1 (MQ=255)
                                           aaaGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTgg        <  2:267568/68‑1 (MQ=255)
                                               ggCGGGGCAAATGCCAGCAATACCTTTAGCGGCAAGTTGACGGTGACGGTTGACAAGGTACTGGTaaa    >  2:174233/1‑68 (MQ=255)
                                                gCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGCCCAGGTACTGGTCAAc   >  1:181129/1‑68 (MQ=255)
                                                gCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAAc   <  1:321871/68‑1 (MQ=255)
                                                 cGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACg  <  1:621422/68‑1 (MQ=255)
                                                     |                                                               
CTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACG  >  NC_000913/1135922‑1136038

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