BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	121,082	A→G	21.0%	A41A (GCT→GCC)	speE ←	spermidine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	121,082	0	A	G	21.0%	143.6 / 29.1	104	A41A (GCT→GCC)	speE	spermidine synthase
Reads supporting (aligned to +/- strand): ref base A (58/24); new base G (14/8); total (72/32)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

ATAGATAAACTCGTCGCGCTCGGTGGTTTGTACTACGCCATCCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTT  >  minE/121016‑121138

                                 tacGCCATCCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGt                       >  1:1721069/1‑69 (MQ=255)

                                 tacGCCATCCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGt                       >  1:2702233/1‑69 (MQ=255)

                                 tacGCCATCCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGt                       >  1:2705987/1‑69 (MQ=255)

                                      cATCCAGCGCCATTACGGGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCtttt                  <  1:2830266/69‑1 (MQ=255)

                                      cATCCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCtttt                  <  1:220534/69‑1 (MQ=255)

                                      cATCCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCtttt                  <  1:67202/69‑1 (MQ=255)

                                      cATCCAGCGCCAATACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCtttt                  <  1:3402965/69‑1 (MQ=255)

                                          cAGCGCCATTACGCGACCAAATGCAGCGTTCTCaa                                                >  1:882218/1‑35 (MQ=255)

                                          cAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCAt               >  1:991407/1‑68 (MQ=255)

                                          cAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCAt               >  1:1986174/1‑68 (MQ=255)

                                          cAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCAt               >  1:476477/1‑68 (MQ=255)

                                          cAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCAt               >  1:607714/1‑68 (MQ=255)

                                              gcCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATAc          <  1:2492269/69‑1 (MQ=255)

                                              gcCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATAc          <  1:1346236/69‑1 (MQ=255)

                                                cATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAg        <  1:822201/69‑1 (MQ=255)

                                                cATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAg        <  1:953524/69‑1 (MQ=255)

                                                cATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAg        <  1:1070557/69‑1 (MQ=255)

                                                cATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAg        <  1:1124668/69‑1 (MQ=255)

                                                cATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAg        <  1:1159100/69‑1 (MQ=255)

                                                cATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAg        <  1:1315894/69‑1 (MQ=255)

                                                cATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGCTGATCGGTCTTTTCATGATACAg        <  1:1182570/69‑1 (MQ=255)

                                                    aggcgACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCg                         >  1:1572872/3‑48 (MQ=255)

                                                    aggcgACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCg                         >  1:1577270/3‑48 (MQ=255)

                                                    aCGCGACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACg    <  1:1392477/69‑1 (MQ=255)

                                                    aCGCGACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACg    <  1:3409426/69‑1 (MQ=255)

                                                    aCGCGACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACg    <  1:3296256/69‑1 (MQ=255)

                                                    aCGCGACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACg    <  1:2361023/69‑1 (MQ=255)

                                                    aCGCGACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACg    <  1:2326154/69‑1 (MQ=255)

                                                    aCGCGACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACg    <  1:1327693/69‑1 (MQ=255)

                                                    aCGCGACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACg    <  1:209059/69‑1 (MQ=255)

                                                    aCGCGACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACg    <  1:1871306/69‑1 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAg                                      >  1:1338240/1‑35 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCg                         >  1:3033204/1‑48 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCg                         >  1:2496731/1‑48 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCg                         >  1:3129567/1‑48 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCg                         >  1:819971/1‑48 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCg                         >  1:2878814/1‑48 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCg                         >  1:2466128/1‑48 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCg                         >  1:230180/1‑48 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCg                         >  1:2737432/1‑48 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGaa      >  1:2966907/1‑67 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:351218/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:3092194/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:1463287/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:3351960/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:2703197/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:1303040/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:2730708/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:1600480/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:2819572/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:1223163/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:2903847/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:3135221/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:3007330/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:3079852/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:382050/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:1849901/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:1085711/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:2129245/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:2154676/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:1113892/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:219664/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:1343192/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:77351/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:230405/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:1395897/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:231612/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:1719068/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:1713582/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:530638/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:2473965/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:382235/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCACGATACAGAAc     >  1:3168541/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTACATGATACAGAAc     >  1:2521095/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTCTCAAAAATGATCAGAACCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:869856/1‑68 (MQ=255)

                                                    aCGCGACCAAATGCAGCGTTATCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAAc     >  1:1946455/1‑68 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGTTGATCGGTCTTTTc                 >  1:2941366/1‑53 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCg                         >  1:2511320/1‑45 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTc                 >  1:3003383/1‑53 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTc                 >  1:104104/1‑53 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTc                 >  1:1998020/1‑53 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTc                 >  1:1109991/1‑53 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTc                 >  1:869578/1‑53 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTc                 >  1:2157025/1‑53 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTc                 >  1:116826/1‑53 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTc                 >  1:2533068/1‑53 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTc                 >  1:1675646/1‑53 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTc                 >  1:2788399/1‑53 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTc                 >  1:2798289/1‑53 (MQ=255)

                                                       cgACCAAATGCGGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTc                 >  1:3122133/1‑53 (MQ=255)

                                                       cgACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGtt  <  1:1316126/68‑1 (MQ=255)

                                                       cgACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGtt  <  1:1416139/68‑1 (MQ=255)

                                                       cgACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGtt  <  1:2403724/68‑1 (MQ=255)

                                                       cgACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGtt  <  1:2307958/68‑1 (MQ=255)

                                                       cgACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGtt  <  1:1347123/68‑1 (MQ=255)

                                                       cgACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGtt  <  1:2248441/68‑1 (MQ=255)

                                                       cgACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGtt  <  1:2161947/68‑1 (MQ=255)

                                                       cgACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGtt  <  1:2566218/68‑1 (MQ=255)

                                                       cgACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGtt  <  1:862632/68‑1 (MQ=255)

                                                       cgACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGtt  <  1:1721116/68‑1 (MQ=255)

ATAGATAAACTCGTCGCGCTCGGTGGTTTGTACTACGCCATCCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTT  >  minE/121016‑121138

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