Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,266,995	A→T	25.0%	E864V (GAA→GTA)	fimV →	motility protein FimV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,266,995	0	A	T	25.0%	22.5 / 6.2	16	E864V (GAA→GTA)	fimV	motility protein FimV
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base T (2/2);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CCGAGCTGGACAAGCTGTCGCAGAGCCTCGAAACGCCATCGCTGGAGCCGCACTTCACCAGCGAAGACGCGGTGGCGCAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTC  >  NC_002947/2266860‑2267105                                                                                                                                        |                                                                                                               ccGAGCTGGACAAGCTGTCGCAGAGCCTCGAAACGCCATCGCTGGAGCCGCACTTCACCAGCGAAGACGCGGTGGCGCAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGa                                                                                                                >  2:149268/1‑136 (MQ=255)       tGGACAAGCTGTCGCAGAGCCTCGAAACGCCATCGCTGGAGCCGCACTTCACCAGCGAAGACGCGGTGGCGCAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGc                                                                                                          <  1:61283/136‑1 (MQ=255)            aaGCTGTCGCAGAGCCTCGAAACGCCATCGCTGGAGCGGCAGTTCACCAGCGTAGAGGCGGTGGCGCAACGAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCTGTGACGAAGTCGccacct                                                                                                     <  2:369682/136‑2 (MQ=255)                      agagCCTCGAAACGCCATCGCTGGAGCCGCACTTCACCAGCGAAGACGCGGTGGCGCAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGTAGTCGCCACCAAGCTCGACCt                                                                                           >  1:42248/1‑136 (MQ=255)                      agagCCTCGAAACGCCATCGCTGGAGCCGCACTTCACCAGCGAAGACGCGGTGGCGCAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGTAGTCGCCACCAAGCTCGACCt                                                                                           >  1:44214/1‑136 (MQ=255)                             cGAAACGCCATCGCTGGAGCCGCACTTCACCAGCGAAGACGCGGTGGCGCAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGt                                                                                    >  1:239755/1‑136 (MQ=255)                                       tCGCTGGAGCCGCACTTCACCAGCGAAGACGCGGTGGCGCAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGTAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCg                                                                          <  2:44214/136‑1 (MQ=255)                                       tCGCTGGAGCCGCACTTCACCAGCGAAGACGCGGTGGCGCAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGTAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCg                                                                          <  2:42248/136‑1 (MQ=255)                                              aGCCGCACTTCACCAGCGAAGACGCGGTGGCGCAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATggg                                                                   <  2:200076/136‑1 (MQ=255)                                                      ttCACCAGCGAAGACGCGGTGGCGCAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGcc                                                           >  1:218980/1‑136 (MQ=255)                                                                           gcgcAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTcg                                      >  2:23608/1‑136 (MQ=255)                                                                           gcgcAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTcg                                      >  2:28566/1‑136 (MQ=255)                                                                                                    acgacCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGcgacga             >  2:292565/1‑136 (MQ=255)                                                                                                          tGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGAcagcca       >  1:120794/1‑136 (MQ=255)                                                                                                           ggACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGAcagccag      >  1:208553/1‑136 (MQ=255)                                                                                                               ttcgacttcTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTc  <  2:208553/136‑1 (MQ=255)                                                                                                               ttcgacttcTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTc  <  2:239755/136‑1 (MQ=255)                                                                                                                                        |                                                                                                               CCGAGCTGGACAAGCTGTCGCAGAGCCTCGAAACGCCATCGCTGGAGCCGCACTTCACCAGCGAAGACGCGGTGGCGCAACCAGAGCCTGAGCCTTTGGACGACCTGGACTTCGACTTCTTCTCCGGCAGTGACGAAGTCGCCACCAAGCTCGACCTGGCCCGTGCCTATATCGACATGGGTGACAGCCAGGGCGCGCGCGACATCCTCGACGAAGTGGTCAAGGATGGCGACGACAGCCAGCGTC  >  NC_002947/2266860‑2267105

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