Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,387,862	C→T	E189K (GAA→AAA)	yfbN ←	uncharacterized protein YfbN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,387,862	0	C	T	88.9%	19.7 / ‑3.5	9	E189K (GAA→AAA)	yfbN	uncharacterized protein YfbN
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (4/4);  total (4/5)
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.95e-01

AACCTCTCTTGCAAGTTGGTCACAATGTATTATAAGCATTTTGAATATGTCATCCTCTGTAATCTCGTGAAAATTTTCGGCCTTGTCTTTTAGCTCATCATATCGGATTCTTGCAGTAACCAT  >  NC_000913/2387785‑2387907                                                                              |                                              aaCCTCTCTTGCAAGTTGGTCACAATGTATTATAAGCATTTTGAATATGTCATCCTCTGTAATCTCGTGAAAATTTTTgg                                             <  2:407066/80‑1 (MQ=255)  aCCTCTCTTGCAAGTTGGTCACAATGTATTATAAGCATTTTGAATATGTCATCCTCTGTAATCTCGTGAAAATTTTTGGcc                                           <  1:407067/81‑1 (MQ=255)         ttGCAAGTTGGTCACAATGTATTATAAGCATTTTGAATATGTCATCCTCTGTAATCTCGTGAAAATTTTCGGCCTTGTc                                      <  2:407068/79‑1 (MQ=255)          tGCAAGTTGGTCACAATGTATTATAAGCATTTTGAATATGTCATCCTCTGTAATCTCGTGAAAATTTTTGGCCTTGTCttt                                   >  1:407084/1‑81 (MQ=255)                      acaATGTATTATAAGCATTTTGAATATGTCATCCTCTGTAATCTCGTGAAAATTTTTGGCCTTGTCTTTTAGCTCATCata                       >  2:407080/1‑81 (MQ=255)                          tGTATTATAAGCATTTTGAATATGTCATCCTCTGTAATCTCGTGAAAATTTTTGGCCTTGTCTTTTAGCTCATCATATCg                    >  2:407081/1‑80 (MQ=255)                                 taAGCATTTTGAATATGTCATCCTCTGTAATCTCGTGAAAATTTTTGGCCTTGTCTTTTAGCTCATCATATCGGATTCTTg            >  1:407090/1‑81 (MQ=255)                                      aTTTTGAATATGTCATCCTCTGTAATCTCGTGAAAATTTTTGGCCTTGTCTTTTAGCTCATCATATCGGATTCTTGCAGTa       <  1:407069/81‑1 (MQ=255)                                            aaTATGTCATCCTCTGTAATCTCGTGAAAATTTTTGGCCTTGTCTTTTAGCTCATCATATCGGATTCTTGCAGTAACCAt  <  1:407070/80‑1 (MQ=255)                                                                              |                                              AACCTCTCTTGCAAGTTGGTCACAATGTATTATAAGCATTTTGAATATGTCATCCTCTGTAATCTCGTGAAAATTTTCGGCCTTGTCTTTTAGCTCATCATATCGGATTCTTGCAGTAACCAT  >  NC_000913/2387785‑2387907

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