BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	494,622	C→T	41.4%	S224N (AGC→AAC)	galM ←	galactose‑1‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	494,622	0	C	T	41.4%	13.5 / 25.8	41	S224N (AGC→AAC)	galM	galactose‑1‑epimerase
Reads supporting (aligned to +/- strand): ref base C (8/16); new base T (3/14); total (11/30)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAAT  >  minE/494560‑494674

cGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGa                                                <  1:3179997/69‑1 (MQ=255)

cGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGa                                                <  1:1109572/69‑1 (MQ=255)

cGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGa                                                <  1:2330096/69‑1 (MQ=255)

cGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGa                                                <  1:1492405/69‑1 (MQ=255)

cGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGa                                                <  1:1637358/69‑1 (MQ=255)

cGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGa                                                <  1:1695321/69‑1 (MQ=255)

cGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGa                                                <  1:1842817/69‑1 (MQ=255)

 gTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGa                                                <  1:361939/68‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAATAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATCAAAAGAGGtt                                  <  1:2478732/69‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAAGAAACTTACTGGCGATGATTTTGGCGTTGCGGAAATCAAAAGAGGtt                                  <  1:515897/69‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATCATAAGAGGtt                                  <  1:3187992/69‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATCAAAAGAGGtt                                  <  1:1156817/69‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATCAAAAGAGGtt                                  <  1:2624177/69‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATCAAAAGAGGtt                                  <  1:21268/69‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATCAAAAGAGGtt                                  <  1:204821/69‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATCAAAAGAGGtt                                  <  1:2028164/69‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATCAAAAGAGGtt                                  <  1:1870651/69‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATCAAAAGAGGtt                                  <  1:1789988/69‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATCAAAAGAGGtt                                  <  1:3411504/69‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATCAAAAGAGGtt                                  <  1:1388792/69‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATCAAAAGAGGtt                                  <  1:501275/69‑1 (MQ=255)

              ttGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATCAAAAGAGGtt                                  <  1:1315619/69‑1 (MQ=255)

                        cgtcgtcgGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAg                        <  1:2560651/69‑1 (MQ=255)

                        cgtcgtcgGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAg                        <  1:3159395/69‑1 (MQ=255)

                        cgtcgtcgGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAg                        <  1:3155382/69‑1 (MQ=255)

                        cgtcgtcgGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAg                        <  1:2899412/69‑1 (MQ=255)

                        cgtcgtcgGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAg                        >  1:2884546/1‑69 (MQ=255)

                        cgtcgtcgGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAg                        <  1:2664941/69‑1 (MQ=255)

                        cgtcgtcgGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAg                        <  1:2589992/69‑1 (MQ=255)

                         gtcgtcgGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAg                        <  1:1467962/68‑1 (MQ=255)

                                    gAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATc                                            >  1:2933573/1‑37 (MQ=255)

                                    gAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATc                                            >  1:476784/1‑37 (MQ=255)

                                    gAAACTCACTGGCGATGATTTTGGCGTTGCGGAAATc                                            >  1:1333033/1‑37 (MQ=255)

                                      aaCTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCcgtcg          <  1:876945/69‑1 (MQ=255)

                                              ggCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGaa   >  1:700229/1‑68 (MQ=255)

                                              ggCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAAt  >  1:3107556/1‑69 (MQ=255)

                                              ggCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAAt  >  1:2758481/1‑69 (MQ=255)

                                              ggCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAAt  >  1:1865409/1‑69 (MQ=255)

                                              ggCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAAt  >  1:1450168/1‑69 (MQ=255)

                                              ggCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAAt  >  1:1244484/1‑69 (MQ=255)

                                              ggCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAAt  >  1:680337/1‑69 (MQ=255)

CGTAACCTTTCACTTTGCGCTGATCGTCGTCGGCAAGAAACTCACTGGCGATGATTTTGGCGCTGCGGAAATCAAAAGACGTTCCGGCGACAGATTTCAGGCCGTCGTGCGGAAT  >  minE/494560‑494674

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