Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,021,506	C→T	T512I (ACT→ATT)	xdhD →	fused predicted xanthine/hypoxanthine oxidase, molybdopterin‑binding subunit and Fe‑S binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,021,506	0	C	T	97.1%	106.4 / ‑3.8	36	T512I (ACT→ATT)	xdhD	fused predicted xanthine/hypoxanthine oxidase, molybdopterin‑binding subunit and Fe‑S binding subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (17/17);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGTCGC  >  W3110S.gb/3021452‑3021562                                                       |                                                         gAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAAc                                          <  1:3531863/71‑1 (MQ=255) gAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAAc                                          <  1:1051290/71‑1 (MQ=255) gAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAAc                                          <  1:3126252/71‑1 (MQ=255) gAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAAc                                          <  1:2937935/71‑1 (MQ=255) gAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAAc                                          <  1:497581/71‑1 (MQ=255)                               tGAAGATGGATTTCCGCGCCATCATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACgg                <  1:3934650/67‑1 (MQ=255)                               tGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACgg                <  1:2223965/67‑1 (MQ=255)                               tGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACgg                <  1:3441650/67‑1 (MQ=255)                               tGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACgg                <  1:3584532/67‑1 (MQ=255)                               tGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACgg                <  1:3865427/67‑1 (MQ=255)                               tGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACgg                <  1:2539934/67‑1 (MQ=255)                               tGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACgg                <  1:1813491/67‑1 (MQ=255)                               tGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACgg                <  1:759203/67‑1 (MQ=255)                               tGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACgg                <  1:1573693/67‑1 (MQ=255)                               tGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACgg                <  1:1558154/67‑1 (MQ=255)                               tGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACgg                <  1:1384214/67‑1 (MQ=255)                               tGAAGATGGATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACgg                <  1:1226701/67‑1 (MQ=255)                                gAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCgg           >  1:1618014/1‑71 (MQ=255)                                      ggATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTg      >  1:3727505/1‑70 (MQ=255)                                      ggATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTg      >  1:3384003/1‑70 (MQ=255)                                      ggATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGt     >  1:3898209/1‑71 (MQ=255)                                      ggATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGt     >  1:3555070/1‑71 (MQ=255)                                      ggATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGt     >  1:3485978/1‑71 (MQ=255)                                      ggATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGt     >  1:3490662/1‑71 (MQ=255)                                       gATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGTACCGGCGCTGTc    >  1:801373/1‑71 (MQ=255)                                       gATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGTc    >  1:325972/1‑71 (MQ=255)                                       gATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGTc    >  1:3087369/1‑71 (MQ=255)                                       gATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGTc    >  1:2033823/1‑71 (MQ=255)                                       gATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGTc    >  1:68537/1‑71 (MQ=255)                                       gATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGTc    >  1:1726844/1‑71 (MQ=255)                                       gATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGTc    >  1:1521770/1‑71 (MQ=255)                                       gATTTCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGATGTc    >  1:2059728/1‑71 (MQ=255)                                       gATTTCCGCGCCAACAATGGCCCTTACGGCAACCACTCACACACCGTACCGTGTAACg                 >  1:3209321/1‑58 (MQ=255)                                         tttCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGTCg   >  1:2572631/1‑70 (MQ=255)                                         tttCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGTCGc  >  1:2108961/1‑71 (MQ=255)                                         tttCCGCGCCAACATTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGTCGc  >  1:656381/1‑71 (MQ=255)                                                       |                                                         GAGCGAAAAAAGATGGTCGCCTGACGGCAGTGAAGATGGATTTCCGCGCCAACACTGGCCCTTACGGCAACCACTCACTCACCGTACCGTGTAACGGACCGGCGCTGTCGC  >  W3110S.gb/3021452‑3021562

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