Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,015,192 A→G 8.3% S105G (AGC→GGC)  fliG → flagellar motor switching and energizing component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,015,1920AG 8.3% 119.5 / 3.6 49S105G (AGC→GGC) fliGflagellar motor switching and energizing component
Reads supporting (aligned to +/- strand):  ref base A (26/18);  new base G (2/2);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

CTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTC  >  NC_000913/2015045‑2015336
                                                                                                                                                   |                                                                                                                                                
cTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAg                                                                                                                                                 <  2:935747/149‑1 (MQ=255)
           tGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCTGGCGGTAtcgcga                                                                                                                                      <  1:360639/149‑6 (MQ=255)
            gAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCTGGCGGTAtcgcgag                                                                                                                                     <  2:622672/149‑7 (MQ=255)
            gAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAAc                                                                                                                                     >  1:829587/1‑149 (MQ=255)
              gCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGGCAGGCGGCTGGAAGATATGCTCGTAACTCGAGATTCAGCCAGCGGGATTGAAACGc                                                                                                                                   >  2:695035/1‑149 (MQ=255)
                 agaagCTGAACAGTTTGCCGCACTGAATTTCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCTGCATGCTGGATGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCa                                                                                                                                >  2:128351/1‑149 (MQ=255)
                 agaagCTGAACAGTTTGCCGCACGTAATATCAACGCGAACGATTATCTGCGCTAGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCa                                                                                                                                <  2:575784/149‑1 (MQ=255)
                    agCTGAACAGTTTGGCGCACTGAATATCAACGGCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTGTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACt                                                                                                                             <  2:79338/149‑1 (MQ=255)
                       tGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTa                                                                                                                          >  2:11103/1‑149 (MQ=255)
                        taaCAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCAAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTAt                                                                                                                         <  1:907328/148‑1 (MQ=255)
                           cAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGa                                                                                                                      <  1:328487/149‑1 (MQ=255)
                               ttGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCca                                                                                                                  <  1:289094/149‑1 (MQ=255)
                                tGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCcac                                                                                                                 <  1:412202/149‑1 (MQ=255)
                                 gCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCCGCGGTATTGAACCGCTCAACTTTATGGAGCcaca                                                                                                                >  2:803819/1‑149 (MQ=255)
                                   cGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACaga                                                                                                              >  2:836257/1‑149 (MQ=255)
                                      aCTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAGGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAgcg                                                                                                           >  2:360874/1‑149 (MQ=255)
                                              tCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGAt                                                                                                   <  1:27682/149‑1 (MQ=255)
                                                  ggCCAACGATTATCTGCGCTCTGTATTGGTCAAAGCTCTGGGTGAAGAACGGGCCGCGAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGa                                                                                               <  1:743385/148‑1 (MQ=255)
                                                  cGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGGGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGa                                                                                               <  1:71016/149‑1 (MQ=255)
                                                      aaCGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTcg                                                                                           >  2:524489/1‑149 (MQ=255)
                                                         gATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGa                                                                                        >  2:551376/1‑149 (MQ=255)
                                                             aTCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAg                                                                                    >  1:895411/1‑149 (MQ=255)
                                                                 gcgcTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGAGCCGCCGCTCTGATTCGCGATGAGCATc                                                                                >  1:756485/1‑149 (MQ=255)
                                                                   gcTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCg                                                                              >  1:234494/1‑149 (MQ=255)
                                                                        gTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGCGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAAt                                                                         >  1:30680/1‑149 (MQ=255)
                                                                               tCACAGCTCTGGGCGAATAACGTGCCGCCAGCCTGCTTGAATATCTTCTCGCAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGcc                                                                  <  2:971436/149‑1 (MQ=255)
                                                                                 aaaGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGccac                                                                >  1:865428/1‑149 (MQ=255)
                                                                                 aaaGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCAGATCTGATTCGCGATGACCAGCCGCAAATTATCGccac                                                                >  1:989058/1‑149 (MQ=255)
                                                                                        tGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTg                                                         >  1:47573/1‑149 (MQ=255)
                                                                                          ggTGAGGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGt                                                       <  1:360874/149‑1 (MQ=255)
                                                                                           gTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTg                                                      >  2:953796/1‑149 (MQ=255)
                                                                                              aagaaCGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAATTTATCGCCACCATTCTGGTGCAt                                                   <  1:551376/149‑1 (MQ=255)
                                                                                               tgaaCGTGCCGCGAGCCTGCTGGAAGATATTCTCGAAACTCGCTATACCGCCAGCGGTATTGAAACGCTCATCTTTATGGAGACACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATc                                                  <  1:968339/148‑1 (MQ=255)
                                                                                                          ccAGCCTGCTGGAAGATATTCTCGAAACTCGCGATAGCGCCAGCGGTTTTGAACCGCTCAACTTTATGGAGACACAGAGCGCCGCCGATCTGATTCGCGATGCGCATCCGCAAATTATCGCCACAATACTGGTGCATCTGAAGCGCGcc                                       >  2:205310/1‑149 (MQ=255)
                                                                                                          ccAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCAGCTGAAGCGCGcc                                       >  2:665372/1‑149 (MQ=255)
                                                                                                             gcatgctgGAAGATATTCTCTAAACTCGCGATACCGACATAGGTATTGAAACTCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAATCGCGCCCaa                                    <  1:403422/146‑1 (MQ=255)
                                                                                                              ccTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAg                                   >  2:629832/1‑149 (MQ=255)
                                                                                                              ccTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCAACAGCGCGCCGCAGATCTGATTCGCGATGAGCATCCGCAAATAATCGCCACCATTCTGGTGCATCTGCAGCGCGCCCAAg                                   >  2:674423/1‑149 (MQ=255)
                                                                                                                  ctgGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAgccgc                               >  1:390252/1‑149 (MQ=255)
                                                                                                                       aGATATTCTCGAATCTCGCGATACCGCCAGCGGTATTGAAACGCACAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCGGAAGCGCGCCCAAGCCGCCGata                          >  1:827892/1‑149 (MQ=255)
                                                                                                                           atTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCt                      >  1:130238/1‑149 (MQ=255)
                                                                                                                            ttCTCGACACTCGCGATACCGCCAGCGGTATGGAAACGCTCAACTTTAGGGAGCCACAGAGCGCCGCCGATCTGAATCGCGATGAGCCTCCGCAAATTCTCGCCACCATTCTGGTGCATCTGCAGCGCGCCCAAGCCGCCGATATTCTg                     <  1:465594/149‑1 (MQ=255)
                                                                                                                                cGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCgt                 >  2:968676/1‑149 (MQ=255)
                                                                                                                                    aCTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGACACCATTATGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGAgttgtt             >  2:33203/1‑149 (MQ=255)
                                                                                                                                      tCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCg           <  2:521475/149‑1 (MQ=255)
                                                                                                                                       ttcaaTACCGCTGGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGGATCTGAAGCGCGCCCAAGCCGCCGATTTTCTGGCGTTGTTCGa          >  1:417568/5‑149 (MQ=255)
                                                                                                                                        gcgATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGAt         <  1:1010724/149‑1 (MQ=255)
                                                                                                                                         caaTACCGCTGGCGGTATTGAAACGCTCAACTTTATGGATCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATg        >  1:915963/3‑149 (MQ=255)
                                                                                                                                               cGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTc  >  1:515393/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTGGCGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGAACGTC  >  NC_000913/2015045‑2015336

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: ‑