Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,326,279	G→T	T33N (ACC→AAC)	phnG ←	carbon‑phosphorus lyase complex subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,326,279	0	G	T	79.2%	80.8 / 26.9	48	T33N (ACC→AAC)	phnG	carbon‑phosphorus lyase complex subunit
Reads supporting (aligned to +/- strand):  ref base G (7/3);  new base T (38/0);  total (45/3)
Fisher's exact test for biased strand distribution p-value = 6.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAG  >  W3110S.gb/4326236‑4326334                                            |                                                        caCCAGGCCAGTTTCATCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2729445/2‑59 (MQ=255) caCCAGGCCAGTTTCATCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:1312875/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGtt                                                  >  1:118655/2‑51 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGt                                                   >  1:1398750/2‑50 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:3058396/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:238568/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:242253/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2491288/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2494442/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:1054139/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2557674/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2677047/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:290454/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2909585/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2930417/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2950406/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2998977/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2546985/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:469216/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:473128/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:546563/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:547451/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:636555/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:685951/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:755054/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:853262/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:960348/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:991530/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:235667/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:1867664/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2115271/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2157655/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2270107/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:1156590/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2274207/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:2316772/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTCAGCGCGt                                          >  1:231950/2‑59 (MQ=255) caCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAATCGGCGTTGATGTTAAGCGCGt                                          >  1:916582/2‑59 (MQ=255)   ccAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGtt                                         <  1:80745/58‑1 (MQ=255)   ccAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGtt                                         <  1:740838/58‑1 (MQ=255)   ccAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGtt                                         <  1:1780953/58‑1 (MQ=255)                    gcgcGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGc                    >  1:567627/1‑62 (MQ=255)                                      tCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAg  >  1:1385721/1‑62 (MQ=255)                                      tCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAg  >  1:3116764/1‑62 (MQ=255)                                      tCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAg  >  1:1780638/1‑62 (MQ=255)                                      tCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAg  >  1:1096283/1‑62 (MQ=255)                                      tCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAg  >  1:297124/1‑62 (MQ=255)                                      tCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAg  >  1:2961999/1‑62 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGcag          >  1:2802904/1‑50 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:574795/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:2350156/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:98918/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:1109666/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:2410404/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:2533484/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:2684174/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:1166852/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:2709767/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:1179031/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:659850/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:1182897/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:2102103/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:1216625/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:1358766/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:325581/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:2823455/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:2093660/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:202102/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCtg          >  1:1866423/1‑52 (MQ=255)                                        ggcggTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCt           >  1:1218862/1‑51 (MQ=255)                                            |                                                        TACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAG  >  W3110S.gb/4326236‑4326334

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