BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,010,698	A→G	50.0%	Q14Q (CAA→CAG)	yghZ →	aldo‑keto reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,010,698	0	A	G	50.0%	‑4.2 / 59.5	66	Q14Q (CAA→CAG)	yghZ	aldo‑keto reductase
Reads supporting (aligned to +/- strand): ref base A (19/14); new base G (16/17); total (35/31)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAATTGAACGGTGGAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGTTTATGGCACAATTT  >  minE/2010630‑2010763

gAAATTGAACGGTGGAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCaa                                                                   >  1:831801/1‑69 (MQ=255)

gAAATTGAACGGTGGAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCaa                                                                   >  1:1944214/1‑69 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTCCGGGCAGATGCAATACCGCTATTGCg                                                      >  1:1457691/1‑68 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGCTGCAGTACCGCTATTGCg                                                      >  1:351644/1‑68 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGc                                                       >  1:902925/1‑67 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGc                                                       >  1:1429459/1‑67 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGc                                                       >  1:159445/1‑67 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCg                                                      >  1:976645/1‑68 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCg                                                      >  1:320521/1‑68 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCg                                                      >  1:3302132/1‑68 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCg                                                      >  1:3387239/1‑68 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCg                                                      >  1:2476489/1‑68 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCg                                                      >  1:475187/1‑68 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCg                                                      >  1:509259/1‑68 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCg                                                      >  1:2162959/1‑68 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCg                                                      >  1:749345/1‑68 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCg                                                      >  1:860029/1‑68 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCg                                                      >  1:1840344/1‑68 (MQ=255)

              gAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCg                                                      >  1:1671457/1‑68 (MQ=255)

                         tCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGtt                                          <  1:2991454/69‑1 (MQ=255)

                             ggTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGTAAAAAGTGGTTTAc                                       <  1:1355398/68‑1 (MQ=255)

                             ggTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:517526/68‑1 (MQ=255)

                             ggTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:2910626/68‑1 (MQ=255)

                             ggTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:459612/68‑1 (MQ=255)

                             ggTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:1536766/68‑1 (MQ=255)

                             ggTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:1762838/68‑1 (MQ=255)

                             ggTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:2193360/68‑1 (MQ=255)

                             ggTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:1796746/68‑1 (MQ=255)

                             ggTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:680592/68‑1 (MQ=255)

                             ggTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:79302/68‑1 (MQ=255)

                             ggTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:2927487/68‑1 (MQ=255)

                             ggTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:808102/68‑1 (MQ=255)

                             ggTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGATATTGCGGAAAAAGTGGTTTAc                                       <  1:1861879/68‑1 (MQ=255)

                             ggTCTGGTTACCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:2365488/68‑1 (MQ=255)

                                 tGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTg                                  <  1:3369482/69‑1 (MQ=255)

                                 tGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTg                                  <  1:2046376/69‑1 (MQ=255)

                                  ggTTAGCGAATCCCGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:1420417/63‑1 (MQ=255)

                                              ccGAACGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:2081702/51‑1 (MQ=255)

                                              ccGAACGTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCg                      >  1:3383971/1‑68 (MQ=255)

                                              ccGAACGTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGc                     >  1:799046/1‑69 (MQ=255)

                                              ccGAACGTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGc                     >  1:1249503/1‑69 (MQ=255)

                                              ccGAACGTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGc                     >  1:208448/1‑69 (MQ=255)

                                              ccGAACGTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGc                     >  1:2548728/1‑69 (MQ=255)

                                              ccGAACGTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGc                     >  1:2121166/1‑69 (MQ=255)

                                                   cGTTACGGGCAGATGCAGTACCGCTATTGCGGAAAAAGTGGTTTAc                                       <  1:2871112/46‑1 (MQ=255)

                                                   cGTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTg                                  <  1:420975/51‑1 (MQ=255)

                                                    gTTACGTGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGtt               <  1:2338408/69‑1 (MQ=255)

                                                    gTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTg                                  <  1:2788113/50‑1 (MQ=255)

                                                    gTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGtt               <  1:2849089/69‑1 (MQ=255)

                                                    gTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGtt               <  1:2113484/69‑1 (MQ=255)

                                                    gTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGtt               <  1:894723/69‑1 (MQ=255)

                                                    gTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGtt               <  1:2514985/69‑1 (MQ=255)

                                                    gTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGtt               <  1:77819/69‑1 (MQ=255)

                                                    gTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGtt               <  1:49488/69‑1 (MQ=255)

                                                    gTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGtt               <  1:2189200/69‑1 (MQ=255)

                                                    gTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGtt               <  1:2152558/69‑1 (MQ=255)

                                                                tGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGTTTATGGCACAAtt   >  1:143604/1‑69 (MQ=255)

                                                                tGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGTTTATGGCACAAtt   >  1:856273/1‑69 (MQ=255)

                                                                tGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGTTTATGGCACAAtt   >  1:1050403/1‑69 (MQ=255)

                                                                tGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGTTTATGGCACAAtt   >  1:1425673/1‑69 (MQ=255)

                                                                tGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGTTTATGGCACAAtt   >  1:3117582/1‑69 (MQ=255)

                                                                tGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGTTTATGGCACAAtt   >  1:1847085/1‑69 (MQ=255)

                                                                tGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGTTTATGGCACAAtt   >  1:701578/1‑69 (MQ=255)

                                                                tGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGTTTATGGCACAAtt   >  1:2419589/1‑69 (MQ=255)

                                                                tGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGTTTATGGCACAAtt   >  1:2916231/1‑69 (MQ=255)

                                                                tGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGTTTATGGCACAAt    >  1:1477787/1‑68 (MQ=255)

                                                                tGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGACCGCGTTATCGCTCGGTTTATGGCACAAtt   >  1:447172/1‑69 (MQ=255)

                                                                tGCAAT‑CCCCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGTTTATGGCACAAttt  >  1:440652/1‑69 (MQ=255)

GAAATTGAACGGTGGAAAGGAGAGGTCATGGTCTGGTTAGCGAATCCCGAACGTTACGGGCAGATGCAATACCGCTATTGCGGAAAAAGTGGTTTACGCCTGCCCGCGTTATCGCTCGGTTTATGGCACAATTT  >  minE/2010630‑2010763

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