Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,799,697	A→G	11.9%	K395K (AAA→AAG)	yicM →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,799,697	0	A	G	11.9%	99.3 / 6.5	42	K395K (AAA→AAG)	yicM	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (15/22);  new base G (3/2);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

AACTTTTGCATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATCGGCGGTTAT  >  W3110S.gb/3799628‑3799763                                                                      |                                                                   aaCTTTTGCATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAg                                                                   >  1:3580530/1‑71 (MQ=255)   cTTTTGCATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGcc                                                                 <  1:3799612/71‑1 (MQ=255)     tttGCATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGc                                                                  >  1:665013/1‑68 (MQ=255)     tttGCATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGc                                                                  >  1:4110799/1‑68 (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:2238825/60‑1 (MQ=255)                         ttGGTCAACGTGGATCACTCGCTCATTGGCCGATCAGGCAGAAAAGGCTGGATCTATTCAg                                                     <  1:3079009/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:1876901/70‑1 (MQ=255)                                     gATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCt                                >  1:1069836/1‑70 (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: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)                                                                      |                                                                   AACTTTTGCATTGGTTCCCGTCGGCTGGTCAACGTGGATCACCCGCTCGCTGGCCGATCAGGCAGAAAAAGCCGGGTCTATTCAGGTGGCGGTTATTCAGCTTGCTAATACCTGTGGCGCGGCAATCGGCGGTTAT  >  W3110S.gb/3799628‑3799763

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