Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,573,836	A→G	48.5%	V200A (GTG→GCG)	yihN ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,573,836	0	A	G	48.5%	‑3.2 / 45.0	33	V200A (GTG→GCG)	yihN	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (8/9);  new base G (7/9);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

AGCCAGATAGTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTACGGC  >  W3110S.gb/3573767‑3573897                                                                      |                                                              aGCCAGATAGTTTTGTCTTTCAGGACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCa                                                               <  1:1810701/70‑1 (MQ=255) aGCCAGATAGTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCa                                                               <  1:1971224/70‑1 (MQ=255) aGCCAGATAGTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCa                                                               <  1:2832015/70‑1 (MQ=255) aGCCAGATAGTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCa                                                               <  1:2332108/70‑1 (MQ=255)      gATAGTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCGCCTcga                                                         <  1:3334216/71‑3 (MQ=255)      gATAGTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCGCCTCaa                                                         <  1:1309900/71‑1 (MQ=255)      gATAGTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCGCCTCaa                                                         <  1:1068949/71‑1 (MQ=255)      gATAGTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCGCCTCaa                                                         <  1:716441/71‑1 (MQ=255)      gATAGTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCGCCTCaa                                                         <  1:3035383/71‑1 (MQ=255)      gATAGTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCGCCTCaa                                                         <  1:2832232/71‑1 (MQ=255)      gATAGTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCGCCTCaa                                                         <  1:2382480/71‑1 (MQ=255)             ttGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCACCTCAacgga                                                     <  1:1821918/68‑1 (MQ=255)                                     cTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGa                           >  1:3302721/1‑70 (MQ=255)                                     cTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGa                           >  1:2800660/1‑70 (MQ=255)                                     cTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGa                           >  1:2163729/1‑70 (MQ=255)                                                 gATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat              <  1:93212/71‑1 (MQ=255)                                                 gATGCTCCGTCTTCTTTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat              <  1:2007111/71‑1 (MQ=255)                                                 gATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat              <  1:3239731/71‑1 (MQ=255)                                                 gATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat              <  1:1613989/71‑1 (MQ=255)                                                 gATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat              <  1:1652781/71‑1 (MQ=255)                                                 gATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat              <  1:68731/71‑1 (MQ=255)                                                 gATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat              <  1:1759041/71‑1 (MQ=255)                                                 gATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat              <  1:2879412/71‑1 (MQ=255)                                                 gATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat              <  1:2736884/71‑1 (MQ=255)                                                 gATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAaataat              <  1:210340/71‑1 (MQ=255)                                                   tGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTc                                   >  1:3178777/1‑48 (MQ=255)                                                   tGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTc                                   >  1:1714544/1‑48 (MQ=255)                                                   tGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTc                                   >  1:3335682/1‑48 (MQ=255)                                                   tGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTc                                   >  1:2575488/1‑48 (MQ=255)                                                   tGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTc                                   >  1:907955/1‑48 (MQ=255)                                                              cttTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTACGGc  >  1:1871932/1‑70 (MQ=255)                                                              cttTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTACGGc  >  1:2000235/1‑70 (MQ=255)                                                              cttTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTACGGc  >  1:2208950/1‑70 (MQ=255)                                                              cttTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTACGGc  >  1:3304657/1‑70 (MQ=255)                                                              cttTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTACGGc  >  1:1671788/1‑70 (MQ=255)                                                              cttTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTACGGc  >  1:535701/1‑70 (MQ=255)                                                              cttTTTTCGCCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTACGGc  >  1:84773/1‑70 (MQ=255)                                                                      |                                                              AGCCAGATAGTTTTGTCTTTCAGCACCGAGGTCATACTGGTGTTTTTCGATGCTCCGTCTTCTTTTTTCACCTCAACGGACGGTGCCTCTTCTTTGTCATTCAGGACAAAGAAAATAATAATGCCTACGGC  >  W3110S.gb/3573767‑3573897

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