Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,130,879:1	+A	15.4%	coding (275/1227 nt)	yhfW →	predicted mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,130,879	1	.	A	15.4%	115.1 / 21.8	39	coding (275/1227 nt)	yhfW	predicted mutase
Reads supporting (aligned to +/- strand):  ref base . (12/21);  new base A (4/2);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTGGCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAGCAGGCATT  >  W3110S.gb/4130821‑4130948                                                            |                                                                      ggCGTGGCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGccc                                                           <  1:3968452/71‑1 (MQ=255)   cGTGGCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTT‑Agg                                                                    >  1:1206369/1‑60 (MQ=255)   cGTGGCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTT‑Agg                                                                    >  1:724556/1‑60 (MQ=255)       gCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAAGGCACtcgc                                                             >  1:1145049/1‑64 (MQ=255)       gCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAAGGCACGCGCCCGTTAc                                                      >  1:4066532/1‑71 (MQ=255)        cAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTAcc                                                     >  1:2149468/1‑70 (MQ=255)           agCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAAGGCACGCGCCCGt                                                         >  1:2843713/1‑64 (MQ=255)              tGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAAGGCACGCTcc                                                            <  1:3115121/58‑1 (MQ=255)              tGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCt                                               >  1:537777/1‑70 (MQ=255)              tGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTg                                              >  1:494077/1‑71 (MQ=255)               gCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGc                                             >  1:4049065/1‑71 (MQ=255)                cAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCg                                            >  1:88312/1‑71 (MQ=255)                cAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCg                                            <  1:3043488/71‑1 (MQ=255)                cAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCg                                            <  1:918931/71‑1 (MQ=255)                 aaCATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTTAAGGCACGCGCCCGTTACCGCCGCTGCg                                            <  1:3241171/71‑1 (MQ=255)                                aTACCTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGtt                                                        <  1:713512/43‑1 (MQ=255)                                aTACCTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGtt                                                        <  1:93026/43‑1 (MQ=255)                                   ccTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgtg                         <  1:231225/71‑1 (MQ=255)                                   ccTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgtg                         <  1:722196/71‑1 (MQ=255)                                   ccTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgtg                         <  1:2666874/71‑1 (MQ=255)                                   ccTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgtg                         <  1:2917889/71‑1 (MQ=255)                                   ccTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgtg                         <  1:808389/71‑1 (MQ=255)                                   ccTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgt                          <  1:868748/70‑1 (MQ=255)                                   ccTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgt                          <  1:939539/70‑1 (MQ=255)                                   ccTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgt                          <  1:3406845/70‑1 (MQ=255)                                   ccTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgt                          <  1:1272248/70‑1 (MQ=255)                                    cTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGAtgt                          <  1:4081856/69‑1 (MQ=255)                                           gggCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCg                 <  1:1100273/71‑1 (MQ=255)                                            ggCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCgt                <  1:2199153/71‑1 (MQ=255)                                            ggCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCgt                <  1:1346311/71‑1 (MQ=255)                                                  aGGAAATTTTAAGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGtt             >  1:1737887/1‑69 (MQ=255)                                                  aGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCgtgt              >  1:835542/1‑67 (MQ=255)                                                  aGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGa           >  1:3250214/1‑70 (MQ=255)                                                   ggAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAg          <  1:2393158/70‑1 (MQ=255)                                                    gAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCgt                <  1:61442/63‑1 (MQ=255)                                                    gAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCgt                <  1:659176/63‑1 (MQ=255)                                                     aaaTTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGcc                                      >  1:354948/1‑40 (MQ=255)                                                     aaaTTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAg          >  1:3882708/1‑68 (MQ=255)                                                     aaaTTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAg          >  1:1514783/1‑68 (MQ=255)                                                      atttta‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAGCAgg      >  1:1716367/3‑71 (MQ=255)                                                      atttta‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAGCAgg      >  1:1984491/3‑71 (MQ=255)                                                        tttt‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCgtg               <  1:2617571/60‑1 (MQ=255)                                                         ttt‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAGCAGGc     <  1:1648307/69‑1 (MQ=255)                                                           t‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAGCAGGCAtt  <  1:1791213/70‑1 (MQ=255)                                                            |                                                                      GGCGTGGCAGAGCTGCAACATGAAGGTGGCGATACCTTTATGGGGCATCAGGAAATTTT‑AGGCACGCGCCCGTTACCGCCGCTGCGGATGCCTTTTCGCGATGTGATTGACCGTGTTGAGCAGGCATT  >  W3110S.gb/4130821‑4130948

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