Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,015,328	G→A	40.0%	H238Y (CAT→TAT)	yhhI ←	predicted transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,015,328	0	G	A	40.0%	24.1 / 24.3	25	H238Y (CAT→TAT)	yhhI	predicted transposase
Reads supporting (aligned to +/- strand):  ref base G (9/6);  new base A (7/3);  total (16/9)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACTGTCATGCTCTGGATTATTTAA  >  W3110S.gb/4015298‑4015397                               |                                                                      tAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGAtt                               >  1:3753907/1‑71 (MQ=34)   aGTTCATCAGGGACATCGCAAACAATATAAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGc                             >  1:3448386/1‑71 (MQ=21)     ttCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGc                             <  1:3665298/69‑1 (MQ=35)       catcaGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTa                          <  1:2981478/70‑1 (MQ=35)       catcaGGGACATCGCAAACAATATAAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTaa                         >  1:164495/1‑71 (MQ=21)       catcaGGGACATCGCAAACAATATAAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTa                          <  1:3114651/70‑1 (MQ=25)        atcaGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTaa                         <  1:914261/70‑1 (MQ=35)         tcaGGGACATCGCAAACAATATAAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACt                       >  1:1938433/1‑71 (MQ=21)           aGGGACATCGCAAACAATATAAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACTGt                     >  1:5057098/1‑71 (MQ=21)             ggACATCGCAAACAATATGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACTGtaa                   <  1:3718714/71‑3 (MQ=21)             ggACATCGCAAACAATATGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACTGTCa                   <  1:1722315/71‑1 (MQ=34)                  tCGCAAACAATATAAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTg                                  <  1:4379946/51‑1 (MQ=9)                    gCAAACAATGTAAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGAt                                >  1:4675888/1‑51 (MQ=11)                    gCAAACAATATGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACTGTCATGCTCTg            >  1:1672663/1‑71 (MQ=34)                    gCAAACAATATGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACTGTCATGCTCTg            >  1:4683237/1‑71 (MQ=34)                    gCAAACAATATGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACTGTCATGCTCTg            >  1:4619482/1‑71 (MQ=34)                    gCAAACAATATGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTAACTGATTGCGTAACTGTCATGCTCTg            >  1:190132/1‑71 (MQ=21)                    gCAAACAATATAAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCCCTGATTGCGTAACTGTCATGCTATg            <  1:2614148/71‑1 (MQ=11)                    gCAAACAATATAAAGACGGATTTCTTCTCTGCC‑TGACTCTTTTCACTGATTGCGTAACTGTCATGCTCTg            >  1:1297376/1‑70 (MQ=11)                     cAAACAATATGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACTGTCATGCTCTg            >  1:4421853/1‑70 (MQ=35)                       aaCAATATGAAGACGGATCTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACTGTCATGCTCTGGat         >  1:85790/1‑71 (MQ=255)                         cAATATAAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACTGTCATGCTCTGGatta       >  1:1856634/1‑71 (MQ=21)                             atGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGAtt                               >  1:4587018/1‑43 (MQ=38)                             atGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACTGTCATGCTCTGGATTAttt    <  1:4400070/70‑1 (MQ=35)                              tGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACTGTCATGCTCTGGATTATTTaa  >  1:3975462/1‑71 (MQ=34)                               |                                                                      TAAGTTCATCAGGGACATCGCAAACAATATGAAGACGGATTTCTTCTCTGCCGTGACTCTTTTCACTGATTGCGTAACTGTCATGCTCTGGATTATTTAA  >  W3110S.gb/4015298‑4015397

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