BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,697,162	T→C	37.2%	V38V (GTT→GTC)	yfhR →	predicted peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,697,162	0	T	C	37.2%	20.5 / 24.5	43	V38V (GTT→GTC)	yfhR	predicted peptidase
Reads supporting (aligned to +/- strand): ref base T (16/11); new base C (8/8); total (24/19)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

AAAATTCTGTTTATTCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATAAAATTTACGGTCCCGATC  >  minE/1697106‑1697211

aaaaTTCTGTTTATTCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTCGCCATCaa                                           <  1:1127042/65‑1 (MQ=255)

aaaaTTCTGTTTATTCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTCGCCATCaa                                           <  1:1189062/65‑1 (MQ=255)

aaaaTTCTGTTTATTCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTCGCCATCaa                                           <  1:774313/65‑1 (MQ=255)

aaaaTTCTGTTTATTCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTCGCCATCaa                                           <  1:2837147/65‑1 (MQ=255)

aaaaTTCTGTTTATTCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTCGCCATCaa                                           <  1:222046/65‑1 (MQ=255)

aaaaTTCTGTTTATTCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTCGCCATCaa                                           <  1:218374/65‑1 (MQ=255)

aaaaTTCTGTTTATTCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTCGCCATCaa                                           <  1:2144557/65‑1 (MQ=255)

aaaaTTCTGTTTATTCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTCGCCATCaa                                           <  1:1972123/65‑1 (MQ=255)

      cTGTTTATTCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCtat                                 <  1:983174/69‑1 (MQ=255)

      cTGTTTATTCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCtat                                 <  1:2363072/69‑1 (MQ=255)

      cTGTTTATTCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCtat                                 <  1:1951036/69‑1 (MQ=255)

      cTGTTTATTCTCTTAGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCtat                                 <  1:1912989/69‑1 (MQ=255)

             ttCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATccc                           <  1:54660/68‑1 (MQ=255)

                    tgttgttgCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGAt                     >  1:3118918/1‑67 (MQ=255)

                    tgttgttgCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATaa                   >  1:1888481/1‑69 (MQ=255)

                    tgttgttgCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATaa                   >  1:3036102/1‑69 (MQ=255)

                    tgttgttgCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATaa                   >  1:732831/1‑69 (MQ=255)

                    tgttgttgCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATaa                   >  1:1759636/1‑69 (MQ=255)

                    tgttgttgCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATaa                   >  1:2108990/1‑69 (MQ=255)

                    tgttgttgCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATaa                   >  1:943531/1‑69 (MQ=255)

                        gttgCCTTCTGCGTTTATTTAGTGCCGGGCGTTGCCATCAACTTCTTCTATTATCCCGACGATAAAAtt               >  1:411914/1‑69 (MQ=255)

                        gttgCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATAAAAtt               >  1:755655/1‑69 (MQ=255)

                        gttgCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATAAAAtt               >  1:120525/1‑69 (MQ=255)

                        gttgCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATAAAAtt               >  1:224516/1‑69 (MQ=255)

                        gttgCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATAAAAtt               >  1:3126520/1‑69 (MQ=255)

                                 tGCGTTTATTTAGTGCCGCGCGTCGCCATCAACTTCTTCTATTATCCTGACGATAAAATATACGGTccc      >  1:886833/1‑69 (MQ=255)

                                 tGCGTTTATTTAGTGCCGCGCGTCGCCATCAACTTCTTCTATTATCCTGACGATAAAATATACGGTccc      >  1:1109659/1‑69 (MQ=255)

                                 tGCGTTTATTTAGTGCCGCGCGTCGCCATCAACTTCTTCTATTATCCTGACGATAAAATATACGGTccc      >  1:2790804/1‑69 (MQ=255)

                                 tGCGTTTATTTAGTGCCGCGCGTCGCCATCAACTTCTTCTATTATCCTGACGATAAAATATACGGTccc      >  1:265987/1‑69 (MQ=255)

                                 tGCGTTTATTTAGTGCCGCGCGTCGCCATCAACTTCTTCTATTATCCTGACGATAAAATATACGGTccc      >  1:254073/1‑69 (MQ=255)

                                 tGCGTTTATTTAGTGCCGCGCGTCGCCATCAACTTCTTCTATTATCCTGACGATAAAATATACGGTccc      >  1:2378316/1‑69 (MQ=255)

                                 tGCGTTTATTTAGTGCCGCGCGTCGCCATCAACTTCTTCTATTATCCTGACGATAAAATATACGGTccc      >  1:204537/1‑69 (MQ=255)

                                 tGCGTTTATTTAGTGCCGCGCGTCGCCATCAACTTCTTCTATTATCCTGACGATAAAATATACGGTccc      >  1:1966291/1‑69 (MQ=255)

                                  gCGTTTATTTAGTGCCGCGCGTTGCCATCAACttattc                                    >  1:3393775/1‑38 (MQ=255)

                                  gCGTTTATTTAGTGCCGCGCGTTGCCATCAACTtcttct                                   >  1:2351942/1‑39 (MQ=255)

                                  gCGTTTATTTAGTGCCGCGCGTTGCCATCAACTtcttct                                   >  1:1776420/1‑39 (MQ=255)

                                  gCGTTTATTTAGTGCCGCGCGTTGCCATCAACTtcttct                                   >  1:2109870/1‑39 (MQ=255)

                                     tttatttaGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATAAAATTTACGGTCCCg     <  1:9612/66‑1 (MQ=255)

                                     tttatttaGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATAAAATTTACGGTCCCGATc  <  1:1208331/69‑1 (MQ=255)

                                     tttatttaGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATAAAATTTACGGTCCCGATc  <  1:2558498/69‑1 (MQ=255)

                                     tttatttaGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATAAAATTTACGGTCCCGATc  <  1:1237241/69‑1 (MQ=255)

                                     tttatttaGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATAAAATTTACGGTCCCGATc  <  1:1608336/69‑1 (MQ=255)

                                     tttatttaGTGCCGCGCGTTGCCATCAACTTCTCCTATTATCCCGACGATAAAATTTACGGTCCCGATc  <  1:1858538/69‑1 (MQ=255)

AAAATTCTGTTTATTCTCTTTGTTGTTGCCTTCTGCGTTTATTTAGTGCCGCGCGTTGCCATCAACTTCTTCTATTATCCCGACGATAAAATTTACGGTCCCGATC  >  minE/1697106‑1697211

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`