Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	182,263	C→G	L460L (CTC→CTG)	degP →	serine endoprotease (protease Do), membrane‑associated

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	182,263	0	C	G	98.2%	190.5 / ‑4.1	57	L460L (CTC→CTG)	degP	serine endoprotease (protease Do), membrane‑associated
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (56/0);  total (56/1)
Fisher's exact test for biased strand distribution p-value = 1.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CGACAGCAAACCGTCTGTGCTGGCACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGC  >  W3110S.gb/182236‑182303                            |                                         cGACAGCAAACCGTCTGTGGTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCt          >  1:454890/1‑60 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAATATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1470701/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACGTTCAGCGCGGCGACAGCACCATCTACCTg         >  1:2645893/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAgcg                               >  1:882311/1‑39 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCTCCAt                >  1:2521590/1‑54 (MQ=39) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCAcc                  >  1:416983/1‑52 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTa             >  1:379461/1‑57 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCt          >  1:2489397/1‑60 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCt          >  1:1024043/1‑60 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCt          >  1:2291207/1‑60 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCt          >  1:2622763/1‑60 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCt          >  1:1805299/1‑60 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCt          >  1:656607/1‑60 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCt          >  1:668649/1‑60 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCt          >  1:1312052/1‑60 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCt          >  1:807937/1‑60 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:2827447/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:2417218/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:856128/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:828417/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:260008/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:816973/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:2783126/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:2814702/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:614348/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:333968/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:732772/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:405197/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:451393/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:635331/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:550226/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1743904/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1190574/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1231624/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1317388/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1361929/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1411005/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1440555/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1466088/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1526737/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1613283/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1710619/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:2506403/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1791337/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1797086/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1848732/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1916118/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:2020788/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:2106619/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:2132735/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:2210077/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:2233935/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:2344755/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAACATTCAGAGCGGCGACAGCACCATCTACCTg         >  1:1540984/1‑61 (MQ=255) cGACAGCAAACCGTCTGTGCTGGCACTGAAAATTCAGCGCGGCGACAg                      >  1:247886/1‑48 (MQ=38) cGACAGCAAACCGTATGTGCTGGCACTGAACATTCAGCGCGGCGACAGCACCATCTACCTg         >  1:1301161/1‑61 (MQ=255)       cAAACCGTCTGTGCTGGCACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGc  <  1:366061/62‑1 (MQ=255)                            |                                         CGACAGCAAACCGTCTGTGCTGGCACTCAACATTCAGCGCGGCGACAGCACCATCTACCTGTTAATGC  >  W3110S.gb/182236‑182303

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