Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,742,822	A→G	80.4%	S400G (AGC→GGC)	ydhC →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,742,822	0	A	G	80.4%	101.7 / 18.7	46	S400G (AGC→GGC)	ydhC	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (2/7);  new base G (22/15);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 6.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

GGTTACATGATGCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATAGCGAATCACACTGATCTATATCGATATACTTATACTTAGGCTGCTAACAAAATTTTGTTGTATGATT  >  W3110S.gb/1742756‑1742889                                                                   |                                                                    ggTTACATGATGCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCg                                                                  <  1:568921/70‑1 (MQ=255) ggTTACATGATGCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCg                                                                  <  1:3594476/70‑1 (MQ=255)     aCATGATGCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATca                                                             >  1:2587974/1‑71 (MQ=255)         gatgCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCacac                                                          >  1:8240/1‑70 (MQ=255)         gatgCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACt                                                         >  1:3285453/1‑71 (MQ=255)         gatgCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACt                                                         >  1:705399/1‑71 (MQ=255)         gatgCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACt                                                         >  1:131201/1‑71 (MQ=255)         gatgCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACt                                                         >  1:3050855/1‑71 (MQ=255)         gatgCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACt                                                         >  1:3326926/1‑71 (MQ=255)         gatgCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACt                                                         >  1:3356955/1‑71 (MQ=255)         gatgCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACt                                                         >  1:2785675/1‑71 (MQ=255)         gatgCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACt                                                         >  1:43863/1‑71 (MQ=255)          atgCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTg                                                        <  1:1088075/71‑1 (MQ=255)          atgCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTg                                                        <  1:441103/71‑1 (MQ=255)          atgCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTg                                                        <  1:1720987/71‑1 (MQ=255)           tgCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGa                                                       >  1:3397357/1‑71 (MQ=255)              aaCGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCg                                                                  <  1:2279124/57‑1 (MQ=255)               aCGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCta                                                   <  1:1751327/71‑1 (MQ=255)               aCGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCta                                                   <  1:3077199/71‑1 (MQ=255)               aCGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCta                                                   <  1:3167943/71‑1 (MQ=255)               aCGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATAGCGAATCACACTGATCta                                                   <  1:2166266/71‑1 (MQ=255)               aCGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATAGCGAATCACACTGATCta                                                   <  1:1759818/71‑1 (MQ=255)               aCGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATAGCGAATCACACTGATCta                                                   <  1:3303875/71‑1 (MQ=255)                     tgAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCTATATCGa                                             <  1:2440800/71‑1 (MQ=255)                           agTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCTATATCGAtatact                                       <  1:3083495/71‑1 (MQ=255)                           agTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCTATATCGAtatact                                       <  1:1147767/71‑1 (MQ=255)                               ggCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCTATATCGAtatacttata                                   >  1:2116797/1‑71 (MQ=255)                               ggCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCTATATCGAtatacttata                                   >  1:3508917/1‑71 (MQ=255)                                gCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCTATATCGAtatacttata                                   <  1:3518528/70‑1 (MQ=255)                                          tCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCTATATCGAtatactta                                     >  1:945330/1‑58 (MQ=255)                                          tCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCTATATCGAtatactta                                     >  1:3640915/1‑58 (MQ=255)                                          tCATGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCTATATCGAtatactta                                     >  1:333703/1‑58 (MQ=255)                                            aTGGCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCTATATCGATATACTTATACTTAGGCTGCTaa                      >  1:1074312/1‑71 (MQ=255)                                               gCAATGCCGAAGTCGCTCATGGCGAATCACACTGATCTATATCGATATACTTATACTTAGGCTGCTAACaa                   <  1:701865/71‑1 (MQ=255)                                                     ccGAAGTCGCTCATGGCGAATCACACTGATCTATATCGATATACTTATACTTAGGCTGCTAACAAAAttt              >  1:3471985/1‑70 (MQ=255)                                                     ccGAAGTCGCTCATAGCGAATCACACTGATCTATATCGa                                             >  1:3451877/1‑39 (MQ=255)                                                     ccGAAGTCGCTCATAGCGAATCACACTGATCTATATCGATATACTTATACTTAGGCTGCTAACAAAAttt              >  1:1800532/1‑70 (MQ=255)                                                         aGTCGCTCATGGCGAATCACACTGATCTATATCGATATACTTATACTTAGGCTGCTAACAAAATTttgttg         <  1:230836/71‑1 (MQ=255)                                                            cGCTCATGGCGAATCACACTGATCTATATCGAtata                                         >  1:3515330/1‑36 (MQ=255)                                                             gCTCATGGCGAATCACACTGATCTATATCGATATACTTATACTTAGGCTGCTAACAAAATTTTGTTGTAt      >  1:2384990/1‑70 (MQ=255)                                                             gCTCATGGCGAATCACACTGATCTATATCGATATACTTATACTTAGGCTGCTAACAAAATTTTGTTGTAt      >  1:1164981/1‑70 (MQ=255)                                                                cATGGCGAATCACACTGATCTATATCGATATACTTATACTTAGGCTGCTAACAAAATTTTGTTGTATGAtt  >  1:820623/1‑71 (MQ=255)                                                                 aTAGCGAATCACACTGATCTATATCGAtatacttatac                                  <  1:2772556/38‑1 (MQ=255)                                                                 aTAGCGAATCACACTGATCTATATCGAtatacttatac                                  <  1:2925887/38‑1 (MQ=255)                                                                 aTAGCGAATCACACTGATCTATATCGAtatacttatac                                  <  1:2961790/38‑1 (MQ=255)                                                                 aTAGCGAATCACACTGATCTATATCGAtatacttatac                                  <  1:3077370/38‑1 (MQ=255)                                                                   |                                                                    GGTTACATGATGCAACGTTGTGAAGAAGTTGGCTGCCAGAATCATGGCAATGCCGAAGTCGCTCATAGCGAATCACACTGATCTATATCGATATACTTATACTTAGGCTGCTAACAAAATTTTGTTGTATGATT  >  W3110S.gb/1742756‑1742889

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