Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,894,068	T→C	27.9%	M86V (ATG→GTG)	yiaO ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,894,068	0	T	C	27.9%	62.6 / 28.6	43	M86V (ATG→GTG)	yiaO	predicted transporter
Reads supporting (aligned to +/- strand):  ref base T (18/13);  new base C (2/10);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGGACATCAAGCAAGTTCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTGCGCGTTACCGAG  >  W3110S.gb/3894002‑3894119                                                                   |                                                    aaGGGACATCAAGCAAGTTCATCAC‑‑GTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCACATCGAt                                               <  1:521996/71‑1 (MQ=255)   gtgACATCAAGCAAGTTCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCACATCGAt                                               <  1:2856407/69‑1 (MQ=255)   gggACATCAAGCAAGTTCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCACATCGAt                                               <  1:119939/71‑1 (MQ=255)   gggACATCAAGCAAGTTCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCACATCGAt                                               <  1:3083547/71‑1 (MQ=255)   gggACATCAAGCAAGTTCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCACATCGAt                                               <  1:3122463/71‑1 (MQ=255)   gggACATCAAGCAAGTTCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCACATCGAt                                               <  1:1909600/71‑1 (MQ=255)   gggACATCAAGCAAGTTCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCACATCGAt                                               <  1:461619/71‑1 (MQ=255)   gggACATCAAGCAAGTTCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCACATCGAt                                               <  1:1988586/71‑1 (MQ=255)   gggACATCAAGCAAGTTCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCACATCGAt                                               <  1:2127472/71‑1 (MQ=255)           aagcaagTTCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCATATCGAt                                               >  1:499006/1‑63 (MQ=255)                  ttCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTa                                 >  1:2643461/1‑70 (MQ=255)                  ttCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTAc                                >  1:2119181/1‑71 (MQ=255)                  ttCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTAc                                >  1:449667/1‑71 (MQ=255)                  ttCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTAc                                >  1:718375/1‑71 (MQ=255)                                tAACCCGGCAAAGTTATTCGAGCCGGACATTTCCACATCGATGGTGCCGCCACGTACGCCGCTGATCATCg                  >  1:394111/1‑71 (MQ=255)                                tAACCCGGCAAAGTTATTCGAGCCGGACATTTCCACATCGATGGTGCCGCCACGTACGCCGCTGATCATCg                  >  1:367060/1‑71 (MQ=255)                                 aaCCCGGCAAAGTTATTCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCg                  >  1:2098479/1‑70 (MQ=255)                                 aaCCCGGCAAAGTTATTCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGc                 >  1:827563/1‑71 (MQ=255)                                 aaCCCGGCAAAGTTATTCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGc                 >  1:1110364/1‑71 (MQ=255)                                 aaCCCGGCAAAGTTATTCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGc                 >  1:907703/1‑71 (MQ=255)                                 aaCCCGGCAAAGTTATTCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGc                 >  1:2493850/1‑71 (MQ=255)                                     cGGCAAAGTTATTCGAGCCGGACATTTCCACATCGAt                                               <  1:828175/37‑1 (MQ=255)                                               aTTCGAGCCGGACATTTACATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTGCGCGTTAcc     >  1:212318/1‑69 (MQ=255)                                                ttCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTACGCCGCTGatcatc                   <  1:1801189/54‑1 (MQ=255)                                                ttCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTACGCCGCTGatcatc                   <  1:2152244/54‑1 (MQ=255)                                                ttCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTACGCCGCTGatcatc                   <  1:204798/54‑1 (MQ=255)                                                ttCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTACGCCGCTGagcatc                   <  1:1845437/54‑1 (MQ=255)                                                ttCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTGCGCGTTACCgag  <  1:1969240/71‑1 (MQ=255)                                                             tttCCATATCGATTGTGCCGCCACGTACGCCGCTGATCATCGCCTgcgc           <  1:285313/49‑1 (MQ=255)                                                             tttCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTgcgc           <  1:642165/49‑1 (MQ=255)                                                             tttCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTgcgc           <  1:79511/49‑1 (MQ=255)                                                             tttCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTgcgc           <  1:3305903/49‑1 (MQ=255)                                                             tttCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTgcgc           <  1:3173916/49‑1 (MQ=255)                                                             tttCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTgcgc           <  1:3144594/49‑1 (MQ=255)                                                             tttCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTgcgc           <  1:2800344/49‑1 (MQ=255)                                                             tttCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTgcgc           <  1:2551512/49‑1 (MQ=255)                                                               tCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTGCGCGt         >  1:1354850/1‑49 (MQ=255)                                                               tCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTGCGCGt         >  1:1293281/1‑49 (MQ=255)                                                               tCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTGCGCGt         >  1:587825/1‑49 (MQ=255)                                                               tCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTGCGCGt         >  1:1876114/1‑49 (MQ=255)                                                               tCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTGCGCGt         >  1:2098620/1‑49 (MQ=255)                                                               tCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTGCGCGt         >  1:3025235/1‑49 (MQ=255)                                                               tCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTGCGCGt         >  1:3015159/1‑49 (MQ=255)                                                                   |                                                    AAGGGACATCAAGCAAGTTCATCACTGGTGATAACCCGGCAAAGTTATTCGAGCCGGACATTTCCATATCGATGGTGCCGCCACGTACGCCGCTGATCATCGCCTGCGCGTTACCGAG  >  W3110S.gb/3894002‑3894119

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