BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,183,814	T→A	27.1%	V129E (GTG→GAG)	yegU →	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,183,814	0	T	A	27.1%	80.2 / 36.1	52	V129E (GTG→GAG)	yegU	predicted hydrolase
Reads supporting (aligned to +/- strand): ref base T (15/22); new base A (5/9); total (21/31)
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.73e-01

CCTCGAAGATTGCGCTTAACGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTC  >  W3110S.gb/2183748‑2183880

ccTCGAAGATTGCGCTTAACGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGAcc                                                                <  1:377272/71‑1 (MQ=255)

          tGCGCTTAACGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAg                                                      >  1:2786660/1‑71 (MQ=255)

          tGCGCTTAACGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAg                                                      >  1:2961496/1‑71 (MQ=255)

          tGCGCTTAACGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGAGTGACCAACGGCGCAg                                                      >  1:2487231/1‑71 (MQ=255)

            cgcTTAACGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGc                                                     <  1:2036290/70‑1 (MQ=255)

            cgcTTAACGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGAGACCAACGGCGCAGCg                                                    <  1:2869782/71‑1 (MQ=255)

            cgcTTAACGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGAGACCAACGGCGCAGCg                                                    <  1:2250912/71‑1 (MQ=255)

               ttAACGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGAGACCAACGGCGCAGCGATg                                                 >  1:3642746/1‑71 (MQ=255)

                 aaCGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATgc                                                >  1:2103301/1‑70 (MQ=255)

                 aaCGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATgc                                                >  1:3333928/1‑70 (MQ=255)

                 aaCGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATg                                                 <  1:1980511/69‑1 (MQ=255)

                  aCGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATgcg                                               <  1:3114728/70‑1 (MQ=255)

                  aCGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATgcg                                               <  1:734477/70‑1 (MQ=255)

                         tCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATgcg                                               <  1:2775705/63‑1 (MQ=255)

                          cgcgACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGc                                       >  1:906900/1‑70 (MQ=255)

                          cgcgACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGc                                       >  1:1223113/1‑70 (MQ=255)

                          cgcgACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGAGACCAACGGCGCAGCGATGCGCGTCTCGc                                       >  1:2432460/1‑70 (MQ=255)

                          cgcgACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGAGACCAACGGCGCAGCGATGCGCGTCTCGc                                       >  1:4065772/1‑70 (MQ=255)

                          cgcgACGGTAAACCCGTTGCTGAACTGGAAAACAAAGGCGAGACCAACGGCGCAGCGATGCGCGTCTCGc                                       >  1:2979213/1‑70 (MQ=255)

                             gACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGAGACCAACGGCGCAGCGATGCGCGTCTCGCCAtt                                   <  1:3646549/71‑1 (MQ=255)

                              aCGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTa                                  <  1:3097414/71‑1 (MQ=255)

                              aCGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTa                                  <  1:2466514/71‑1 (MQ=255)

                               cGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAg                                 >  1:1522537/1‑71 (MQ=255)

                                   aaaCCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGgttg                             >  1:2316815/1‑71 (MQ=255)

                                   aaaCCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGgttg                             >  1:2978200/1‑71 (MQ=255)

                                   aaaCCCGTTGCTGAACTGGAAAACAACGGCGCGACCAACGGCGCAGCGATGCGCGTCTCGACATTAGgttg                             >  1:2394229/1‑71 (MQ=255)

                                   aaaCCCGTTGCTGAACTGGAAAACAACGGCGAGACAAACGGCGCAGCGATGCGCGTCTCGCCATTAGgttg                             >  1:3884343/1‑71 (MQ=255)

                                    aaCCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGgttg                             <  1:2861669/70‑1 (MQ=255)

                                    aaCCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGgttg                             <  1:3426503/70‑1 (MQ=255)

                                    aaCCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGgttg                             <  1:579815/70‑1 (MQ=255)

                                          ttGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGc                                       <  1:1562180/54‑1 (MQ=255)

                                          ttGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGc                                       <  1:86974/54‑1 (MQ=255)

                                            gCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTcc                    >  1:3970753/1‑71 (MQ=255)

                                            gCTGAACTGGAAAACAACGGCGTGACCAACGGCGAAGCGa                                                   >  1:3653315/1‑40 (MQ=255)

                                                  cTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGc                        <  1:2992121/61‑1 (MQ=255)

                                                    ggAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGt             >  1:1948350/1‑70 (MQ=255)

                                                     gAAAACAACGGCGAGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGt             <  1:32071/69‑1 (MQ=255)

                                                       aaaCAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATg         <  1:2507097/71‑1 (MQ=255)

                                                        aacaacGGCGTGACCAACGGCGCAGCGATTCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATg         <  1:2411860/70‑1 (MQ=255)

                                                         acaacGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGtt       <  1:1532432/71‑1 (MQ=255)

                                                         acaacGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGtt       <  1:2404532/71‑1 (MQ=255)

                                                         acaacGGCGAGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGtt       <  1:3283858/71‑1 (MQ=255)

                                                         acaacGGCGAGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGtt       <  1:1928682/71‑1 (MQ=255)

                                                         acaacGGCGAGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGtt       <  1:1832014/71‑1 (MQ=255)

                                                         acaacGGCGAGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGtt       <  1:2894274/71‑1 (MQ=255)

                                                          caacGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGtt       >  1:2797731/1‑70 (MQ=255)

                                                           aacGGCGAGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGtt       <  1:3223591/69‑1 (MQ=255)

                                                             cGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCt                       <  1:3388465/51‑1 (MQ=255)

                                                             cGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCt                       <  1:3383886/51‑1 (MQ=255)

                                                             cGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCt                       <  1:405957/51‑1 (MQ=255)

                                                             cGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCt                       <  1:2772329/51‑1 (MQ=255)

                                                              ggCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTc  >  1:2446193/1‑71 (MQ=255)

CCTCGAAGATTGCGCTTAACGCCATTCGCGACGGTAAACCCGTTGCTGAACTGGAAAACAACGGCGTGACCAACGGCGCAGCGATGCGCGTCTCGCCATTAGGTTGTTTGCTTCCGGCGCGTGATGTTGATTC  >  W3110S.gb/2183748‑2183880

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