Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,803,912	G→T	27.2%	H125N (CAT→AAT)	thrS ←	threonyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,803,912	0	G	T	27.2%	74.1 / 13.7	33	H125N (CAT→AAT)	thrS	threonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (10/14);  new base T (3/6);  total (13/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATC  >  W3110S.gb/1803855‑1803966                                                          |                                                       cTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCg                                           <  1:2123441/71‑1 (MQ=255)   tCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAg                                         >  1:423459/1‑71 (MQ=255)    cGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAg                                         <  1:609601/70‑1 (MQ=255)    cGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAg                                         <  1:1751085/70‑1 (MQ=255)    cGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAg                                         <  1:3934266/70‑1 (MQ=255)    cGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAg                                         <  1:865181/70‑1 (MQ=255)    cGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAg                                         <  1:2876041/70‑1 (MQ=255)    cGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAg                                         <  1:278940/70‑1 (MQ=255)         cAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTc                                   <  1:38866/71‑1 (MQ=255)           gCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATTCATCCGCTTCTCGAGTGCTTCGa                                 <  1:3584785/71‑1 (MQ=255)           gCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATTCATCCGCTTCTCGAGTGCTTCGa                                 <  1:3802862/71‑1 (MQ=255)           gCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATTCATCCGCTTCTCGAGTGCTTCGa                                 <  1:3011043/71‑1 (MQ=255)           gCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATTCATCCGCTTCTCGAGTGCTTCGa                                 <  1:1765261/71‑1 (MQ=255)               aCTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACa                               <  1:1862103/69‑1 (MQ=255)                   tcttctTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATTCATCCGCTTCTCGAGTGCTTCGACATCTTcc                         >  1:3406245/1‑71 (MQ=255)                   tcttctTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTcc                         >  1:782212/1‑71 (MQ=255)                   tcttctTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTcc                         >  1:3960668/1‑71 (MQ=255)                   tcttctTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTc                          >  1:2775309/1‑70 (MQ=255)                    cttcttAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTcc                         <  1:4019851/70‑1 (MQ=255)                     ttcttAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACAt                              >  1:1847855/1‑64 (MQ=255)                                gtcgtAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAg                                         <  1:181513/42‑1 (MQ=255)                                   gtAGTTTTTCTCAGCAAGCTCATTCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGt         >  1:2018913/1‑71 (MQ=255)                                   gtAGTTTTTCTCAGCAAGCTAATTCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGt         >  1:1419606/1‑71 (MQ=255)                                      gTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCa       <  1:53550/70‑1 (MQ=255)                                      gTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCa       <  1:656386/70‑1 (MQ=255)                                      gTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCa       <  1:1568946/70‑1 (MQ=255)                                        ttttCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGCCATCTTCCTGGGTTAACGTGCGGTc        >  1:3135907/1‑67 (MQ=255)                                        ttttCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTc        >  1:3521552/1‑67 (MQ=255)                                        ttttCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTc        >  1:3201782/1‑67 (MQ=255)                                         tttCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGa    >  1:255552/1‑70 (MQ=255)                                           tctcAGCAAGCTCATTCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATc  <  1:4020780/70‑1 (MQ=255)                                           tctcAGCAAGCTCATTCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATc  <  1:3923450/70‑1 (MQ=255)                                                       cATGCATCCGCTTCTCGAGTGCTTCGACATCTTcc                         >  1:2730489/1‑35 (MQ=255)                                                          |                                                       CTTCGTGCCAGCTGACTTTCTTCTTAATGACGTCGTAGTTTTTCTCAGCAAGCTCATGCATCCGCTTCTCGAGTGCTTCGACATCTTCCTGGGTTAACGTGCGGTCAAGATC  >  W3110S.gb/1803855‑1803966

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