Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,020,842 C→T R252R (CGC→CGT fliM → flagellar motor switch protein FliM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,020,8420CT83.3% 16.2 / ‑0.8 13R252R (CGC→CGTfliMflagellar motor switch protein FliM
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (6/4);  total (9/4)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

ATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGC  >  NC_000913/2020657‑2021083
                                                                                                                                                                                         |                                                                                                                                                                                                                                                 
atCACCACCTCGCCGCACGACATTGTGGTTAACACGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCGACTTTAATATACGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGTCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGc                                                                                                                                                                                   >  1:173736/1‑250 (MQ=255)
  cccccccccgccaaacgaagttGTGGTTACCCCGCCTTTCCTGGTGAATATTGGCCACCTGACCGGCGAATTTAATATCTCCCGCCCTTTCAGCATGCTCGAGCCGCTACGGGAATTGTGGTTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGTCAGGTGCAGCATGCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAg                                                                                                                                                                                 <  1:293977/230‑1 (MQ=255)
                                      gTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCCTTCAGCATGATCGAGCCGCTAAGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAAACCCGGCGACGTCCTGCCGAt                                                                                                                                              >  1:207345/1‑249 (MQ=255)
                                                                  tgCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCAGCTACGGAAATTGTAGGTTAACACGCAGCTGGAAAACTCGCGTAATGAAGATCAGAAATGGCGCGATAACCTGGTGCGACAGGTGCAGCCTTCCCAGCTGGCGATGGTCGCCAACTTTGCCGATATCTCGCTACGACTGTCg                                                                                                                                                                                    >  1:231773/2‑183 (MQ=255)
                                                                                            naGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGTCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGCTATCTCGCTACGCCTGTCGCAGATTTTAAAAc                                                                                                                                                                       >  1:261087/2‑170 (MQ=255)
                                                                                            cAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGTCAGGTGCAGCATTCACAGCTGGAGCTGGTCGcc                                                                                                                                                                                                                  >  1:9586/1‑127 (MQ=255)
                                                                                                                          ggTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATAGGAACTGGCGCGATAACCTGGTGCGTCAGGTGCAGCCTTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCg                                                                                                                                                                                    >  1:237499/1‑127 (MQ=255)
                                                                                                                           gataccnnGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGTCAGGTGCAGCATTCACAGCTGGAGCTGGTCGcc                                                                                                                                                                                                                  <  2:9586/91‑1 (MQ=255)
                                                                                                                                              aCTCGCGTAATGAAGATCATAACTGGCGCGATAACCTGGTGCGTCAGGTGCAGCATTCACAGCTGGGGCTGGTCGCCAACTTTGCGGATATCGCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGGCGTCCTGCCGATAGAAAATCGCGATCGCATCATCGCTCAGGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCGCGCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGAtt                                    <  1:280791/251‑1 (MQ=255)
                                                                                                                                               cTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAACACTGAACCCCGGCGACGTCCTGCCg                                                                                                                                                >  1:298913/1‑142 (MQ=255)
                                                                                                                                                             aTCAGAACTGGCGCGATAACCTGGTGCGTCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTcccn                                                                                                   <  1:61926/173‑2 (MQ=255)
                                                                                                                                                                         gcgATAACCTGGTGCGTCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATCTCTCGCTACGCCTGTAGCAGATTTTaa                                                                                                                                                                          >  1:41645/1‑90 (MQ=255)
                                                                                                                                                                                ccTGGTGCGTCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGc  >  1:310886/1‑251 (MQ=255)
                                                                                                                                                                                         |                                                                                                                                                                                                                                                 
ATCACCACCTCGCCGAACGACATTGTGGTTAACACGCCGTTCCATGTGGAGATTGGCAACCTGACCGGCGAATTTAATATCTGCCTGCCATTCAGCATGATCGAGCCGCTACGGGAATTGTTGGTTAACCCGCCGCTGGAAAACTCGCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGC  >  NC_000913/2020657‑2021083

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