Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,013,101 T→G intergenic (‑87/‑128) fliE ← / → fliF flagellar basal‑body component/flagellar basal‑body MS‑ring and collar protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,013,1010TG98.0% 182.3 / ‑4.0 49intergenic (‑87/‑128)fliE/fliFflagellar basal‑body component/flagellar basal‑body MS‑ring and collar protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (20/28);  minor base A (1/0);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.19e-01

GCACTCATCGCCGTAGCCTGTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTTATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGAATGCGACTGCAGCCC  >  NC_000913/2012953‑2013247
                                                                                                                                                    |                                                                                                                                                   
gcacTCATCGCCGTAGCCTGTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGAt                                                                                                                                                   >  2:33489/1‑151 (MQ=255)
     cATCGCCGTAGCCTGTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCg                                                                                                                                              >  2:115376/1‑151 (MQ=255)
         gCCGTAGCCTGTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGAt                                                                                                                                           >  2:397084/1‑150 (MQ=255)
            gTAGCCTGTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGaa                                                                                                                                      <  1:314004/151‑1 (MQ=255)
                 cTGTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGaaaaaa                                                                                                                                  >  1:492361/1‑150 (MQ=255)
                   gTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAAc                                                                                                                                <  2:720472/150‑1 (MQ=255)
                   gTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACg                                                                                                                               >  1:708829/1‑151 (MQ=255)
                     aaCTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGa                                                                                                                              <  1:460468/150‑1 (MQ=255)
                         ggCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAg                                                                                                                           <  2:666498/149‑1 (MQ=255)
                          gCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGcac                                                                                                                        >  2:16428/1‑151 (MQ=255)
                              aTAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACaa                                                                                                                      <  2:417149/149‑1 (MQ=255)
                              aTAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAaat                                                                                                                    <  2:64793/151‑1 (MQ=255)
                                 aCCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAaataat                                                                                                                 >  2:840432/1‑151 (MQ=255)
                                        cAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCg                                                                                                          >  2:672926/1‑151 (MQ=255)
                                        cAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCg                                                                                                          <  1:592049/151‑1 (MQ=255)
                                           tCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTc                                                                                                        >  1:373341/1‑150 (MQ=255)
                                            ccccTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCa                                                                                                       <  1:190267/150‑1 (MQ=255)
                                                   aTCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCg                                                                                               <  1:384893/151‑1 (MQ=255)
                                                      gCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGt                                                                                            <  1:803764/151‑1 (MQ=255)
                                                         gACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGttt                                                                                          <  1:922567/150‑1 (MQ=255)
                                                         gACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGttt                                                                                          <  1:269630/150‑1 (MQ=255)
                                                         gACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGttt                                                                                          <  1:274796/150‑1 (MQ=255)
                                                          aCATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGc                                                                                        >  1:99160/1‑151 (MQ=255)
                                                           cATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCt                                                                                       <  1:282250/151‑1 (MQ=255)
                                                           cATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCt                                                                                       >  2:554967/1‑151 (MQ=255)
                                                                 cGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGAc                                                                                    <  2:165702/148‑1 (MQ=255)
                                                                      cccGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCggg                                                                               <  2:745301/148‑1 (MQ=255)
                                                                       ccGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGgagtg                                                                           <  2:815717/151‑1 (MQ=255)
                                                                          gATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGgagtgagt                                                                        >  1:760610/1‑151 (MQ=255)
                                                                           aTAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTc                                                                       <  1:703448/151‑1 (MQ=255)
                                                                                        gCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTg                                                          <  2:747468/151‑1 (MQ=255)
                                                                                           aaGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGATTGAGTCTTGTTCCACTTTGcc                                                        >  1:227683/1‑150 (MQ=255)
                                                                                                 aCCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAAc                                                  <  1:397093/150‑1 (MQ=255)
                                                                                                  ccAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGc                                                <  2:778845/151‑1 (MQ=255)
                                                                                                    aGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCg                                              <  1:518631/151‑1 (MQ=255)
                                                                                                      tAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGt                                             >  1:843092/1‑151 (MQ=255)
                                                                                                             aTAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATc                                     >  2:497312/1‑151 (MQ=255)
                                                                                                               aaGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATc                                     >  1:866288/1‑149 (MQ=255)
                                                                                                                  aCAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCa                                <  2:302770/151‑1 (MQ=255)
                                                                                                                        gcgcTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTTTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGt                          <  2:856984/151‑1 (MQ=255)
                                                                                                                            tAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTgcgc                      >  2:100403/1‑151 (MQ=255)
                                                                                                                             aaaTAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTgcg                       <  2:704321/149‑1 (MQ=255)
                                                                                                                             aaaTAGCAACAAAAAAACGGGTTAATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGcgcg                     >  1:36031/1‑151 (MQ=255)
                                                                                                                                 aGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGa                 <  1:143039/151‑1 (MQ=255)
                                                                                                                                  gCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGaa                <  1:440771/151‑1 (MQ=255)
                                                                                                                                      caaAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGAATGCg            >  1:657641/1‑151 (MQ=255)
                                                                                                                                      caaAAAAACGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGAATGCg            <  2:795941/151‑1 (MQ=255)
                                                                                                                                            aaCGGGTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGAATGCGACTGc       >  1:760043/1‑150 (MQ=255)
                                                                                                                                                ggTTGATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGAATGCGACTGCAGccc  >  2:480393/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCACTCATCGCCGTAGCCTGTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTTATTGGCGGATAGAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGAATGCGACTGCAGCCC  >  NC_000913/2012953‑2013247

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