Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,018,530	(C)6→7	coding (145/1128 nt)	fliK →	flagellar hook‑length control protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,018,530	1	.	C	96.3%	87.5 / ‑2.5	27	coding (145/1128 nt)	fliK	flagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (11/15);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

GGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCC‑AGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAACCACCGCTCAGACGATGGCGTTGGCTGCGGTGGCTGACAAAAATACGACAAAAGACGAAAAAGCGGATGATCTGAATGAAGACGTCACCGCAAGCCTGAG  >  NC_000913/2018295‑2018780                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                           ggCAACCCAGGGCGGTACTGCTTGAAAAAACCCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGGCAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCa                                                                                                                                                                                                                                              <  1:135296/251‑1 (MQ=255)    aaTCCACGGCGGCACAGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGGGGCacc                                                                                                                                                                                                                                            >  2:85516/1‑248 (MQ=255)                              aaCCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGcc                                                                                                                                                                                                                 <  1:1304/251‑1 (MQ=255)                                ccGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGc                                                                                                                                                                                                               <  1:130983/251‑1 (MQ=255)                                  gCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCCCGACAAAAGGCGAGCCGCTg                                                                                                                                                                                                             >  2:139480/1‑251 (MQ=255)                                      cGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGAt                                                                                                                                                                                                           <  1:78170/249‑1 (MQ=255)                                            gAAAAAGATGGATGAGTTCGCCCCGCGCGCCGCCATGAGGAACCCTGAATGATTCGCTTAGCGCCCTTGCTTACCGACGACGTTGACACCCCCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGAACAGCGat                                                                                                                                                                                                    >  2:257860/1‑250 (MQ=255)                                                     ggatgAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTcc                                                                                                                                                                                           <  1:267490/250‑1 (MQ=255)                                                          aGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGAcgcg                                                                                                                                                                                     >  2:76063/1‑251 (MQ=255)                                                                           gcaaggagcaacccacaaggATTCGCTTCGCGCCCTTGATTACCGCCGCCGTTTACACCGCCTCATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTGGTTTAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCATCAGATATGCCCACGACAAAATGCTAGCCGCTGATCAGCGATATTGTTTCCGACACGCAACAAGCTAATTTACt                                                                                                                                                                    <  2:15612/232‑1 (MQ=255)                                                                                     aaCCTGACTGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGAACAAGATTTTCACGCGTTGTTGAGAGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACCAGCTACATTACTGATCACtgtg                                                                                                                                                          >  2:173924/1‑251 (MQ=255)                                                                                                 tCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACAcc                                                                                                                                               <  1:178630/250‑1 (MQ=255)                                                                                                    cTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCt                                                                                                                                           >  1:188959/1‑251 (MQ=255)                                                                                                       aGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGtc                                                                                                                                        <  1:192634/251‑1 (MQ=255)                                                                                                                  ttACCGCCGACGTTGACACCACCACATTGCCTGGCGGAAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCACCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAacgacg                                                                                                                              >  1:284438/1‑250 (MQ=255)                                                                                                                           aCGTTGACACCACCACATTGCGTGGCGGCAAAGCCAGCGATGCTGCACACGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCATCAGCTAAGCCCACGACAACAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGACCAATCcac                                                                                                                    >  2:268041/1‑251 (MQ=255)                                                                                                                               tGACACCACCCCATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCACCAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATc                                                                                                                 <  1:101472/250‑1 (MQ=255)                                                                                                                                     caccacATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCg                                                                                                          <  1:76063/251‑1 (MQ=255)                                                                                                                                             ggccgggcggcAAACCCAGCGATGCGGCACAAGATTTTCGCGCGTGGTGGAGCGAAGCATTAGCAGGCGAGACAACTCCCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTTTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAaccac                                                                                                  <  2:203472/246‑1 (MQ=255)                                                                                                                                                     ggcAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGa                                                                                                                                                                                         <  1:287126/156‑1 (MQ=255)                                                                                                                                                     ggcAAAGCCAGCGATGCTGCACAAGATTTTATCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGa                                                                                                                                                                                         >  2:287126/1‑156 (MQ=255)                                                                                                                                                     ggcAAAGCCAGCGATGCTGCACAAAATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCCACCCA‑AGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGCTCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAACCACCGCTc                                                                                             >  1:99887/1‑248 (MQ=255)                                                                                                                                                                                                  gaggCATTAGCAGGCGAGACAACTACCGAAAAAGCGGCC‑CCCCCAGTTGCTGGGGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATTTTGTTTCCGACGCGCATTAAGCTAATTTATTGATCCCTTGGGATGAAACTCCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAACCACCGCTTAGACGATTGCGTTGGCTGCGGTGGCTGACAAAAATACGACAAAAGACg                                             <  1:313391/248‑1 (MQ=255)                                                                                                                                                                                                    aGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAACCACCGCTCAGACGATGGCGTTGGCTGCGGTGGCTGACAAAAATACGACAAAAGACGaa                                           <  1:78677/251‑1 (MQ=255)                                                                                                                                                                                                    aGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAACCACCGCTCAGACGATGGCGTTGGCTGCGGTGGCTGACAAAAATACGACAAAAGACGa                                            >  1:17325/1‑250 (MQ=255)                                                                                                                                                                                                              ggCGAGACAACTACCGACAAAGCGGCC‑CCCCCAGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAACCACCGCTCAGACGATGGCGTTGGCTGCGGTGGCTGACAAAAATACGACAAAAGCCGAAAAAGCGgatg                                 >  1:230381/1‑251 (MQ=255)                                                                                                                                                                                                                cGAGACAACTCGCGACAAAGCGGCC‑CCCCCAGTTGCTGTGGGCAACAGATACGCCCACTACAAAAGTCGAGCCTCTGATCAGCGATATTGTTTCCTACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCACCGACGTACAATCCACATCAACACCGTTAACCACCGCTCAGACGATGGCGTTGGCTGCGGTGGCTGACAAAAATACGACAAAAGACGAAAAAGCGgatgat                               <  2:212997/251‑1 (MQ=255)                                                                                                                                                                                                                                             c‑AGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAACCACCGCTCAGACGATGGCGTTGGCTGCGGTGGCTGACAAAAATACTACAAAAGACGAAAAAGCGGATGATCTGAATGAAGACGTCACCGCAAGCCTGAg  <  1:139480/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                           GGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGGATGAGTTCGCCCAGCGCGCCGCCATGAGGAAACCTGAATGATTCGCTTAGCGCCCTTGATTACCGCCGACGTTGACACCACCACATTGCCTGGCGGCAAAGCCAGCGATGCTGCACAAGATTTTCTCGCGTTGTTGAGCGAAGCATTAGCAGGCGAGACAACTACCGACAAAGCGGCC‑CCCC‑AGTTGCTGGTGGCAACAGATAAGCCCACGACAAAAGGCGAGCCGCTGATCAGCGATATTGTTTCCGACGCGCAACAAGCTAATTTACTGATCCCTGTGGATGAAACACCGCCTGTCATCAACGACGAACAATCCACATCAACACCGTTAACCACCGCTCAGACGATGGCGTTGGCTGCGGTGGCTGACAAAAATACGACAAAAGACGAAAAAGCGGATGATCTGAATGAAGACGTCACCGCAAGCCTGAG  >  NC_000913/2018295‑2018780

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