Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,854,548	+A	coding (1146/1359 nt)	ydjE ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,854,545	1	.	A	98.4%	170.8 / ‑3.1	61	coding (1149/1359 nt)	ydjE	predicted transporter
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (0/60);  total (0/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GAACCGCGCAGGCGT‑AAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGC  >  W3110S.gb/1854531‑1854591                |                                                gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:589597/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:272893/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:289626/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:294492/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:33553/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:379646/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:495634/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:501090/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:505287/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:52492/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:533956/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:534308/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:538178/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:558065/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:588964/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:186164/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:633586/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:99195/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:659777/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:687222/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:713820/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:767783/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:800721/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:848934/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:864466/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:952162/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:953788/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:956866/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:963188/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:970877/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1262001/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1014487/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1048959/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1049167/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1059933/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1083023/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1091711/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1104475/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1110634/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1122832/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1140318/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1143071/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1163092/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1169137/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1257686/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1263657/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1330882/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1378087/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1379051/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1414390/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1477271/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1477325/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1488956/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1522618/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1532991/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1540453/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1556529/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:1607433/62‑1 (MQ=255) gAACCGCGCAGGCGTAAAA‑TGCGTTGGCCAAAGCTCAGGGATATAAACCGCCGACGCGAAGc  <  1:1036169/62‑1 (MQ=255) gAACCGCGCAGGCGT‑AAATTGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:644335/62‑1 (MQ=255)          aGGCGTAAAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGc  <  1:993655/53‑1 (MQ=255)                |                                                GAACCGCGCAGGCGT‑AAA‑TGCGTTGGCCAAAGCTCCGGGATATAAACCGCCGACGCGAAGC  >  W3110S.gb/1854531‑1854591

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: ‑