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

GCACTCATCGCCGTAGCCTGTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTTATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGAATGCGACTGCAGCCC  >  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)
            gTAGCCTGTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑aa                                                                                                                                      <  1:314004/151‑1 (MQ=255)
                 cTGTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑aaaaaa                                                                                                                                  >  1:492361/1‑150 (MQ=255)
                   gTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAAc                                                                                                                                <  2:720472/150‑1 (MQ=255)
                   gTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACg                                                                                                                               >  1:708829/1‑151 (MQ=255)
                     aaCTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGa                                                                                                                              <  1:460468/150‑1 (MQ=255)
                         ggCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAg                                                                                                                           <  2:666498/149‑1 (MQ=255)
                          gCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGcac                                                                                                                        >  2:16428/1‑151 (MQ=255)
                              aTAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACaa                                                                                                                      <  2:417149/149‑1 (MQ=255)
                              aTAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAaat                                                                                                                    <  2:64793/151‑1 (MQ=255)
                                 aCCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAaataat                                                                                                                 >  2:840432/1‑151 (MQ=255)
                                        cAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCg                                                                                                          >  2:672926/1‑151 (MQ=255)
                                        cAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCg                                                                                                          <  1:592049/151‑1 (MQ=255)
                                           tCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTc                                                                                                        >  1:373341/1‑150 (MQ=255)
                                            ccccTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCa                                                                                                       <  1:190267/150‑1 (MQ=255)
                                                   aTCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCg                                                                                               <  1:384893/151‑1 (MQ=255)
                                                      gCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGt                                                                                            <  1:803764/151‑1 (MQ=255)
                                                         gACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGttt                                                                                          <  1:922567/150‑1 (MQ=255)
                                                         gACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGttt                                                                                          <  1:269630/150‑1 (MQ=255)
                                                         gACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGttt                                                                                          <  1:274796/150‑1 (MQ=255)
                                                          aCATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGc                                                                                        >  1:99160/1‑151 (MQ=255)
                                                           cATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCt                                                                                       <  1:282250/151‑1 (MQ=255)
                                                           cATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCt                                                                                       >  2:554967/1‑151 (MQ=255)
                                                                 cGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGAc                                                                                    <  2:165702/148‑1 (MQ=255)
                                                                      cccGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCggg                                                                               <  2:745301/148‑1 (MQ=255)
                                                                       ccGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGgagtg                                                                           <  2:815717/151‑1 (MQ=255)
                                                                          gATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGgagtgagt                                                                        >  1:760610/1‑151 (MQ=255)
                                                                           aTAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTc                                                                       <  1:703448/151‑1 (MQ=255)
                                                                                        gCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTg                                                          <  2:747468/151‑1 (MQ=255)
                                                                                           aaGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGATTGAGTCTTGTTCCACTTTGcc                                                        >  1:227683/1‑150 (MQ=255)
                                                                                                 aCCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAAc                                                  <  1:397093/150‑1 (MQ=255)
                                                                                                  ccAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGc                                                <  2:778845/151‑1 (MQ=255)
                                                                                                    aGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCg                                              <  1:518631/151‑1 (MQ=255)
                                                                                                      tAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAGAAAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGt                                             >  1:843092/1‑151 (MQ=255)
                                                                                                             aTAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATc                                     >  2:497312/1‑151 (MQ=255)
                                                                                                               aaGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATc                                     >  1:866288/1‑149 (MQ=255)
                                                                                                                  aCAAAGGCGCTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCa                                <  2:302770/151‑1 (MQ=255)
                                                                                                                        gcgcTAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTTTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGt                          <  2:856984/151‑1 (MQ=255)
                                                                                                                            tAAATAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTgcgc                      >  2:100403/1‑151 (MQ=255)
                                                                                                                             aaaTAGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTgcg                       <  2:704321/149‑1 (MQ=255)
                                                                                                                             aaaTAGCAACAAAAAAACGGGTTAATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGcgcg                     >  1:36031/1‑151 (MQ=255)
                                                                                                                                 aGCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGa                 <  1:143039/151‑1 (MQ=255)
                                                                                                                                  gCAACAAAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGaa                <  1:440771/151‑1 (MQ=255)
                                                                                                                                      caaAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGAATGCg            >  1:657641/1‑151 (MQ=255)
                                                                                                                                      caaAAAAACGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGAATGCg            <  2:795941/151‑1 (MQ=255)
                                                                                                                                            aaCGGGTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGAATGCGACTGc       >  1:760043/1‑150 (MQ=255)
                                                                                                                                                ggTTGATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGAATGCGACTGCAGccc  >  2:480393/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCACTCATCGCCGTAGCCTGTAACTGGCTGATAACCCCTTCAATCCCCTGTATCGCTGACATTCTCGTCTCCCGGATAATTTCTGGTAGCAAAGCCTACCAGTAAGTCAATAAGACAAAGGCGCTAAATAGCAACAAAAAAACGGGTTTATTGGCGGATAG‑AAAAAAACGAAAGCACAAATAATGGGAGCGTCAATTTTTCGAGTTTGCTGACCCGGGAGTGAGTCTTGTTCCACTTTGCCAATAACGCCGTCCATAATCAGCCACGAGGTGCGCGATGAATGCGACTGCAGCCC  >  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: ‑