BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,529,224	A→G	37.9%	G66G (GGT→GGC)	yfcL ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,529,224	0	A	G	37.9%	24.2 / 29.2	50	G66G (GGT→GGC)	yfcL	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (12/19); new base G (12/7); total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGCTGAAACAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGCCTGAGCGGAGTG  >  minE/1529157‑1529280

ccTGCTGAAACAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCAc                                                         <  1:2449057/69‑1 (MQ=255)

 ctgctgAAACAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCAcc                                                        <  1:390189/69‑1 (MQ=255)

         aCAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCtt                                                >  1:2991741/1‑69 (MQ=255)

         aCAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCtt                                                >  1:2307704/1‑69 (MQ=255)

         aCAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCtt                                                >  1:1816981/1‑69 (MQ=255)

         aCAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCtt                                                >  1:1724266/1‑69 (MQ=255)

         aCAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACAAAt                                                     >  1:2147162/1‑64 (MQ=255)

         aCAGGTGTTCCCACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCtt                                                >  1:1128/1‑69 (MQ=255)

               tttCCCACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCa                                          >  1:2414813/1‑69 (MQ=255)

                   ccACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGAc                                       <  1:1685874/68‑1 (MQ=255)

                   ccACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGAc                                       <  1:873819/68‑1 (MQ=255)

                   ccACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGAc                                       <  1:1054599/68‑1 (MQ=255)

                   ccACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGAc                                       <  1:1895979/68‑1 (MQ=255)

                       atatCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTg                                   >  1:375978/1‑68 (MQ=255)

                       atatCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTg                                   >  1:693037/1‑68 (MQ=255)

                       atatCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTg                                   >  1:243717/1‑68 (MQ=255)

                        tatCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCAcc                                                        <  1:966844/46‑1 (MQ=255)

                        tatCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCAcc                                                        <  1:9073/46‑1 (MQ=255)

                        tatCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCAcc                                                        <  1:3323239/46‑1 (MQ=255)

                        tatCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCAcc                                                        <  1:122858/46‑1 (MQ=255)

                        tatCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCAcc                                                        <  1:2153958/46‑1 (MQ=255)

                        tatCGGTCACCAGCCCCTGGTCACGCGGCGACAATTCACCAGCAcc                                                        <  1:1476116/46‑1 (MQ=255)

                        tatCGGTCACCAGCACCTGGTCACGCGGCGACAATTCACCAGCAcc                                                        <  1:3359985/46‑1 (MQ=255)

                                   aGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTa                      >  1:1717842/1‑69 (MQ=255)

                                    gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATgg                                                   >  1:3022572/1‑39 (MQ=37)

                                    gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                      >  1:1898416/1‑68 (MQ=255)

                                    gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                      >  1:2197430/1‑68 (MQ=255)

                                    gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                      >  1:1987156/1‑68 (MQ=255)

                                    gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                      >  1:2040647/1‑68 (MQ=255)

                                    gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                      >  1:3301481/1‑68 (MQ=255)

                                    gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                      >  1:1141408/1‑68 (MQ=255)

                                    gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                      >  1:2831688/1‑68 (MQ=255)

                                    gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                      >  1:20817/1‑68 (MQ=255)

                                    gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                      >  1:2118921/1‑68 (MQ=255)

                                    gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                      >  1:1744965/1‑68 (MQ=255)

                                    gcgcCTGGTCACGCGGTGACAATTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTa                      >  1:1619735/1‑68 (MQ=255)

                                     cgcCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATg                    <  1:2778017/69‑1 (MQ=255)

                                     cgcCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATg                    <  1:3346363/69‑1 (MQ=255)

                                      gcCTGGTCACGCGGCGACAATTCACCAGCACCATTGGCTTTTTCCAGAGTCTGGCTAACGGTCGTATg                    <  1:618649/68‑1 (MQ=255)

                                         tGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATGCAcc                <  1:534212/69‑1 (MQ=255)

                                          ggTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATGCAcc                <  1:2724270/68‑1 (MQ=255)

                                             cACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGc              >  1:2519614/1‑67 (MQ=255)

                                                   gcgACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGCCTGAGCgg      <  1:874159/69‑1 (MQ=255)

                                                        aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg  <  1:624053/68‑1 (MQ=255)

                                                        aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg  <  1:483136/68‑1 (MQ=255)

                                                        aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg  <  1:821134/68‑1 (MQ=255)

                                                        aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg  <  1:435577/68‑1 (MQ=255)

                                                        aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg  <  1:3431122/68‑1 (MQ=255)

                                                        aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg  <  1:1179411/68‑1 (MQ=255)

                                                        aaTTCACCCGCGCCAATGGCTTTTTCCAGACTCTGGCTAACGGTTGTATGTACCGCCTGAGCGGAGTg  <  1:977265/68‑1 (MQ=255)

CCTGCTGAAACAGGTTTTCCCACATATCGGTCACCAGCGCCTGGTCACGCGGCGACAATTCACCAGCACCAATGGCTTTTTCCAGACTCTGGCTAACGGTCGTATGCACCGCCTGAGCGGAGTG  >  minE/1529157‑1529280

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