Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,799,703	G→A	11.6%	G397G (GGG→GGA)	yicM →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,799,703	0	G	A	11.6%	116.7 / 6.0	43	G397G (GGG→GGA)	yicM	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (14/24);  new base A (3/2);  total (17/26)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATCGGCGGTTATGCGCT  >  W3110S.gb/3799636‑3799768                                                                    |                                                                  cATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTc                                                                 <  1:2129597/70‑1 (MQ=255) cATTGGTTCCCGTCGGCTGGGCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTc                                                                 <  1:3159034/70‑1 (MQ=255)      gTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCa                                                           <  1:2224550/71‑1 (MQ=255)               ggCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGt                                                  <  1:4324399/71‑1 (MQ=255)               ggCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGt                                                  <  1:1309878/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:1690663/66‑1 (MQ=255)                         cGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTc                                             <  1:1860058/66‑1 (MQ=255)                           tGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATT‑AGGTGGCGGTTATTCAGCTTGCt                                     <  1:3534901/71‑1 (MQ=255)                             gATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCt                                     <  1:2622604/70‑1 (MQ=255)                             gATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCt                                     >  1:1069836/1‑70 (MQ=255)                             gATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCt                                     <  1:1876901/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:2196084/71‑1 (MQ=255)                                                          aaaaaGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATCGGCGGTTAt       >  1:2158322/1‑71 (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)                                                                    gTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTTGCGCGGCAATcg               <  1:1044666/53‑1 (MQ=255)                                                                    gTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATcg               <  1:2220677/53‑1 (MQ=255)                                                                    |                                                                  CATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATCGGCGGTTATGCGCT  >  W3110S.gb/3799636‑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: ‑