BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,665,986	T→A	61.9%	V257V (GTA→GTT)	guaB ←	IMP dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,665,986	0	T	A	61.9%	24.7 / 35.8	42	V257V (GTA→GTT)	guaB	IMP dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (7/9); new base A (8/18); total (15/27)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

GCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTC  >  minE/1665921‑1666054

gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                                     <  1:2510776/67‑1 (MQ=255)

gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                                     <  1:2929705/67‑1 (MQ=255)

gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                                     <  1:840657/67‑1 (MQ=255)

gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                                     <  1:2545189/67‑1 (MQ=255)

gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                                     <  1:620743/67‑1 (MQ=255)

gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                                     <  1:1732689/67‑1 (MQ=255)

                 gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                     >  1:2878472/2‑66 (MQ=255)

                 gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                     >  1:3211438/2‑66 (MQ=255)

                 gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                     >  1:1834286/2‑66 (MQ=255)

                 gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                     >  1:485990/2‑66 (MQ=255)

                 gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                     >  1:1704295/2‑66 (MQ=255)

                 gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                     >  1:1356041/2‑66 (MQ=255)

                 gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                     >  1:858483/2‑66 (MQ=255)

                 gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                     >  1:997326/2‑66 (MQ=255)

                    ttGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTggg                                                <  1:2114374/68‑1 (MQ=255)

                            ccGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCg                                        >  1:2786754/1‑68 (MQ=255)

                                   tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcaggag                                >  1:2963779/1‑69 (MQ=255)

                                   tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                                >  1:378671/1‑69 (MQ=255)

                                   tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                                >  1:2591836/1‑69 (MQ=255)

                                   tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                                >  1:367638/1‑69 (MQ=255)

                                   tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                                >  1:612765/1‑69 (MQ=255)

                                   tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                                >  1:1066890/1‑69 (MQ=255)

                                         aCGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTc                          <  1:3001855/69‑1 (MQ=255)

                                          cGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTc                          <  1:1158920/68‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTTCGGCAAc           <  1:2369249/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:3092124/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:1248622/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:941699/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:1350574/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:1648556/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:1813324/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:1816757/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:1897878/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:2058640/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:2635992/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:3071324/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:3068670/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:3014893/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:3001844/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:2529081/69‑1 (MQ=255)

                                                        aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:2606107/69‑1 (MQ=255)

                                                              cAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc           <  1:3076647/63‑1 (MQ=255)

                                                                 tACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTc  <  1:1871025/69‑1 (MQ=255)

                                                                 tACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTc  <  1:2299277/69‑1 (MQ=255)

                                                                 tACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTc  <  1:1276177/69‑1 (MQ=255)

GCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTC  >  minE/1665921‑1666054

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