Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,560,998	A→G	13.9%	F44S (TTT→TCT)	yihX ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,560,998	0	A	G	13.9%	104.2 / 10.7	43	F44S (TTT→TCT)	yihX	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (17/20);  new base G (2/4);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

GAGCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCTCCA  >  W3110S.gb/3560930‑3561062                                                                     |                                                                 gagCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGaaa                                                                 >  1:2181182/1‑70 (MQ=255) gagCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAAc                                                                <  1:1990574/71‑1 (MQ=255) gagCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAAc                                                                <  1:749618/71‑1 (MQ=255)   gCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCAACGCTCATGCTGATGAAACGc                                                              >  1:2454513/1‑71 (MQ=255)    ccATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGcc                                                             >  1:1014221/1‑71 (MQ=255)    ccATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGcc                                                             >  1:510742/1‑71 (MQ=255)      aTCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTc                                                           >  1:771181/1‑71 (MQ=255)      aTCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTc                                                           >  1:1596875/1‑71 (MQ=255)        ctcATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTccc                                                         >  1:298999/1‑71 (MQ=255)         tcATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGaaa                                                                 >  1:2004746/1‑62 (MQ=255)         tcATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGaaa                                                                 >  1:1952296/1‑62 (MQ=255)          cATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTccac                                                        >  1:2105927/1‑68 (MQ=255)                 caGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaa                                                <  1:1965638/71‑1 (MQ=255)                 caGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaa                                                <  1:413181/71‑1 (MQ=255)                   gcgcCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaa                                                >  1:1441920/1‑69 (MQ=255)                      ccTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAActc                                           <  1:2703009/71‑1 (MQ=255)                         ctGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACtcttc                                        >  1:1271020/1‑71 (MQ=255)                          tGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGaaaa                                              >  1:141426/1‑64 (MQ=255)                          tGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTcccc                                                        <  1:1778543/54‑1 (MQ=255)                          tGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaaaa                                              >  1:510209/1‑64 (MQ=255)                          tGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaa                                                <  1:1435149/62‑1 (MQ=255)                           gCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTcccc                                                        <  1:1581722/53‑1 (MQ=255)                            cGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaaaa                                              >  1:1920013/1‑62 (MQ=255)                            cGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaa                                                <  1:475303/60‑1 (MQ=255)                            cGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTcttctt                                      >  1:1110174/1‑70 (MQ=255)                               aCGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                  <  1:1601268/71‑1 (MQ=255)                                cGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaa                                                <  1:597003/56‑1 (MQ=255)                                cGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACtcttc                                        >  1:312436/1‑64 (MQ=255)                                 gCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGa                                <  1:1408113/71‑1 (MQ=255)                                  cTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCACCCATATGAAAACTCTTCTTAAGCGAt                               >  1:801117/1‑71 (MQ=255)                                     cgtcgCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCa                                                      <  1:520685/45‑1 (MQ=255)                                          cTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCgg                       <  1:37980/71‑1 (MQ=255)                                             aTTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCgg                       <  1:1522648/68‑1 (MQ=255)                                             aTTTCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCgg                       <  1:588747/68‑1 (MQ=255)                                             aTTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGACAGCGGaa                     <  1:1396601/70‑1 (MQ=255)                                              tttCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGaa                     <  1:172437/69‑1 (MQ=255)                                              tttCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaa                                                <  1:1485303/42‑1 (MQ=255)                                                tCCCCACGCTCATGCTGATGAGACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATAcg                 >  1:2236165/1‑71 (MQ=255)                                                tCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGaaaa                                              >  1:2113383/1‑42 (MQ=255)                                                  cccACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGaa                     <  1:1740267/65‑1 (MQ=255)                                                      cGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTa            >  1:1058017/1‑70 (MQ=255)                                                            tGCTTATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGc                                  <  1:1411482/42‑1 (MQ=255)                                                              cTGATGAAACGCCTCCCCCATATGAAAACTcttctt                                      <  1:552670/36‑1 (MQ=255)                                                              cTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCTcc   <  1:1083195/71‑1 (MQ=255)                                                                 atgaAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTaaa          >  1:1947795/1‑61 (MQ=255)                                                                 atgaAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCTCCa  <  1:1471329/69‑1 (MQ=255)                                                                     |                                                                 GAGCCATCTCATGACACAGCGCCTCTGCGAACGCTTCGTCGCTAATTTCCCCACGCTCATGCTGATGAAACGCCTCCCCCATATGAAAACTCTTCTTAAGCGATGCCAGCGGAATACGCGTTAAATCGCTCCA  >  W3110S.gb/3560930‑3561062

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