Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,182,652	C→T	37.5%	V464M (GTG→ATG)	gspD ←	general secretory pathway component, cryptic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,182,652	0	C	T	37.5%	25.4 / 31.4	32	V464M (GTG→ATG)	gspD	general secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (10/10);  new base T (8/4);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTGGCGAG  >  W3110S.gb/4182603‑4182715                                                  |                                                                  ttAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGa                                            >  1:430955/1‑71 (MQ=255) ttAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGa                                            >  1:42708/1‑71 (MQ=255)   aaaGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCATATCCTGCCC‑CCACGTTGAAGGAc                                           >  1:317210/1‑71 (MQ=255)   aaaGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCATATCCTGCCC‑CCACGTTGAAGGAc                                           >  1:30204/1‑71 (MQ=255)   aaaGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCATATCCTG‑cc‑ccc                                                      >  1:926730/1‑58 (MQ=255)   aaaGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCATATCCTG‑CCACCACGTTGAAGGAc                                           >  1:1215581/1‑71 (MQ=255)   aaaGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCATATCCTG‑CC‑CCACGTTGAAGGAcg                                          >  1:70912/1‑71 (MQ=255)   aaaGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCATATCCTG‑CC‑CCACGTTGAAGGAcg                                          >  1:1471479/1‑71 (MQ=255)   aaaGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCATATCCTG‑CC‑CCACGTTGAAGGAcg                                          >  1:1568997/1‑71 (MQ=255)          ttATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAg                                              <  1:1136531/60‑1 (MQ=255)          ttATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAg                                              <  1:605136/60‑1 (MQ=255)          ttATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAg                                              <  1:348120/60‑1 (MQ=255)          ttATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAg                                              <  1:114196/60‑1 (MQ=255)          ttATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAg                                              <  1:127520/60‑1 (MQ=255)          ttATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAg                                              <  1:132210/60‑1 (MQ=255)          ttATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGACGCGAGttt                                   <  1:1098122/71‑1 (MQ=255)          ttATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGACGCGAGttt                                   <  1:1492063/71‑1 (MQ=255)           tATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGACGCGAGttt                                   <  1:1916453/70‑1 (MQ=255)                            tgtgACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGAt                 >  1:955458/1‑71 (MQ=255)                            tgtgACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGAt                 >  1:948101/1‑71 (MQ=255)                            tgtgACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGAt                 >  1:1355656/1‑71 (MQ=255)                                   cGGATAGCACCGGCATATCCTG‑CC‑CCACGTTGAAGGAcgc                                         >  1:486531/1‑40 (MQ=255)                                       tAGCACCGGCACATCCTG‑CC‑CCCCGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTgg      >  1:248931/1‑71 (MQ=255)                                       tAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTgg      >  1:452000/1‑71 (MQ=255)                                       tAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTgg      >  1:69179/1‑71 (MQ=255)                                       tAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTgg      >  1:1573302/1‑71 (MQ=255)                                       tAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTgg      >  1:272207/1‑71 (MQ=255)                                           aCCGGCATATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTGGCGAg  <  1:1214740/71‑1 (MQ=255)                                           aCCGGCATATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTGGCGAg  <  1:747453/71‑1 (MQ=255)                                           aCCGGCATATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTGGCGAg  <  1:786350/71‑1 (MQ=255)                                           aCCGGCATATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTGGCGAg  <  1:168665/71‑1 (MQ=255)                                           aCCGGCACATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTGGCGAg  <  1:591984/71‑1 (MQ=255)                                                  |                                                                  TTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTG‑CC‑CCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTGGCGAG  >  W3110S.gb/4182603‑4182715

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