Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	131,249	A→G	14.3%	T4T (ACT→ACC)	yacH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	131,249	0	A	G	14.3%	102.1 / 9.6	42	T4T (ACT→ACC)	yacH	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (17/19);  new base G (4/2);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAG  >  W3110S.gb/131192‑131309                                                          |                                                             ggcggcACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATaa                                                 >  1:1118646/1‑71 (MQ=255)  gcggcACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATaa                                                 <  1:2862234/70‑1 (MQ=255)  gcggcACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATaa                                                 <  1:800423/70‑1 (MQ=255)    ggcACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAAtt                                               >  1:3123096/1‑70 (MQ=255)                   aCTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAGGTCATTTTCATAATTCGGATCTCAAGGAAAt                               >  1:1120850/1‑71 (MQ=255)                   aCTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGAt                                          >  1:91628/1‑60 (MQ=255)                   aCTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGAt                                          >  1:3067791/1‑60 (MQ=255)                      gCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGaaa                                >  1:629035/1‑67 (MQ=255)                      gCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCg                             >  1:1116352/1‑70 (MQ=255)                      gCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTAAAGGAAATCg                             >  1:3645397/1‑70 (MQ=255)                       cAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGc                            <  1:1621202/70‑1 (MQ=255)                       cAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGc                            <  1:3378956/70‑1 (MQ=255)                          aTTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCAt                                                   <  1:1261830/44‑1 (MQ=255)                          aTTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGc                            <  1:1414028/67‑1 (MQ=255)                          aTTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGc                            <  1:2230024/67‑1 (MQ=255)                          aTTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGc                            <  1:2296413/67‑1 (MQ=255)                             aGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCg                                             >  1:804376/1‑47 (MQ=255)                                 ccAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCg                 <  1:954533/71‑1 (MQ=255)                                  cAGCACATGGGGTTTAAACGGCAAGGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCg                 <  1:3220391/70‑1 (MQ=255)                                  cAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGa                <  1:1270292/71‑1 (MQ=255)                                   aGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGaa               >  1:3956555/1‑71 (MQ=255)                                     cacaTGGGGTTTAAACGGCAAGGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACt             >  1:255228/1‑71 (MQ=255)                                      acaTGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTg            >  1:3046669/1‑71 (MQ=255)                                        aTGGGGTTTAAACGGCAAGGTCATTTTCATAATTCGGATCTCAAGGAAATCGc                            <  1:3585190/53‑1 (MQ=255)                                        aTGGGGTTTAAACGGCAAAGTCATTTTCATAATTCgg                                            <  1:3821764/37‑1 (MQ=255)                                        aTGGGGTTTAAACGGCAAAGTCATTTTCATAATTCgg                                            <  1:1233373/37‑1 (MQ=255)                                        aTGGGGTTTAAACGGCAAAGTCATTTTCATAATTCgg                                            <  1:1290986/37‑1 (MQ=255)                                         tGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCa         >  1:1333490/1‑71 (MQ=255)                                          ggggTTTAAACGGCAAGGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCAc        >  1:454135/1‑71 (MQ=255)                                          ggggTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGc                            <  1:1207534/51‑1 (MQ=255)                                              tttAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGc                            <  1:2396900/47‑1 (MQ=255)                                              tttAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGc                            <  1:2217009/47‑1 (MQ=255)                                              tttAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGc                            <  1:1267771/47‑1 (MQ=255)                                              tttAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCACTGGTCGGCGAACTGCCACCCGc    <  1:3670232/71‑1 (MQ=255)                                                tAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTc                     >  1:3449211/1‑52 (MQ=255)                                                tAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGaa               >  1:3370656/1‑58 (MQ=255)                                                tAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGaa               >  1:2852288/1‑58 (MQ=255)                                                tAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGaa               >  1:2345595/1‑58 (MQ=255)                                                tAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCa   <  1:2886094/70‑1 (MQ=255)                                                tAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAg  >  1:1761142/1‑71 (MQ=255)                                                tAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAg  >  1:1458607/1‑71 (MQ=255)                                                  aaCGGCAAGGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAAc              >  1:684295/1‑57 (MQ=255)                                                          |                                                             GGCGGCACAAAGCCCGGCACTGCAAATTAGCGCCAGCACATGGGGTTTAAACGGCAAAGTCATTTTCATAATTCGGATCTCAAGGAAATCGCAATGGTCGGCGAACTGCCACCCGCAG  >  W3110S.gb/131192‑131309

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