Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,898,954	+T	coding (120/474 nt)	yiaI →	predicted hydrogenase, 4Fe‑4S ferredoxin‑type component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,898,951	1	.	T	82.8%	164.4 / 30.9	58	coding (117/474 nt)	yiaI	predicted hydrogenase, 4Fe‑4S ferredoxin‑type component
Reads supporting (aligned to +/- strand):  ref base . (10/0);  new base T (48/0);  total (58/0)
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.36e-01

GCGCAGTGTCGCATCATGAGAATCAGGATTGCGCTGCGTTGTCACCAGACGAG‑TTTATTTCCCGTATTCGTGTCATTAAAGAC  >  W3110S.gb/3898899‑3898981                                                      |                               gcgcAGTGTCGCATCATGAGAATCAGGATTGCGCTGCGTTGTCACCAGACGAG‑TTTATTTCCCGt                    >  1:2315539/1‑65 (MQ=255) gcgcAGTGTCGCATCATGAGAATCAGGATTGCGCTGCGTTGTCACCAGACGAG‑TTTATTTCCCGt                    >  1:1420446/1‑65 (MQ=255) gcgcAGTGTCGCATCATGAGAATCAGGATTGCGCTGCGTTGTCACCAGACGAG‑TTTATTTCCCGt                    >  1:1417554/1‑65 (MQ=255) gcgcAGTGTCGCATCATGAGAATCAGGATTGCGCTGCGTTGTCACCAGACGAG‑TTTATTTCCCGt                    >  1:1012328/1‑65 (MQ=255) gcgcAGTGTCGCATCATGAGAATCAGGATTGCGCTGCGTTGTCACCAGACGAG‑TTTATTTCCCGt                    >  1:1319945/1‑65 (MQ=255) gcgcAGTGTCGCATCATGAGAATCAGGATTGCGCTGCGTTGTCACCAGACGAG‑TTTATTTCCCGt                    >  1:1793240/1‑65 (MQ=255) gcgcAGTGTCGCATCATGAGAATCAGGATTGCGCTGCGTTGTCACCAGACGAG‑TTTATTTCCAGt                    >  1:2229126/1‑65 (MQ=255)            catcatGAGAATCAGGATTGCGCTGCGTTGTCACCAGACGAG‑TTTATTTCCCGTATTCGTGTCa          >  1:554174/1‑64 (MQ=255)            catcatGAGAATCAGGATTGCGCTGCGTTGTCACCAGACGAG‑TTTATTTCCCGTATTCGTGTCa          >  1:723207/1‑64 (MQ=255)            catcatGAGAATCAGGATTGCGCTGCGTTGTCACCAGACGAG‑TTTATTTCCCGTATTCGAGTCa          >  1:1978878/1‑64 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTccc                      >  1:291030/1‑45 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCGGTATTCGTGTCATTAAAGAc  >  1:1282253/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTAt                  >  1:1320302/1‑49 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTaaa     >  1:2233254/1‑62 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:2204572/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:2217139/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:2019080/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:368871/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:2394733/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:2420346/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:24537/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:2520255/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:317024/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:332816/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1878794/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:394976/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:454987/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:482176/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:678848/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:732530/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:827753/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:873801/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:979218/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1505173/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1072928/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1082060/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1306182/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:130946/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1323348/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1327990/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1368224/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1460536/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1476683/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1503299/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1974747/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1558813/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1563982/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1568074/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1591895/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1598923/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1608754/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1713632/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1717040/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1762274/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1771701/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1885789/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1932756/1‑65 (MQ=255)                    gaATCAGGATTGCGCTGCGTTGTCACCAGACGAG‑tt                              >  1:554111/1‑36 (MQ=255)                    gaATCAGGAATGCGCTGCGTTGTCACCAGACGAGTTTTATTTCCCGTATTCGTGTCATTAAAGAc  >  1:1483310/1‑65 (MQ=255)                                                      |                               GCGCAGTGTCGCATCATGAGAATCAGGATTGCGCTGCGTTGTCACCAGACGAG‑TTTATTTCCCGTATTCGTGTCATTAAAGAC  >  W3110S.gb/3898899‑3898981

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