Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,227,635	+T	coding (216/3684 nt)	metH →	homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,227,632	1	.	T	77.8%	111.3 / 33.4	45	coding (213/3684 nt)	metH	homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/10);  new base T (35/0);  total (35/10)
Fisher's exact test for biased strand distribution p-value = 3.13e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TAAACCGGAAGTGATCGCCGCTATCCACAACGCCTAC‑TTTGAAGCGGGCGCGGATATCATCGAAACCAACA  >  W3110S.gb/4227596‑4227666                                      |                                   tAAACCGGAAGTGATCGCCGCTATCCACAACGCCTAC‑TTTGAAGCGGGCGCGGATatcat             <  1:1670933/60‑1 (MQ=255) tAAACCGGAAGTGATCGCCGCTATCCACAACGCCTAC‑TTTGAAGCGGGCGCGGATatcat             <  1:2356438/60‑1 (MQ=255) tAAACCGGAAGTGATCGCCGCTATCCACAACGCCTAC‑TTTGAAGCGGGCGCGGATatcat             <  1:251526/60‑1 (MQ=255) tAAACCGGAAGTGATCGCCGCTATCCACAACGCCTAC‑TTTGAAGCGGGCGCGGATatcat             <  1:2215842/60‑1 (MQ=255) tAAACCGGAAGTGATCGCCGCTATCCACAACGCCTAC‑TTTGAAGCGGGCGCGGATatcat             <  1:2771913/60‑1 (MQ=255) tAAACCGGAAGTGATCGCCGCTATCCACAACGCCTAC‑TTTGAAGCGGGCGCGGATatcat             <  1:485353/60‑1 (MQ=255) tAAACCGGAAGTGATCGCCGCTATCCACAACGCCTAC‑TTTGAAGCGGGCGCGGATatcat             <  1:1377200/60‑1 (MQ=255) tAAACCGGAAGTGATCGCCGCTATCCACAACGCCTAC‑TTTGAAGCGGGCGCGGATatcat             <  1:947907/60‑1 (MQ=255) tAAACCGGAAGTGATCGCCGCTATCCACAACGCCTAC‑TTTGAAGCGGGCGCGGATatcat             <  1:87346/60‑1 (MQ=255) tAAACCGGAAGTGATCGCCGCTATCCACAACGCCTAC‑TTTGAAGCGGGCGCGGATatcat             <  1:1369992/60‑1 (MQ=255)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGa                   >  1:1164999/1‑35 (MQ=37)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATatcat             >  1:2783981/1‑41 (MQ=38)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGa          >  1:2035230/1‑44 (MQ=38)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:2241357/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:866985/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:428385/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:2878139/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:540372/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:618327/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:2752254/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:2718395/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:807692/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:2685594/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:894953/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:836581/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:1707880/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:1036043/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:1131242/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:1155593/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:1293456/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:1332753/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:1449000/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:146503/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:1589727/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:1696161/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:2382315/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:1901978/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:1942074/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:1953604/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:2150210/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAaca  >  1:2270501/1‑52 (MQ=39)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAac   >  1:371648/1‑51 (MQ=38)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAac   >  1:2735232/1‑51 (MQ=38)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAac   >  1:2209593/1‑51 (MQ=38)                     cTATCCACAACGCCTACTTTTGAAGCGGGCGCGGATATCATCGAAACCAac   >  1:1000724/1‑51 (MQ=38)                                      |                                   TAAACCGGAAGTGATCGCCGCTATCCACAACGCCTAC‑TTTGAAGCGGGCGCGGATATCATCGAAACCAACA  >  W3110S.gb/4227596‑4227666

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