BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,966,357	A→G	50.0%	Q155Q (CAA→CAG)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,966,357	0	A	G	50.0%	‑4.7 / 64.1	70	Q155Q (CAA→CAG)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand): ref base A (10/25); new base G (9/26); total (19/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAG  >  minE/2966295‑2966423

cgcgCTAATCTGCCGGGTTACCTCGGTAACTTCCACTCTTCCGGTACGGTCATTCTGGATCAACTggg                                                               >  1:278124/1‑68 (MQ=255)

cgcgCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTggg                                                               >  1:3387352/1‑68 (MQ=255)

         cTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGcac                                                     <  1:489703/69‑1 (MQ=255)

         cTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGcac                                                     <  1:1029239/69‑1 (MQ=255)

         cTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGcac                                                     <  1:2472517/69‑1 (MQ=255)

         cTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGcac                                                     <  1:2997941/69‑1 (MQ=255)

         cTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGcac                                                     <  1:219796/69‑1 (MQ=255)

         cTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGcac                                                     <  1:212006/69‑1 (MQ=255)

         cTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGcac                                                     <  1:2048196/69‑1 (MQ=255)

         cTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGcac                                                     <  1:1875029/69‑1 (MQ=255)

         cTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGcac                                                     <  1:885410/69‑1 (MQ=255)

         cTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGcac                                                     <  1:1611769/69‑1 (MQ=255)

         cTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGACcac                                                     <  1:2657961/69‑1 (MQ=255)

         cTGCCGGGTTACCTCGGTAACTGCCACTCTTCAGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGcac                                                     <  1:557249/69‑1 (MQ=255)

              gggTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGaa                                                <  1:425637/69‑1 (MQ=255)

                  tACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAAcc                                            >  1:2629946/1‑69 (MQ=255)

                  tACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAAcc                                            >  1:2610107/1‑69 (MQ=255)

                  tACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAAcc                                            >  1:3260599/1‑69 (MQ=255)

                   tccTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:506843/67‑1 (MQ=255)

                       cGGTAACTGCCACTCTTCCGGTACGGTCATTCTTGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:2277161/64‑1 (MQ=255)

                       cGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:144346/64‑1 (MQ=255)

                       cGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:1645590/64‑1 (MQ=255)

                       cGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:2738312/64‑1 (MQ=255)

                       cGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:1880282/64‑1 (MQ=255)

                       cGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:2641628/64‑1 (MQ=255)

                       cGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:801743/64‑1 (MQ=255)

                       cGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:3274833/64‑1 (MQ=255)

                       cGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:3118874/64‑1 (MQ=255)

                       cGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:2465464/64‑1 (MQ=255)

                       cGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:2876225/64‑1 (MQ=255)

                       cGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:2107529/64‑1 (MQ=255)

                       cGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAAcc                                            <  1:352041/64‑1 (MQ=255)

                            aCTGCCACTCTTCCGGTACGGTCATTCTTGATCAACTGGGCGAAGAGCACATg                                                  <  1:2551116/53‑1 (MQ=255)

                            aCTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATg                                                  <  1:679524/53‑1 (MQ=255)

                            aCTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATg                                                  <  1:800000/53‑1 (MQ=255)

                             cTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATg                                                  <  1:1675050/52‑1 (MQ=255)

                                  aCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:771203/68‑1 (MQ=255)

                                  aCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:2953964/68‑1 (MQ=255)

                                  aCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:2031229/68‑1 (MQ=255)

                                  aCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:2390338/68‑1 (MQ=255)

                                  aCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:2198260/68‑1 (MQ=255)

                                  aCTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:2780803/68‑1 (MQ=255)

                                   ctctTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCta                           <  1:593155/69‑1 (MQ=255)

                                           gTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCTAGCCGATTTTCTATACCTCCg                   >  1:2964020/1‑69 (MQ=255)

                                           gTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCg                   >  1:476502/1‑69 (MQ=255)

                                           gTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCg                   >  1:95858/1‑69 (MQ=255)

                                           gTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCg                   >  1:582810/1‑69 (MQ=255)

                                               ggTCATTCTTGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:360473/55‑1 (MQ=255)

                                                  cATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTc            >  1:215792/1‑69 (MQ=255)

                                                  cATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTc            >  1:574502/1‑69 (MQ=255)

                                                  cATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTc            >  1:2085281/1‑69 (MQ=255)

                                                  cATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTc            >  1:2441035/1‑69 (MQ=255)

                                                  cATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTc            >  1:1635294/1‑69 (MQ=255)

                                                  cATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTc            >  1:1515016/1‑69 (MQ=255)

                                                  cATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTc            >  1:1468433/1‑69 (MQ=255)

                                                  cATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTc            >  1:416811/1‑69 (MQ=255)

                                                  cATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTc            >  1:97069/1‑69 (MQ=255)

                                                   aTTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:214835/51‑1 (MQ=255)

                                                            cAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:2774351/42‑1 (MQ=255)

                                                            cAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:3349855/42‑1 (MQ=255)

                                                            cAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:3248342/42‑1 (MQ=255)

                                                            cAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:2849132/42‑1 (MQ=255)

                                                            cAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:2533145/42‑1 (MQ=255)

                                                            cAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:2425710/42‑1 (MQ=255)

                                                            cAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:1087065/42‑1 (MQ=255)

                                                            cAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:1076594/42‑1 (MQ=255)

                                                            cAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAg  <  1:3021173/69‑1 (MQ=255)

                                                            cAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAg  >  1:2966081/1‑69 (MQ=255)

                                                            cAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAg  <  1:1937080/69‑1 (MQ=255)

                                                            cAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAg  <  1:1460052/69‑1 (MQ=255)

                                                             aaCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                             <  1:3381773/41‑1 (MQ=255)

                                                              aCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAg  <  1:617994/67‑1 (MQ=255)

CGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAG  >  minE/2966295‑2966423

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