Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	913,779	A→G	72.4%	L153P (CTG→CCG)	hcp ←	hybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	913,779	0	A	G	72.4%	46.1 / 17.9	29	L153P (CTG→CCG)	hcp	hybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base G (2/19);  total (5/24)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

AGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTG  >  W3110S.gb/913722‑913839                                                          |                                                             agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCgag                                                 <  1:788137/71‑1 (MQ=255) agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCgag                                                 <  1:189332/71‑1 (MQ=255) agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCgag                                                 <  1:1597993/71‑1 (MQ=255) agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCgag                                                 <  1:1297050/71‑1 (MQ=255) agaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCgag                                                 <  1:1485155/71‑1 (MQ=255)  gaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCgag                                                 <  1:1343686/70‑1 (MQ=255)               ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCa                                    <  1:1129203/70‑1 (MQ=255)               ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCa                                    <  1:585987/70‑1 (MQ=255)               ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCa                                    <  1:1540338/70‑1 (MQ=255)               ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCa                                    <  1:1329742/70‑1 (MQ=255)               ttCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCa                                    <  1:1327910/70‑1 (MQ=255)                      taGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccgc                            <  1:1237130/71‑1 (MQ=255)                      taGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccgc                            <  1:443578/71‑1 (MQ=255)                      taGGCCGCCGCACCTTTCAGGCCATACAGGCACACCAGACGCAGGCCGAGAATGTTTTCGCCAATcgccgc                            <  1:876373/71‑1 (MQ=255)                              cgcACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCtt                    <  1:889314/71‑1 (MQ=255)                              cgcACCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCtt                    <  1:470669/71‑1 (MQ=255)                                 aCCTTTCAGGCCATACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:1951560/71‑1 (MQ=255)                                              tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:1347224/58‑1 (MQ=255)                                              tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:1845424/58‑1 (MQ=255)                                              tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:816953/58‑1 (MQ=255)                                              tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:863627/58‑1 (MQ=255)                                              tACAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:1520162/58‑1 (MQ=255)                                              tACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:90015/58‑1 (MQ=255)                                               aCAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGt                 <  1:179575/57‑1 (MQ=255)                                                cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg  >  1:128188/1‑71 (MQ=255)                                                cAGGCACAGCGGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg  >  1:1088730/1‑71 (MQ=255)                                                cAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg  >  1:1932158/1‑71 (MQ=255)                                                cAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg  >  1:1295174/1‑71 (MQ=255)                                                cAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTctg  >  1:1530800/1‑71 (MQ=255)                                                          |                                                             AGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGACGCAGGCCGAGAATGTTTTCGCCAATCGCCGCTTTATCTTTGTTAGGGGTAAATTCTG  >  W3110S.gb/913722‑913839

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