Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,021,000 A→G 9.3% Q305R (CAG→CGG)  fliM → flagellar motor switching and energizing component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,021,0000AG 9.3% 119.9 / 3.7 43Q305R (CAG→CGG) fliMflagellar motor switching and energizing component
Reads supporting (aligned to +/- strand):  ref base A (20/19);  new base G (2/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTG  >  NC_000913/2020852‑2021143
                                                                                                                                                    |                                                                                                                                               
atctCACATCTGGAGCTGGTCGCCAACTTTGCGGTTATCTCGCTACGCCTTTCGCAGATTTTAAAACTGATGGCCGGCGACGTCCTGCCGATAGATAAACCCGATCGCATCATCGGCGATGTTGAAGGCGTCCCGTTGCTGACcagtca                                                                                                                                                 <  2:296849/146‑1 (MQ=255)
                  gTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTa                                                                                                                                              <  5:93559/134‑1 (MQ=255)
                  gTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTa                                                                                                                                              >  6:93559/1‑134 (MQ=255)
                      ccAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAg                                                                                                                           >  7:66712/1‑149 (MQ=255)
                      ccAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAg                                                                                                                           >  3:46417/1‑149 (MQ=255)
                       cAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGt                                                                                                                          <  2:168418/149‑1 (MQ=255)
                                 gATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGc                                                                                                                      <  1:34542/143‑1 (MQ=255)
                                 gATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGc                                                                                                                      >  2:34542/1‑143 (MQ=255)
                                   tatCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGa                                                                                                              <  5:2417/149‑1 (MQ=255)
                                    atCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGCCGGCGTCCCGGTGCTGACCAGTCGGTATGGCACCCTCAACGGTCAGTATGCGTTACGGa                                                                                                              >  5:22021/1‑148 (MQ=255)
                                    atCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCGGTATGGCACCCTCAACGGTCAGTATGCGTTACGGa                                                                                                              >  5:17419/1‑148 (MQ=255)
                                      ctcGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGa                                                                                                              >  2:122762/1‑146 (MQ=255)
                                       tcGCTACTCATGTCGCAGATTTTAAAACTGAACCCCGGCGACGTACTGCCTAAAGAAAAACCCTATCGCATCATCGCCCATGTTGACCGCGTCCCGGTGCTGACCAGTCAGTATGGCAACCTCACCGGTCAGTATGCTTTACGGATAg                                                                                                           >  6:211082/1‑148 (MQ=255)
                                             cGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACAttt                                                                                                    <  1:277555/149‑1 (MQ=255)
                                              gCCTGTCGCAGATTTTAAAACTGAACCCCGGCGCCGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGGCAGTATGCGTTACGGATAGAACAttt                                                                                                    >  5:172769/1‑148 (MQ=255)
                                               ccTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTACTGCCGATAGAAAAACCCGATCGCATAATCGCCCATGTTGACGGCGTCCCGGTTCTGACCAGTCAGTATGTCACCCTCAACGGTCAGTATGAGTTACGGATAGAACATTTGa                                                                                                  >  3:206169/1‑149 (MQ=255)
                                                   tCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGACGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACACTCAACGGTCAGAATGCGTTACGGATAGAACATTTGATTaa                                                                                              >  5:149822/1‑149 (MQ=255)
                                                             tAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGAtttt                                                                                     <  1:122762/148‑1 (MQ=255)
                                                                     caCCCCGGCGACGTCCTGCCGCTAGAACAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATtctc                                                                             <  2:142913/147‑1 (MQ=255)
                                                                          cGGCGACGTCCTGCCGATAGAAAAACCCGATCGAATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAAc                                                                        >  3:182423/1‑148 (MQ=255)
                                                                                  tCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAg                                                               <  3:204618/149‑1 (MQ=255)
                                                                                    cTGCCGATAGAAACACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCGGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGcc                                                             <  6:17419/149‑1 (MQ=255)
                                                                                    cTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCGGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGcc                                                             <  6:22021/149‑1 (MQ=255)
                                                                                    cTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACACGATTTTAAATTCTCTGAACGAGGAACAGcc                                                             >  1:21532/1‑149 (MQ=255)
                                                                                           tAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGa                                                      <  6:172769/149‑1 (MQ=255)
                                                                                           tAGAAAAACCCGATCGCATCATCGCCCATGCTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGAACCCTCAACGGTCAGTATGCGTTACGGATAGAACAATTGATTAACCCGATTATAAATTCTCTGAACGAGGAAAATCCAAAAt                                                        >  2:275684/1‑147 (MQ=255)
                                                                                             gAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGt                                                    <  8:77129/149‑1 (MQ=255)
                                                                                              aaaaaCCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGa                                                                                                              >  4:205072/1‑90 (MQ=255)
                                                                                              aaaaaCCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGa                                                                                                              <  3:205072/90‑1 (MQ=255)
                                                                                              aaaaaCCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGt                                                    >  2:174881/1‑148 (MQ=255)
                                                                                                  aCCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACAt                                               <  4:46417/149‑1 (MQ=255)
                                                                                                  aCCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACAt                                               <  5:161046/149‑1 (MQ=255)
                                                                                                     cGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAg                                                                     >  7:264015/1‑124 (MQ=255)
                                                                                                     cGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAg                                                                     <  8:264015/124‑1 (MQ=255)
                                                                                                            atcatcGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCg                                     >  5:204777/1‑149 (MQ=255)
                                                                                                                    ccATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGa                             >  8:91658/1‑149 (MQ=255)
                                                                                                                      aTGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGaca                           <  1:174881/149‑1 (MQ=255)
                                                                                                                        gTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGacaa                          >  1:239222/1‑148 (MQ=255)
                                                                                                                           gACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACaacaac                      >  7:52157/1‑149 (MQ=255)
                                                                                                                              gtcGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGgc                   <  5:211082/147‑1 (MQ=255)
                                                                                                                                 gTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGgcgc                 >  4:143227/1‑148 (MQ=255)
                                                                                                                                    ccGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGCCAACAACGGCGCAAt              <  5:315593/148‑1 (MQ=255)
                                                                                                                                    ccGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGAGTTACGGATAGAAGATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCAGATGACAACAACGGCGCAAt              <  6:149822/148‑1 (MQ=255)
                                                                                                                                       gTGCTGACAAGTCAGTATGGCACCCTGAAAGGTCAGTATGCGTTTCGGATAGAACTTTTTATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGTCCGATGACAACAACGGCGCAATGGAc          <  1:275684/149‑1 (MQ=255)
                                                                                                                                               cagtcagtATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCtgtg  >  8:39699/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCCAAATGAGTGACATGAATAATCCGGCCGATGACAACAACGGCGCAATGGACGATCTGTG  >  NC_000913/2020852‑2021143

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑