Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,135,9750GC39.3% 27.9 / 25.3 30A138P (GCC→CCC) flgHflagellar protein of basal‑body outer‑membrane L ring
Reads supporting (aligned to +/- strand):  ref base G (12/5);  new base C (0/11);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 2.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACG  >  NC_000913/1135912‑1136038
                                                               |                                                               
tttGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGTCCaa                                                             <  2:113957/68‑1 (MQ=255)
      aaCGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAGCGGAAAGGGCGGGGCCAATGCCAg                                                       >  2:420090/1‑68 (MQ=255)
      aaCGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGccgg                                                       >  1:847404/1‑66 (MQ=255)
      aaCGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAg                                                       >  2:3614/1‑68 (MQ=255)
      aaCGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAg                                                       >  2:328009/1‑68 (MQ=255)
      aaCGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAg                                                       >  2:180700/1‑68 (MQ=255)
       aCGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGCCCAATGCCAGc                                                      <  1:508010/68‑1 (MQ=255)
         gCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGCCCAATGCCAGCaa                                                    <  2:823782/68‑1 (MQ=255)
          cTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAAt                                                   >  2:599509/1‑68 (MQ=255)
                 cGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTa                                            <  2:614674/68‑1 (MQ=255)
                         aaGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACg                                    <  1:1145875/68‑1 (MQ=255)
                         aaGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACg                                    <  2:648884/68‑1 (MQ=255)
                                      tAACACGTTCAACGGAAAGGGCGGGGCCAATGCCTGCAATACCTTTAGCGGCACGTtgactgtgacgg                       <  1:367572/68‑1 (MQ=255)
                                        acacGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGtt                     <  2:1134175/68‑1 (MQ=255)
                                         cacTTTCAACGGAAAGGGCGGGTCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTg                    <  2:852034/68‑1 (MQ=255)
                                         cacGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAAAACCTTTAGCGGCACGTTGACGGTGACGGTTg                    >  2:564410/1‑68 (MQ=255)
                                           cGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGAc                  <  2:252107/68‑1 (MQ=255)
                                           cGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGAc                  <  1:642517/68‑1 (MQ=255)
                                           cGTTCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGAc                  <  1:628109/68‑1 (MQ=255)
                                           cGTGCAACGGAAAGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGAc                  <  2:524323/68‑1 (MQ=255)
                                               cAACGGAATGGGCGGGCCCAATGCCAGCAACACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAgg              <  2:1113956/68‑1 (MQ=255)
                                                aaCGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGCCGGTTGACCAGGt             >  2:573006/1‑68 (MQ=255)
                                                aaCGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGt             <  1:3614/68‑1 (MQ=255)
                                                   ggAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGCCGGTTGACCAGGTACt          >  2:767657/1‑68 (MQ=255)
                                                   ggAAAGGGCGGGGCCAATGCCAGCAAAACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACt          >  2:840931/1‑68 (MQ=255)
                                                     aaaGGGCGGGCCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTgg        <  2:83144/68‑1 (MQ=255)
                                                      acgggCGGGGCCAATGCCAGCAATACCTTAAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGt       <  1:1028178/66‑1 (MQ=255)
                                                        gggCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTg                    >  1:938872/1‑53 (MQ=255)
                                                        gggCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTg                    <  2:938869/53‑1 (MQ=255)
                                                           cGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGCCGGTTGACCAGGTACTGGTCAACg  >  1:2496/1‑68 (MQ=255)
                                                               |                                                               
TTTGGTAACGCTCGTGCCGATGTCGAAGCCTCCGGTGGTAACACGTTCAACGGAAAGGGCGGGGCCAATGCCAGCAATACCTTTAGCGGCACGTTGACGGTGACGGTTGACCAGGTACTGGTCAACG  >  NC_000913/1135912‑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: