BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,374,788	G→T	19.9%	R165R (CGG→CGT)	yeiA →	predicted oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,374,788	0	G	T	19.9%	106.3 / 15.3	56	R165R (CGG→CGT)	yeiA	predicted oxidoreductase
Reads supporting (aligned to +/- strand): ref base G (21/23); new base T (2/9); total (23/33)
Fisher's exact test for biased strand distribution p-value = 9.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATCGGTGATATGTGCGAAGTGGC  >  minE/1374723‑1374850

gCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGt                                                             >  1:2604204/1‑69 (MQ=255)

gCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGt                                                             >  1:195158/1‑69 (MQ=255)

      ggTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACg                                                       >  1:493610/1‑69 (MQ=255)

       gTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACg                                                       >  1:1073808/1‑68 (MQ=255)

       gTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGc                                                      >  1:3063627/1‑69 (MQ=255)

       gTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGc                                                      >  1:2297097/1‑69 (MQ=255)

            gATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGcc                                                  <  1:2042809/68‑1 (MQ=255)

            gATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGCGGTTCCACGCTGcc                                                  <  1:3110154/68‑1 (MQ=255)

                tCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGc                                             >  1:2896411/1‑69 (MQ=255)

                     cAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTgg                                          >  1:2483732/1‑67 (MQ=255)

                     cAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTgg                                          >  1:1797328/1‑67 (MQ=255)

                       aaGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGaa                                      <  1:3430149/69‑1 (MQ=255)

                        aGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGaa                                      <  1:232505/68‑1 (MQ=255)

                          cgcGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGAt                                   <  1:2389255/67‑1 (MQ=255)

                          cccGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGAt                                   <  1:3295720/69‑1 (MQ=255)

                               aGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGc                              <  1:970243/69‑1 (MQ=255)

                               aGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGc                              <  1:2505430/69‑1 (MQ=255)

                               aGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGc                              <  1:254595/69‑1 (MQ=255)

                               aGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGc                              <  1:1500433/69‑1 (MQ=255)

                               aGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGc                              <  1:1360046/69‑1 (MQ=255)

                                   ggTAGAAAAATTTTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGaa                          <  1:2321171/69‑1 (MQ=255)

                                   ggTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGaa                          <  1:279427/69‑1 (MQ=255)

                                   ggTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGaa                          <  1:2680438/69‑1 (MQ=255)

                                   ggTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGaa                          <  1:1959626/69‑1 (MQ=255)

                                       gAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCt                                                     >  1:2677398/1‑38 (MQ=255)

                                       gAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCt                                                     >  1:3432952/1‑38 (MQ=255)

                                       gAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCt                                                     >  1:120395/1‑38 (MQ=255)

                                       gAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCt                                                     >  1:931818/1‑38 (MQ=255)

                                       gAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCt                                                     >  1:2673375/1‑38 (MQ=255)

                                       gAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCt                                                     >  1:2444475/1‑38 (MQ=255)

                                        aaaaaTATTGTCGGGCAGTAAAACGTGGTTCCACGCTGCCAATGCTGGCGAAGATg                                  >  1:2283444/1‑56 (MQ=255)

                                        aaaaaTATTGTCGGGCAGTAAAACGTGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAatattg                     >  1:1453337/1‑67 (MQ=255)

                                              atTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATCGGTGata               <  1:198873/69‑1 (MQ=255)

                                              atTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATCGGTGata               <  1:2462764/69‑1 (MQ=255)

                                               ttGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGAt                                   <  1:1538804/48‑1 (MQ=255)

                                               ttGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGAt                                   <  1:3261171/48‑1 (MQ=255)

                                               ttGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGAt                                   <  1:1015468/48‑1 (MQ=255)

                                               ttGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGAt                                   <  1:436073/48‑1 (MQ=255)

                                                           gAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATCGGTGATATGTTCGAAGTGgc  >  1:1985629/1‑69 (MQ=255)

                                                           gAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATCGGTGATATGTGCGAAGTgg   >  1:2948882/1‑68 (MQ=255)

                                                           gAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATCGGTGATATGTGCGAAGTGgc  >  1:833688/1‑69 (MQ=255)

                                                           gAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATCGGTGATATGTGCGAAGTGgc  >  1:3059468/1‑69 (MQ=255)

                                                           gAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATCGGTGATATGTGCGAAGTGgc  >  1:2815096/1‑69 (MQ=255)

                                                           gAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATCGGTGATATGTGCGAAGTGgc  >  1:1884103/1‑69 (MQ=255)

                                                           aaaaCGTGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATTGGTGATATGTGCGAAGTGgc  <  1:1222602/68‑1 (MQ=255)

                                                           aaaaCGTGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATTGGTGATATGTGCGAAGTGgc  <  1:357449/68‑1 (MQ=255)

                                                           aaaaCGTGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATTGGTGATATGTGCGAAGTGgc  <  1:852017/68‑1 (MQ=255)

                                                           aaaaCGTGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATTGGTGATATGTGCGAAGTGgc  <  1:1595207/68‑1 (MQ=255)

                                                           aaaaCGTGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATTGGTGATATGTGCGAAGTGgc  <  1:1987684/68‑1 (MQ=255)

                                                           aaaaCGTGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATTGGTGATATGTGCGAAGTGgc  <  1:1337201/68‑1 (MQ=255)

                                                           aaaaCGTGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATTGGTGATATGTGCGAAGTGgc  <  1:2149670/68‑1 (MQ=255)

                                                           aaaaCGTGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATTGGTGATATGTGCGAAGTGgc  <  1:2317069/68‑1 (MQ=255)

                                                            aaaCGGGGTTCCACGCTGCCAATGCTGGCGAAGAt                                   <  1:500721/35‑1 (MQ=255)

                                                            aaaCGGGGTTCCACGCTGCCAATGCTGGCGAAGAt                                   <  1:2233499/35‑1 (MQ=255)

                                                             aaCGTGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATTGGTGATATGTGCGAAGTGgc  <  1:3219153/67‑1 (MQ=255)

                                                               cGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGaa                          <  1:2514713/41‑1 (MQ=255)

                                                                 gggTTCCACGCTGCCAATGCTGGCGAAGATAACGCCGaa                          <  1:3016517/39‑1 (MQ=255)

GCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGTGAAACGGGGTTCCACGCTGCCAATGCTGGCGAAGATGACGCCGAATATCGGTGATATGTGCGAAGTGGC  >  minE/1374723‑1374850

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