BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	231,934	A→G	50.0%	Y276Y (TAT→TAC)	hemB ←	porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	231,934	0	A	G	50.0%	‑4.2 / 67.7	74	Y276Y (TAT→TAC)	hemB	porphobilinogen synthase
Reads supporting (aligned to +/- strand): ref base A (15/22); new base G (18/19); total (33/41)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

ACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACG  >  minE/231865‑231997

aCCTAAGCTTTCGATCACGACTTTCTCTTCATCTATAGCACCCGCC‑‑CGCGGCGAACTTAATCATCGCa                                                                 <  1:2378671/68‑1 (MQ=255)

 ccTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCa                                                                 <  1:2707567/69‑1 (MQ=255)

 ccTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCa                                                                 <  1:2944278/69‑1 (MQ=255)

  cTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCa                                                                 <  1:3213985/68‑1 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACt                                                             >  1:2678677/1‑68 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:3089066/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:2557499/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:232984/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:2605953/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:2810868/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:2011012/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:2962708/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:1766916/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:1746643/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:3073275/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:170226/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:1613103/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:1059586/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:3284562/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:989729/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:3399513/1‑69 (MQ=255)

      gCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGAGAACTTAATCATCGCGTACTc                                                            >  1:1294936/1‑69 (MQ=255)

                                cTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATc                                   <  1:1746146/68‑1 (MQ=255)

                                cTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATc                                   <  1:174716/68‑1 (MQ=255)

                                cTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATc                                   <  1:600800/68‑1 (MQ=255)

                                cTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATc                                   <  1:85943/68‑1 (MQ=255)

                                cTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATc                                   <  1:657648/68‑1 (MQ=255)

                                 tataGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCg                                  >  1:1031116/1‑68 (MQ=255)

                                 tataGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCg                                  >  1:1654798/1‑68 (MQ=255)

                                 tataGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCg                                  >  1:1178755/1‑68 (MQ=255)

                                 tataGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCg                                  >  1:728484/1‑68 (MQ=255)

                                 tataGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCg                                  >  1:615173/1‑68 (MQ=255)

                                 tataGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCg                                  >  1:14640/1‑68 (MQ=255)

                                 tataGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCg                                  >  1:3073480/1‑68 (MQ=255)

                                                 gcggcgAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCGAAACGGCAATTCAGTACGTTCa                 >  1:1037594/1‑69 (MQ=255)

                                                 gcggcgAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCa                 >  1:204957/1‑69 (MQ=255)

                                                 gcggcgAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCa                 >  1:2848239/1‑69 (MQ=255)

                                                 gcggcgAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCa                 >  1:1609540/1‑69 (MQ=255)

                                                 gcggcgAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCa                 >  1:1210044/1‑69 (MQ=255)

                                                 gcggcgAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCa                 >  1:2817625/1‑69 (MQ=255)

                                                 gcggcgAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCa                 >  1:1402673/1‑69 (MQ=255)

                                                 gcggcgAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATAGGCAATTCAGTACGTTCa                 >  1:1460176/1‑69 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:1303507/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:2410674/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:1692999/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:575873/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:1727324/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:666940/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:702179/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:808929/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:814016/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:1098597/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:887666/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:1065422/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:2009048/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:2444968/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:2839115/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:2340981/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:2200611/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:2643661/43‑1 (MQ=255)

                                                       aaCTTAATCATCGCGTACTCACCGCTCACCTGATACGCGCCaa                                     <  1:2264512/43‑1 (MQ=255)

                                                        aCTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATc                                   <  1:2444876/44‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:2420898/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:2963719/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:871843/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:2520857/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:2308722/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:2654237/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:1274193/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:641480/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:2705012/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:592864/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:2054190/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:1491740/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:3126402/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:2858761/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:308488/69‑1 (MQ=255)

                                                                atcGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg  <  1:1729473/69‑1 (MQ=255)

ACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACG  >  minE/231865‑231997

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