BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,425,365	C→A	12.9%	L189L (CTG→CTT)	yfcF ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,425,365	0	C	A	12.9%	171.2 / 9.3	54	L189L (CTG→CTT)	yfcF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (39/8); new base A (2/5); total (41/13)
Fisher's exact test for biased strand distribution p-value = 6.37e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

CGCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAGCGCCAGAT  >  W3110S.gb/2425302‑2425423

cGCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGt                                                      <  1:2394166/70‑1 (MQ=255)

cGCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGt                                                      <  1:2309707/70‑1 (MQ=255)

cGCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGt                                                      <  1:1174513/70‑1 (MQ=255)

cGCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGt                                                      <  1:1895767/70‑1 (MQ=255)

cGCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGt                                                      <  1:2689009/70‑1 (MQ=255)

 gCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:284149/1‑71 (MQ=255)

 gCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:2432184/1‑71 (MQ=255)

 gCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:2521943/1‑71 (MQ=255)

 gCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:1998145/1‑71 (MQ=255)

 gCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:2650549/1‑71 (MQ=255)

 gCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:1541214/1‑71 (MQ=255)

 gCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:1500629/1‑71 (MQ=255)

 gCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:1454834/1‑71 (MQ=255)

 gCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:1407905/1‑71 (MQ=255)

  ccGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTCt                                                   >  1:2680237/1‑71 (MQ=255)

           cAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCCTAATCCACCAGGCGTTc                                                    >  1:3184218/1‑61 (MQ=255)

           cAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:1686231/1‑61 (MQ=255)

           cAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:2553133/1‑61 (MQ=255)

           cAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:1083137/1‑61 (MQ=255)

           cAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:7795/1‑61 (MQ=255)

           cAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:2444191/1‑61 (MQ=255)

           cAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTc                                                    >  1:2461093/1‑61 (MQ=255)

              tAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACAAGGCGTTCTgg                                                 <  1:148498/61‑1 (MQ=255)

              tAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACAAGGCGTTCTgg                                                 <  1:2105221/61‑1 (MQ=255)

              tAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACAAGGCGTTCTgg                                                 <  1:2094836/61‑1 (MQ=255)

              tAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACAAGGCGTTCTgg                                                 <  1:659885/61‑1 (MQ=255)

              tAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACAAGGCGTTCTgg                                                 <  1:137982/61‑1 (MQ=255)

                aaaCGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACAAGGCGTTCTGGCa                                               >  1:2088201/1‑61 (MQ=255)

                aaaCGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACAAGGCGTTCTGGCa                                               >  1:3352948/1‑61 (MQ=255)

                    gCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGt                                                      <  1:768606/50‑1 (MQ=255)

                        gACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATgtagt                              <  1:3131980/70‑1 (MQ=255)

                         aCAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGt                                                      <  1:2068357/45‑1 (MQ=255)

                            gACGCTCGCTGCCACTGGAATGTTGCATAATccacc                                                            <  1:2398463/36‑1 (MQ=255)

                            gACGCTCGCTGCCACTGGAATGTTGCATAATccacc                                                            <  1:2900803/36‑1 (MQ=255)

                            gACGCTCGCTGCCACTGGAATGTTGCATAATccacc                                                            <  1:2923089/36‑1 (MQ=255)

                            gACGCTCGCTGCCACTGGAATGTTGCATAATccacc                                                            <  1:3032379/36‑1 (MQ=255)

                            gACGCTCGCTGCCACTGGAATGTTGCATAATccacc                                                            <  1:2062671/36‑1 (MQ=255)

                            gACGCTCGCTGCCACTGGAATGTTGCATAATccacc                                                            <  1:463984/36‑1 (MQ=255)

                            gACGCTCGCTGCCACTGGAATGTTGCATAATccacc                                                            <  1:524764/36‑1 (MQ=255)

                            gACGCTCGCTGCCACTGGAATGTTGCATAATccacc                                                            <  1:1517570/36‑1 (MQ=255)

                            gACGCTCGCTGCCACTGGAATGTTGCATAATccacc                                                            <  1:98994/36‑1 (MQ=255)

                                            ggAATGTTGGATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTAc                            >  1:3049793/1‑52 (MQ=255)

                                            ggAATGTTGCATAATCCACCAGGCGTTCTGGCACCTCATCGCCAt                                   >  1:557753/1‑45 (MQ=255)

                                            ggAATGTTGCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTAc                            >  1:300056/1‑52 (MQ=255)

                                            ggAATGTTGCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTAc                            >  1:2305312/1‑52 (MQ=255)

                                            ggAATGTTGCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTAc                            >  1:2873324/1‑52 (MQ=255)

                                            ggAATGTTGCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTAc                            >  1:1073861/1‑52 (MQ=255)

                                            ggAATGTTGCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTAc                            >  1:1661362/1‑52 (MQ=255)

                                            ggAATGTTGCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTAc                            >  1:2009361/1‑52 (MQ=255)

                                            ggAATGTTGCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAgc         >  1:3056929/1‑71 (MQ=255)

                                            ggAATGTTGCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAgc         >  1:1920232/1‑71 (MQ=255)

                                            ggAATGTTGCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAg          >  1:100245/1‑70 (MQ=255)

                                                    gCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGtt                   >  1:381378/1‑53 (MQ=255)

                                                    gCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAGCGCCAGa   >  1:2088622/1‑69 (MQ=255)

                                                    gCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAGCGCCAGAt  >  1:1980855/1‑70 (MQ=255)

                                                    gCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAGCGCCAGAt  >  1:1990869/1‑70 (MQ=255)

                                                    gCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAGCGCCAGAt  >  1:488708/1‑70 (MQ=255)

                                                    gCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAGCGCCAGAt  >  1:1614228/1‑70 (MQ=255)

                                                    gCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAGCGCCAGAt  >  1:1353841/1‑70 (MQ=255)

                                                    gCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAGCGCCAGAt  >  1:1189187/1‑70 (MQ=255)

                                                    gCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAGCGCCAGAt  >  1:1161963/1‑70 (MQ=255)

                                                    gCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAGCGCCAGAt  >  1:2127167/1‑70 (MQ=255)

                                                    gCAGAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAGCGCCAGAt  >  1:1764536/1‑70 (MQ=255)

CGCCGAAAGTGCAATAAAACGCTGGACAGACGCTCGCTGCCACTGGAATGTTGCATAATCCACCAGGCGTTCTGGCACCTCATCGCCATGTAGTACCAGGCGGTTAATCATTAGCGCCAGAT  >  W3110S.gb/2425302‑2425423

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