Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,594,104	G→A	H138Y (CAT→TAT)	ypfI ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,594,104	0	G	A	90.0%	109.9 / 4.3	40	H138Y (CAT→TAT)	ypfI	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base A (20/16);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GAATGGCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATGCTGGACAAAATGCGGCGTCGCAATAGGGTCAGGGCAATCACTCCAGCGCAGCGAGTC  >  W3110S.gb/2594042‑2594161                                                               |                                                          gaATGGCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACaa                                                   <  1:211601/71‑1 (MQ=255)      gCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATGCTGGACAAAATgc                                              >  1:106150/1‑71 (MQ=255)      gCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATGCTGGACAAAATgc                                              >  1:3828194/1‑71 (MQ=255)      gCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATGCTGGACAAAATgc                                              >  1:143760/1‑71 (MQ=255)      gCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATgc                                              >  1:114515/1‑71 (MQ=255)      gCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATgc                                              >  1:3974955/1‑71 (MQ=255)      gCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATgc                                              >  1:4006900/1‑71 (MQ=255)      gCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATgc                                              >  1:342192/1‑71 (MQ=255)      gCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATgc                                              >  1:338245/1‑71 (MQ=255)      gCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATgc                                              >  1:698353/1‑71 (MQ=255)      gCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATgc                                              >  1:3370617/1‑71 (MQ=255)      gCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATgc                                              >  1:456692/1‑71 (MQ=255)      gCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATgc                                              >  1:584396/1‑71 (MQ=255)              tGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACaa                                                   <  1:1197915/58‑1 (MQ=255)               gccgccAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGc                                      <  1:1192430/70‑1 (MQ=255)               gccgccAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGc                                      <  1:2719367/70‑1 (MQ=255)               gccgccAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGc                                      <  1:247724/70‑1 (MQ=255)               gccgccAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGc                                      <  1:1099626/70‑1 (MQ=255)               gccgccAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGc                                      <  1:1885280/70‑1 (MQ=255)               gccgccAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGc                                      <  1:144289/70‑1 (MQ=255)                         gATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATGCTGGACAAAATGCGGCGTCGCAATAGGGTCAg                           >  1:2133246/1‑71 (MQ=255)                         gATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAg                           >  1:3086915/1‑71 (MQ=255)                         gATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAg                           >  1:562499/1‑71 (MQ=255)                         gATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAg                           >  1:1755712/1‑71 (MQ=255)                         gATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAg                           >  1:3412577/1‑71 (MQ=255)                         gATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAg                           >  1:1107325/1‑71 (MQ=255)                         gATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAg                           >  1:3456658/1‑71 (MQ=255)                         gATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAg                           >  1:3596388/1‑71 (MQ=255)                         gATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAg                           >  1:1005874/1‑71 (MQ=255)                         gATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAg                           >  1:1308868/1‑71 (MQ=255)                         gATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCCCAATAGGGTCAg                           >  1:3728884/1‑71 (MQ=255)                             gCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGt                                         <  1:894788/53‑1 (MQ=255)                             gCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGt                                         <  1:3805672/53‑1 (MQ=255)                             gCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGt                                         <  1:3206326/53‑1 (MQ=255)                               cTCGTTATCCGCCGTAAGTACGCGTTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAGGGCAAt                     >  1:1892443/1‑71 (MQ=255)                                               aGTACGCGGTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAGGGCAATCACTCCAGCGCAGCGa     <  1:392658/71‑1 (MQ=255)                                                   cgcgTTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAGGGCAATCACTCCAGCGCAGCGAGTc  <  1:3655490/70‑1 (MQ=255)                                                   cgcgTTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAGGGCAATCACTCCAGCGCAGCGAGTc  <  1:3299927/70‑1 (MQ=255)                                                   cgcgTTTGAGATACTGGACAAAATGCGGCGTCGCAATAGGGTCAGGGCAATCACTCCAGCGCAGCGAGTc  <  1:1716325/70‑1 (MQ=255)                                                     cgTTTGAGATACTGGACAAAATGCGGCGTTGCAATAGGGTCAGGGCAATCACTCCAGCGCAGCGAGTc  <  1:1344146/68‑1 (MQ=255)                                                               |                                                          GAATGGCTGGTTTTGCCGCCAGAGGATAGCCTCGTTATCCGCCGTAAGTACGCGTTTGAGATGCTGGACAAAATGCGGCGTCGCAATAGGGTCAGGGCAATCACTCCAGCGCAGCGAGTC  >  W3110S.gb/2594042‑2594161

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