BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	511,682	C→T	38.7%	T575T (ACC→ACT)	ybhJ →	predicted hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	511,682	0	C	T	38.7%	22.9 / 29.1	49	T575T (ACC→ACT)	ybhJ	predicted hydratase
Reads supporting (aligned to +/- strand): ref base C (15/15); new base T (11/8); total (26/23)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

GATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGC  >  minE/511622‑511729

gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGtctc                                          >  1:364335/1‑68 (MQ=255)

gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTcg                                         >  1:109350/1‑69 (MQ=255)

gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTcg                                         >  1:1554261/1‑69 (MQ=255)

gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTcg                                         >  1:2325798/1‑69 (MQ=255)

gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTcg                                         >  1:2135118/1‑69 (MQ=255)

             gAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGt                            >  1:1786213/1‑69 (MQ=255)

             gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCAcgccgat                                    >  1:1412195/1‑59 (MQ=255)

             gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAggg                            >  1:3176068/1‑68 (MQ=255)

             gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAgg                             >  1:1107020/1‑68 (MQ=255)

             gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt                            >  1:2297477/1‑69 (MQ=255)

             gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt                            >  1:44516/1‑69 (MQ=255)

             gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt                            >  1:2921522/1‑69 (MQ=255)

             gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt                            >  1:2762203/1‑69 (MQ=255)

             gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt                            >  1:2550664/1‑69 (MQ=255)

             gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt                            >  1:2438641/1‑69 (MQ=255)

             gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt                            >  1:24264/1‑69 (MQ=255)

             gAAACCTCTTCATATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGt                            >  1:1934504/1‑69 (MQ=255)

                aCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTa                          <  1:2224982/68‑1 (MQ=255)

                aCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTa                          <  1:491834/68‑1 (MQ=255)

                aCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTa                          >  1:1179/1‑68 (MQ=255)

                aCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTa                          <  1:2559533/68‑1 (MQ=255)

                aCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTa                          <  1:340364/68‑1 (MQ=255)

                   tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt                      <  1:209880/69‑1 (MQ=255)

                   tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt                      <  1:533397/69‑1 (MQ=255)

                   tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt                      <  1:2189004/69‑1 (MQ=255)

                   tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt                      <  1:2220091/69‑1 (MQ=255)

                   tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt                      <  1:386894/69‑1 (MQ=255)

                   tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt                      <  1:1184506/69‑1 (MQ=255)

                   tcttctTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGtt                      <  1:1644619/69‑1 (MQ=255)

                       cttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAg                    <  1:1677678/67‑1 (MQ=255)

                       cttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAg                    <  1:2797675/67‑1 (MQ=255)

                       cttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAg                    <  1:2108237/67‑1 (MQ=255)

                                        aTTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGTTATGTTGGCAGAAGTAAAGCGActgc  <  1:723923/68‑1 (MQ=255)

                                        aTTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGTTATGTTGGCAGAAGTAAAGCGActgc  <  1:3069254/68‑1 (MQ=255)

                                        aTTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGTTATGTTGGCAGAAGTAAAGCGActgc  <  1:1444596/68‑1 (MQ=255)

                                        aTTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGTTATGTTGGCAGAAGTAAAGCGActgc  <  1:2906559/68‑1 (MQ=255)

                                        aTTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGTTATGTTGGCAGAAGTAAAGCGActgc  <  1:2789935/68‑1 (MQ=255)

                                        aTTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGTTATGTTGGCAGAAGTAAAGCGActgc  <  1:1447213/68‑1 (MQ=255)

                                        aTTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGTTATGTTGGCAGAAGTAAAGCGActgc  <  1:155108/68‑1 (MQ=255)

                                        aTTGGTCTGGCGGAGTTTACTCTGTCACGCCGCGATCCAGGTTATGTTGGCAGAAGTAAAGCGActgc  <  1:2127690/68‑1 (MQ=255)

                                              cTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCg       <  1:29598/57‑1 (MQ=255)

                                                         tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGc        >  1:2869090/1‑45 (MQ=255)

                                                         tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGc        >  1:2749454/1‑45 (MQ=255)

                                                         tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGc        >  1:2569035/1‑45 (MQ=255)

                                                         tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGc        >  1:1828658/1‑45 (MQ=255)

                                                         tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGc        >  1:1175288/1‑45 (MQ=255)

                                                         tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGc        >  1:1076560/1‑45 (MQ=255)

                                                         tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGc        >  1:546524/1‑45 (MQ=255)

                                                         tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGc        >  1:2080624/1‑45 (MQ=255)

GATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGC  >  minE/511622‑511729

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