Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,799,700	C→T	11.3%	A396A (GCC→GCT)	yicM →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,799,700	0	C	T	11.3%	119.0 / 5.6	44	A396A (GCC→GCT)	yicM	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (15/24);  new base T (3/2);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CTTTTGCATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATCGGCGGTTATGCGCT  >  W3110S.gb/3799630‑3799768                                                                       |                                                                     cTTTTGCATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGcc                                                                      <  1:3799612/71‑1 (MQ=255)   tttGCATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCgg                                                                    >  1:881229/1‑71 (MQ=255)       cATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTc                                                                 <  1:2129597/70‑1 (MQ=255)       cATTGGTTCCCGTCGGCTGGGCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTc                                                                 <  1:3159034/70‑1 (MQ=255)            gTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCa                                                           <  1:2224550/71‑1 (MQ=255)                     ggCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGt                                                  <  1:1309878/71‑1 (MQ=255)                     ggCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGt                                                  <  1:4324399/71‑1 (MQ=255)                      gCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGtt                                                 <  1:906901/71‑1 (MQ=255)                       ttGGTCAACGTGGATCACTCGCTCATTGGCCGATCAGGCAGAAAAGGCTGGATCTATTCAg                                                          <  1:3079009/60‑1 (MQ=255)                       ttGGTCAACGTGGATCACTCGCTCATTGGCCGATCAGGCAGAAAAGGCTGGATCTATTCAg                                                          <  1:2238825/60‑1 (MQ=255)                          gTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTAtt                                              >  1:4042757/1‑70 (MQ=255)                             aaCGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAg                                                          >  1:3641384/1‑55 (MQ=255)                             aaCGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGc                                          >  1:858813/1‑71 (MQ=255)                               cGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTc                                             <  1:1860058/66‑1 (MQ=255)                               cGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTc                                             <  1:1690663/66‑1 (MQ=255)                                 tGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATT‑AGGTGGCGGTTATTCAGCTTGCt                                     <  1:3534901/71‑1 (MQ=255)                                   gATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCt                                     >  1:1069836/1‑70 (MQ=255)                                   gATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCt                                     <  1:1876901/70‑1 (MQ=255)                                   gATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCt                                     <  1:2622604/70‑1 (MQ=255)                                      cACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATa                                 >  1:1391790/1‑71 (MQ=255)                                       aCCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATAc                                <  1:506847/71‑1 (MQ=255)                                        cccGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCa                                            <  1:2031256/58‑1 (MQ=255)                                         ccGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTAtt                                                             <  1:1051551/40‑1 (MQ=255)                                          cgctcgctGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTaa                                   >  1:3748134/1‑65 (MQ=255)                                            ctcgctGGCCGATCAGGCAGAAAAAGCCGGGTCTa                                                               <  1:4168180/35‑1 (MQ=255)                                            ctcgctGGCCGATCAGGCAGAAAAAGCCGGGTCTATTc                                                            <  1:864036/38‑1 (MQ=255)                                                ttGGCCGATCAGGCAGAAAAGGCTGGATCTATTCAGGTGGCCGTTATTCAGCTTGCTAACACCTGCggtg                        >  1:1592268/2‑68 (MQ=255)                                                ttGGCCGATCAGGCAGAAAAGGCTGGATCTATTCAGGTGGCCGTTATTCAGCTTGCTAACACCTGCGgtgc                       >  1:3291195/2‑71 (MQ=255)                                                ttGGCCGATCAGGCAGAAAAGGCTGGATCTATTCAGGTGGCCGTTATTCAGCTTGCTAACACATGCggtg                        >  1:1414623/2‑68 (MQ=255)                                                ctGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGgcgc                       >  1:4171094/1‑71 (MQ=255)                                                 tGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGgcgc                       <  1:2466345/70‑1 (MQ=255)                                                   gCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTc                                             >  1:2765324/1‑46 (MQ=255)                                                    ccGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCa                   <  1:4397394/71‑1 (MQ=255)                                                     cGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCaa                  >  1:1364897/1‑71 (MQ=255)                                                         cAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATcgg              >  1:2367112/1‑71 (MQ=255)                                                         cAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATcg               >  1:2141756/1‑70 (MQ=255)                                                          aGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATcgg              <  1:814982/70‑1 (MQ=255)                                                             cAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATcggcgg           <  1:1302382/70‑1 (MQ=255)                                                             cAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATCGGCGGt          <  1:3523890/71‑1 (MQ=255)                                                                aaaaaGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATCGGCGGTTAt       >  1:2158322/1‑71 (MQ=255)                                                                aaaaaGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATCGGCGGTTAt       <  1:2196084/71‑1 (MQ=255)                                                                 aaaaGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATcggc             >  1:4254054/1‑64 (MQ=255)                                                                 aaaaGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATcggc             >  1:357925/1‑64 (MQ=255)                                                                     gCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATCGGCGGTTATGCGCt  <  1:3919657/71‑1 (MQ=255)                                                                     gCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATCGGCGGTTATGCGCt  <  1:3501357/71‑1 (MQ=255)                                                                       |                                                                     CTTTTGCATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATCGGCGGTTATGCGCT  >  W3110S.gb/3799630‑3799768

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