BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	511,688	T→A	33.2%	S577S (TCT→TCA)	ybhJ →	predicted hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	511,688	0	T	A	33.2%	58.1 / 38.6	57	S577S (TCT→TCA)	ybhJ	predicted hydratase
Reads supporting (aligned to +/- strand): ref base T (23/15); new base A (11/8); total (34/23)
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

GATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGC  >  minE/511622‑511753

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

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

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

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

gATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTcg                                                                 >  1:1554261/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: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:2297477/1‑69 (MQ=255)

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

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

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

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

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

                aCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTa                                                  >  1:1179/1‑68 (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:386894/69‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                               gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGc  >  1:183761/1‑69 (MQ=255)

                                                               gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGc  >  1:2155524/1‑69 (MQ=255)

                                                               gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGc  >  1:2203150/1‑69 (MQ=255)

                                                               gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGc  >  1:3408514/1‑69 (MQ=255)

                                                               gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGc  >  1:1252590/1‑69 (MQ=255)

                                                               gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGc  >  1:1815443/1‑69 (MQ=255)

                                                               gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGc  >  1:2530307/1‑69 (MQ=255)

                                                               gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGc  >  1:1616532/1‑69 (MQ=255)

                                                               gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGc  >  1:2685142/1‑69 (MQ=255)

GATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGC  >  minE/511622‑511753

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