Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,137,212	C→T	R313R (CGC→CGT)	flgI →	flagellar P‑ring protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,137,212	0	C	T	95.8%	76.4 / NA	25	R313R (CGC→CGT)	flgI	flagellar P‑ring protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/15);  minor base G (1/0);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGCCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTGATGTCCATACTGCAATCAATGCAAAGTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCAATCACTCAACGAACTAAAGGCGAAAGCGGGCGAA  >  NC_000913/1136968‑1137451                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   aGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCTGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAgc                                                                                                                                                                                                                                              <  1:62670/251‑1 (MQ=255)  ggTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAgcg                                                                                                                                                                                                                                             >  2:7631/1‑251 (MQ=255)         ttCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCCGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTAGGTGGTGATGAATCGCGAAGTGACCCTCGACCGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGGCAGAGCGGCGGtgg                                                                                                                                                                                                                                      >  1:194858/1‑249 (MQ=255)         ttCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGATAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATAAACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGttt                                                                                                                                                                                                                                      >  1:254941/1‑250 (MQ=255)                           aagatGCAGGTTAATGTCACCCCGCAGGACGCGAAAGTAGTGATTAACTCGCGCGCCGGTTCGGTGGTGATGAATCGCGAAGGGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAAGCTCGCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGGGGACAGACTGTGGTTACTCCACAAACGCAGATCGATGTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTc                                                                                                                                                                                                                  <  2:134686/248‑1 (MQ=255)                                    gTTAATGTAACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGAGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGGAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑AGCTGCAAA‑GCGTACGTTCCAGCG‑CCAg                                                                                                                                                                                                        <  2:179388/251‑1 (MQ=255)                                                                        aaCTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTgg                                                                                                                                                                     <  2:120709/250‑1 (MQ=255)                                                                                          gtggtgATGAATCGCGAAGTGACCCTCGAGATCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGGCAGGCCAATGTCAGCCAGCCAGATACGCCGTTGGGTGGTTGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCTGTG‑CGCTGCAAT‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGTTGCTCGCGCTTAATGCGCTTGGTGCTACGCCTATGGAttt                                                                                                                                                  <  2:174140/251‑3 (MQ=255)                                                                                                              gACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTGATGTCCATACTGcaatca                                                                                                                               <  1:135725/250‑1 (MQ=255)                                                                                                                     gACAGCTGCGCGGTAGCGCAGGGGAATCTCTCCGTAAAAGTTAATCGTCCGGCCAATGTAAGCAAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTTTCG‑TTCAAACACG‑ACGTTCCCGCG‑CCCGCCCCAATAACGCTGTGCGGGCGCCCCATGCGGCGGGCGCTAAGACGGTGGATCTGGTGTACCTACCGcaatacatggaaa                                                                                                                        >  1:24121/1‑241 (MQ=255)                                                                                                                     gACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCCGCG‑CCAGCCTCACAAAAGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGGTCTGATGTCCATACTGCAATCAATGCaa                                                                                                                         >  1:71305/1‑249 (MQ=255)                                                                                                                       cAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGACAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTGATGTCCATACTGCAATCAATGCAAAGTg                                                                                                                     <  2:207249/251‑1 (MQ=255)                                                                                                                                      gcAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑CGCTGCCAA‑ACGAACGTTCCAGCGGCCGCCCTCA‑TAACGAGGTAGGCGCGCCCAAAGCGCTGGGCGCTACGCCCATGGATCTGATGGCCATACTTCAATTAAAGCaaaaggccgaacggaagagc                                                                                                       >  1:262576/1‑233 (MQ=255)                                                                                                                                      gcAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCGCAATGCGCTGGGCGCTACGCCGATGGATCTGATGTCCATACTGCAATCAATGCAAAGTGc                                                                                                                    <  2:262576/237‑1 (MQ=255)                                                                                                                                                    cAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACCCCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTGATGTCCATACTGCAATCAATGCAAAGTGCGGGATGTCTGCGGGCAAAACTGGaa                                                                                           <  2:5270/248‑1 (MQ=255)                                                                                                                                                                   gTCAGGCCAATGTCAGCCAGCCAGATACACCGTGTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTGATGTCCATACTGCAAGCAATGCAAAGTGCGGGATGTCTGCGGGAAAAACTGGAAATCATCTGATGATCAGCg                                                                         <  2:293031/251‑1 (MQ=255)                                                                                                                                                                     cAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCGACGCCGATGGATCTGATGTCCATACGGCAATCAATGCAAAGTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGAc                                                                       <  2:274278/251‑1 (MQ=255)                                                                                                                                                                       ggCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTGATGTCCATACTGCAATCAATGCAAAGTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAg                                                                     <  2:252058/251‑1 (MQ=255)                                                                                                                                                                                        cagATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑TTCCGCAAA‑GCGTAAGGTCCAGCG‑CCCGCCCCAAAAAAGGGGTGGGCGCGCCCAATGCGCTGGGCGCCACGCCCGTGGAACTGATGTCCCTAATGCAATCAAAGCAAAGTGCGGGGAGTCCGAGGGCCAAACTGGAAAACATCTGATGATACGCGACCGCAAACTACTGGccaggg                                                    >  1:60152/1‑246 (MQ=255)                                                                                                                                                                                         agATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑TGCTGCCAA‑ACGTACGTTCCCGCG‑CCAGCCCCCATAACGTGGGGCGCGCGCGCAATGCGCCGGGCGCTACGCCCATGGGTCTGATGTCCCCACTTCAATCAATGCAAAGTGCGGGGTGTCTGCGGGCAAACCTGGAAAACATCTGATGATCCGCGACAGCCAACTACTGGCCAGTGc                                                   >  2:184201/1‑251 (MQ=255)                                                                                                                                                                                                      ggtggtggACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTGATGTCCATACTGCAATCAATGCAAAGTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATgc                                       <  1:7631/250‑1 (MQ=255)                                                                                                                                                                                                                aGACTGTGGTGACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTGATGTCCATACTGCAAGCAATGCAAAGTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCAATCACTCa                            <  2:145720/251‑1 (MQ=255)                                                                                                                                                                                                                          tACTCCACAAACGCAGATCGATTTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTGATGTCCATACTGCAATCAATGCAAAGTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCAATCACTCAACGAACTaa                   <  1:36664/250‑1 (MQ=255)                                                                                                                                                                                                                                          aTCGATTTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTGATGTCCATACTGCAATCAATGCAAAGTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCAATCACTCAACGAACTAAAGGCGAAAGCGGGCGaa  <  2:202466/251‑1 (MQ=255)                                                                                                                                                                                                                                              aTTTACGTCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTGATGTCCATACTGCAATCAATGCAAAGTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCAATCACTCAACGAACTAAAGGCGAAAACGGGCGa   >  1:3050/1‑246 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   AGGTCCGCTTCCTTGCCGATATTCAGAATATGCAGGTTAATGTCACCCCGCAGGACGCTAAAGTAGTGATTAACTCGCGCACCGGTTCGGTGGTGATGAATCGCGAAGTGACCCTCGACAGCTGCGCGGTAGCGCAGGGGAATCTCTCAGTAACAGTTAATCGTCAGGCCAATGTCAGCCAGCCAGATACACCGTTTGGTGGTGGACAGACTGTGGTTACTCCACAAACGCAGATCGATTTACGCCAGAGCGGCGGTT‑CGCTGCAAA‑GCGTACGTTCCAGCG‑CCAGCCTCAATAACGTGGTGCGCGCGCTCAATGCGCTGGGCGCTACGCCGATGGATCTGATGTCCATACTGCAATCAATGCAAAGTGCGGGATGTCTGCGGGCAAAACTGGAAATCATCTGATGATCAGCGACAGCAAACTACTGGCAAGTGCGGCCTGGGATGCGCAATCACTCAACGAACTAAAGGCGAAAGCGGGCGAA  >  NC_000913/1136968‑1137451

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