Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,019,245 T→C V287A (GTT→GCT)  fliK → flagellar hook‑length control protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,019,2450TC96.9% 108.8 / NA 32V287A (GTT→GCT) fliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (19/12);  minor base A (0/1);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGT  >  NC_000913/2019122‑2019363
                                                                                                                           |                                                                                                                      
aGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCa                                                                                                                    <  2:527246/128‑1 (MQ=255)
aGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCa                                                                                                                    >  1:527246/1‑128 (MQ=255)
  cATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATc                                                                                                             >  2:339820/1‑133 (MQ=255)
  cATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATc                                                                                                             >  1:46069/1‑133 (MQ=255)
       ttCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCAt                                                                                                        >  2:553392/1‑133 (MQ=255)
       ttCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCAt                                                                                                        >  2:508171/1‑133 (MQ=255)
              ttCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTAcgcg                                                                                                 <  1:553392/133‑1 (MQ=255)
              ttCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTAcgcg                                                                                                 <  1:508171/133‑1 (MQ=255)
               tCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCGGATGGCTGCACCGCATCAGCAAGTACgcgc                                                                                                >  1:276891/1‑133 (MQ=255)
               tCACCCGCAAGGGGCCACCAAGTGCAGAGTAGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACgcgc                                                                                                <  1:200163/133‑1 (MQ=255)
                cACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGcgcc                                                                                               <  2:359806/133‑1 (MQ=255)
                                caaAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAgcgc                                                                               >  2:583917/1‑133 (MQ=255)
                                               cGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGcgca                                                                >  1:627066/1‑133 (MQ=255)
                                                    gCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAg                                                           >  2:411027/1‑133 (MQ=255)
                                                    gCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAg                                                           >  1:139924/1‑133 (MQ=255)
                                                                  tAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGGGCATGCAGCTGGCAGAAAGTgg                                             >  1:232930/1‑133 (MQ=255)
                                                                     gTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCAt                                          >  2:65023/1‑133 (MQ=255)
                                                                     gTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCAt                                          <  2:139924/133‑1 (MQ=255)
                                                                     gTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCAt                                          <  1:339820/133‑1 (MQ=255)
                                                                     gTGAAGTGCAAATCTCCCTCAAAGCGGATGATTACCAGGCGAAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCAt                                          <  2:232930/133‑1 (MQ=255)
                                                                          gTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGt                                     >  2:326934/1‑133 (MQ=255)
                                                                                aTCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGc                               <  1:411027/133‑1 (MQ=255)
                                                                                     ccTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGc                          >  1:528191/1‑133 (MQ=255)
                                                                                          aaGTGGATGATAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACAt                     >  2:430712/1‑133 (MQ=255)
                                                                                                 tgatAACCTGGCGCAGATCCAGGTGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGc              <  1:211595/133‑1 (MQ=255)
                                                                                                 tgatAACCAGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGc              <  2:528191/133‑1 (MQ=255)
                                                                                                      aCCAGGCGTAAATCCTGACGGCTTCACCGCCTCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAg         <  1:163752/133‑1 (MQ=255)
                                                                                                         aGGCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCtt      >  1:519822/1‑133 (MQ=255)
                                                                                                          ggCGCAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCttt     >  1:23524/1‑133 (MQ=255)
                                                                                                            cgcAAATCCAGATGGATTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAg   <  2:276891/133‑1 (MQ=255)
                                                                                                             gcAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCATCATCAGTGGCGAATGCTTTAGt  >  1:207663/1‑133 (MQ=255)
                                                                                                             gcAAATCCAGATGGCTTCACCGCATCAGCATGTACGCGCCGCCCTAGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGt  >  2:447989/1‑133 (MQ=255)
                                                                                                                           |                                                                                                                      
AGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCATCAGCATGTACGCGCCGCCCTGGAAGCAGCGCTGCCGGTACTGCGCACGCAGCTGGCCGAAAGTGGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGT  >  NC_000913/2019122‑2019363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑