Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,020,943 T→C 7.7% I286T (ATA→ACA)  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,020,9430TC 7.7% 151.4 / 3.8 52I286T (ATA→ACA) fliMflagellar motor switching and energizing component
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base C (2/2);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

GCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCC  >  NC_000913/2020803‑2021085
                                                                                                                                            |                                                                                                                                              
gcgTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAAc                                                                                                                                        >  6:96259/1‑149 (MQ=255)
gcgTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAAc                                                                                                                                        <  5:93518/149‑1 (MQ=255)
     aTGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGAt                                                                                                                                   >  5:14317/1‑149 (MQ=255)
      tGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATc                                                                                                                                  <  8:101750/149‑1 (MQ=255)
         agatcagatcTGGCGCGATAACCTGGTGCGCCAGGTGCAGCTTTCACAGCTGGAGCTGGTCGGCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCGGATAGATAAACGCGATCGc                                                                                                                                <  6:38868/148‑1 (MQ=255)
           aTCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGcat                                                                                                                              >  4:30562/1‑148 (MQ=255)
           aTCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGc                                                                                                                                <  5:100128/146‑1 (MQ=255)
           aTCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGc                                                                                                                                >  6:100128/1‑146 (MQ=255)
               gAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCg                                                                                                                         >  4:33104/1‑149 (MQ=255)
               gAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCg                                                                                                                         <  4:28625/149‑1 (MQ=255)
                aaCTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGc                                                                                                                        >  3:8554/1‑149 (MQ=255)
                aaCTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGc                                                                                                                        <  5:198415/149‑1 (MQ=255)
                     gcgcGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATg                                                                                                                   >  2:160231/1‑149 (MQ=255)
                       gcgATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGtt                                                                                                                 >  7:22149/1‑149 (MQ=255)
                       gcgATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGtt                                                                                                                 >  2:198218/1‑149 (MQ=255)
                         gATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGa                                                                                                               >  8:166626/1‑149 (MQ=255)
                                     cgcCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGt                                                                                                   <  3:33104/149‑1 (MQ=255)
                                        cAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCt                                                                                                <  6:14317/149‑1 (MQ=255)
                                              cagcaTTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACca                                                                                           <  5:96259/148‑1 (MQ=255)
                                               agcaTTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACcagt                                                                                         <  1:40185/149‑1 (MQ=255)
                                                gcaTTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACcagt                                                                                         >  7:27441/1‑148 (MQ=255)
                                                            ggagctggTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTa                                                                                    <  7:88020/141‑1 (MQ=255)
                                                            ggagctggTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTa                                                                                    >  8:88020/1‑141 (MQ=255)
                                                             gagctggTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCAcc                                                                             <  8:22149/147‑1 (MQ=255)
                                                                ctggTCGCCAACTTTACCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCaa                                                                        >  3:20028/1‑149 (MQ=255)
                                                                 tggTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAAc                                                                       >  2:41918/1‑149 (MQ=255)
                                                                      gCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGt                                                                    >  5:69476/1‑147 (MQ=255)
                                                                      gCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTc                                                                   <  5:135137/148‑1 (MQ=255)
                                                                      gCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCa                                                                  <  1:160231/149‑1 (MQ=255)
                                                                       ccAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGACAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACcagtca                                                                                       >  1:10963/1‑127 (MQ=255)
                                                                       ccAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGACAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACcagtca                                                                                       >  1:7010/1‑127 (MQ=255)
                                                                       ccAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGACAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACcagtca                                                                                       <  2:10963/127‑1 (MQ=255)
                                                                       ccAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGACAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACcagtca                                                                                       <  2:7010/127‑1 (MQ=255)
                                                                          aCTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTAt                                                              <  5:91530/149‑1 (MQ=255)
                                                                                 cGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTAc                                                       >  3:215923/1‑149 (MQ=255)
                                                                                      tctcGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGa                                                    >  8:237142/1‑147 (MQ=255)
                                                                                          gCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGaa                                               <  4:8554/148‑1 (MQ=255)
                                                                                            tACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACAt                                            >  8:164420/1‑149 (MQ=255)
                                                                                                     cGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAAc                                   <  6:69476/149‑1 (MQ=255)
                                                                                                       cAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAAccc                                 >  6:177508/1‑149 (MQ=255)
                                                                                                           ttttAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGAtt                             >  1:131615/1‑149 (MQ=255)
                                                                                                            tttAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCCCCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGAttt                            >  5:122032/1‑149 (MQ=255)
                                                                                                               aaaaCTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTa                          >  3:81405/1‑148 (MQ=255)
                                                                                                                aaaCTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTaa                         <  1:198218/148‑1 (MQ=255)
                                                                                                                 aaCTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAAc                                                                       >  6:224708/1‑101 (MQ=255)
                                                                                                                 aaCTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAAc                                                                       <  5:224708/101‑1 (MQ=255)
                                                                                                                      aaCCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATtctc                   <  3:30562/148‑1 (MQ=255)
                                                                                                                      aaCCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATtct                    <  4:20028/147‑1 (MQ=255)
                                                                                                                      aaCCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTctct                  <  1:41918/149‑1 (MQ=255)
                                                                                                                           cGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAAc              >  7:181386/1‑148 (MQ=255)
                                                                                                                             gCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAg           <  7:166626/149‑1 (MQ=255)
                                                                                                                                        ccGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGccc  <  7:164420/147‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                              
GCGTAATGAAGATCAGAACTGGCGCGATAACCTGGTGCGCCAGGTGCAGCATTCACAGCTGGAGCTGGTCGCCAACTTTGCCGATATCTCGCTACGCCTGTCGCAGATTTTAAAACTGAACCCCGGCGACGTCCTGCCGATAGAAAAACCCGATCGCATCATCGCCCATGTTGACGGCGTCCCGGTGCTGACCAGTCAGTATGGCACCCTCAACGGTCAGTATGCGTTACGGATAGAACATTTGATTAACCCGATTTTAAATTCTCTGAACGAGGAACAGCCC  >  NC_000913/2020803‑2021085

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