Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,017,557 A→G Q335R (CAG→CGG)  fliI → flagellum‑specific ATP synthase FliI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,017,5570AG94.1% 52.0 / ‑3.8 17Q335R (CAG→CGG) fliIflagellum‑specific ATP synthase FliI
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (6/10);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

ACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCG  >  NC_000913/2017447‑2017676
                                                                                                              |                                                                                                                       
aCCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCggg                                                                                            <  1:285838/140‑1 (MQ=255)
  cGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCGGGCg                                                                                          <  1:112402/140‑1 (MQ=255)
   gTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCggg                                                                                            <  2:401625/137‑1 (MQ=255)
   gTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCggg                                                                                            >  1:401625/1‑137 (MQ=255)
   gTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCGGGCGa                                                                                         >  1:17810/1‑140 (MQ=255)
          tttGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTCCCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGa                                                                                  <  2:164000/140‑1 (MQ=255)
          tttGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGa                                                                                  <  2:17810/140‑1 (MQ=255)
              ccAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGt                                                                              <  2:120635/140‑1 (MQ=255)
                         gCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTc                                                                                    >  2:216025/1‑123 (MQ=255)
                         gCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTc                                                                                    <  1:216025/123‑1 (MQ=255)
                               gTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCg                                                             >  1:219422/1‑140 (MQ=255)
                               gTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGc                                                            >  1:97861/1‑140 (MQ=255)
                               gTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCg                                                             >  2:169759/1‑140 (MQ=255)
                                                  gCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCggg                                          <  1:260194/140‑1 (MQ=255)
                                                   cATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGc                                         >  1:24299/1‑140 (MQ=255)
                                                            cggcggcTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCgg                                <  2:219422/140‑1 (MQ=255)
                                                                                      ccGTGCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGc      <  1:218693/140‑1 (MQ=255)
                                                                                          gCTCACTGAAGGCGATGACCGGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCcgcg  <  2:201770/140‑1 (MQ=255)
                                                                                                              |                                                                                                                       
ACCGTCGGTGTTTGCCAAATTACCGGCACTGGTCGAGCGTGCCGGAAATGGCATTAGCGGCGGCGGCTCGATTACCGCGTTTTATACCGTGCTCACTGAAGGCGATGACCAGCAGGATCCGATTGCCGACTCCGCGCGGGCGATCCTCGACGGTCACATTGTGCTGTCTCGCCGACTGGCGGAAGCCGGGCACTATCCGGCTATCGATATTGAAGCGTCGATCAGCCGCG  >  NC_000913/2017447‑2017676

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