BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	327,396	C→G	53.4%	A291A (GCC→GCG)	hemH →	ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	327,396	0	C	G	53.4%	1.0 / 30.8	30	A291A (GCC→GCG)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand): ref base C (8/6); new base G (10/6); total (18/12)
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

GATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTA  >  minE/327337‑327463

gATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGaa                                                             >  1:1864116/1‑68 (MQ=255)

gATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGa                                                              >  1:1219268/1‑67 (MQ=255)

  ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                                           >  1:1337251/1‑68 (MQ=255)

  ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                                           >  1:1555911/1‑68 (MQ=255)

  ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                                           >  1:1569480/1‑68 (MQ=255)

  ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                                           >  1:462241/1‑68 (MQ=255)

  ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                                           >  1:226273/1‑68 (MQ=255)

  ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                                           >  1:3397521/1‑68 (MQ=255)

  ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                                           >  1:3389383/1‑68 (MQ=255)

  ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                                           >  1:3340371/1‑68 (MQ=255)

  ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                                           >  1:3029049/1‑68 (MQ=255)

  ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGaaaa                                                           >  1:2795573/1‑68 (MQ=255)

            aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGAAAAAATATGAatat                                                <  1:908913/69‑1 (MQ=255)

            aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGAAAAAATATGAatat                                                <  1:1003171/69‑1 (MQ=255)

            aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGAAAAAATATGAatat                                                <  1:524140/69‑1 (MQ=255)

            aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGAAAAAATATGAatat                                                <  1:2888558/69‑1 (MQ=255)

            aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGAAAAAATATGAatat                                                <  1:2769788/69‑1 (MQ=255)

            aCGCTGCAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCGGGCGGGAAAAAATATGAatat                                                <  1:1678927/69‑1 (MQ=255)

                     gagaTTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgc                                        >  1:80764/1‑68 (MQ=255)

                       gaTTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgcgc                                      >  1:3313461/1‑68 (MQ=255)

                              gAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGc                               <  1:2415645/68‑1 (MQ=255)

                                     aCCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCgg                       <  1:3088028/69‑1 (MQ=255)

                                     aCCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCgg                       <  1:238497/69‑1 (MQ=255)

                                               ctcCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgcg                                       <  1:820704/43‑1 (MQ=255)

                                               ctcCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgcg                                       <  1:1370512/43‑1 (MQ=255)

                                                tcCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgcg                                       <  1:2512692/42‑1 (MQ=255)

                                                    tCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAatga        >  1:306521/1‑69 (MQ=255)

                                                    tCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAatga        >  1:3093422/1‑69 (MQ=255)

                                                    tCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAatga        >  1:607756/1‑69 (MQ=255)

                                                    tCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAatg         >  1:698713/1‑68 (MQ=255)

                                                          ccGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTa  >  1:2184820/1‑69 (MQ=255)

GATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTCCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTA  >  minE/327337‑327463

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