BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,912,163	G→A	53.4%	D45D (GAC→GAT)	yjgF ←	ketoacid‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,912,163	0	G	A	53.4%	1.3 / 48.4	58	D45D (GAC→GAT)	yjgF	ketoacid‑binding protein
Reads supporting (aligned to +/- strand): ref base G (7/20); new base A (15/16); total (22/36)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

TTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCAT  >  minE/2912098‑2912222

ttCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCt                                                          >  1:568892/1‑69 (MQ=255)

ttCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCt                                                          >  1:286427/1‑69 (MQ=255)

ttCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCt                                                          >  1:2422315/1‑69 (MQ=255)

ttCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCt                                                          >  1:1602600/1‑69 (MQ=255)

ttCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATACAGCGACTGACGTGCCTGTGCAGCGACGTCt                                                          >  1:3209228/1‑69 (MQ=255)

     gCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCgg                                                     >  1:3064479/1‑69 (MQ=255)

     gCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCg                                                      <  1:1100706/68‑1 (MQ=255)

     gCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCg                                                      <  1:695488/68‑1 (MQ=255)

     gCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCg                                                      <  1:1547172/68‑1 (MQ=255)

        aGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTAc                                                  >  1:2154885/1‑69 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:2193522/69‑1 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:3023834/69‑1 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:308749/69‑1 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:845221/69‑1 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:2850770/69‑1 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:1274088/69‑1 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:921705/69‑1 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:3405521/69‑1 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:2168361/69‑1 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:923136/69‑1 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:957666/69‑1 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:1764430/69‑1 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:1763148/69‑1 (MQ=255)

                      tCGCTTTTACGTTATCCAGCGACTGACGCGCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:1696682/69‑1 (MQ=255)

                                 ttATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCggg                         <  1:2093153/69‑1 (MQ=255)

                                 ttATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCggg                         <  1:22101/69‑1 (MQ=255)

                                 ttATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCggg                         <  1:2297705/69‑1 (MQ=255)

                                 ttATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCggg                         <  1:2024928/69‑1 (MQ=255)

                                 ttATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCggg                         <  1:2533052/69‑1 (MQ=255)

                                 ttATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCggg                         <  1:1476298/69‑1 (MQ=255)

                                 ttATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCggg                         <  1:451206/69‑1 (MQ=255)

                                            cTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGa              <  1:906802/69‑1 (MQ=255)

                                            cTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGa              <  1:3156767/69‑1 (MQ=255)

                                            cTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGa              <  1:2310052/69‑1 (MQ=255)

                                            cTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGa              <  1:98050/69‑1 (MQ=255)

                                             tGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGa              <  1:3287850/68‑1 (MQ=255)

                                                cgcgCCTGTGCATCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:307721/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:1051782/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:762678/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:1130987/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:884817/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:1285626/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:523270/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:451211/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:3355514/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:30197/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:3010111/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:2562529/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:1417978/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:2014460/4‑60 (MQ=255)

                                                cgcgCCTGTGCAGCGACATCTGCCGGGACTTCGCCCGTTTTCGGGTTTACCGGGATCTGa                   >  1:1821308/4‑60 (MQ=255)

                                                 gTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGtg          <  1:3362917/68‑1 (MQ=255)

                                                 gTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGtg          <  1:28578/68‑1 (MQ=255)

                                                 gTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGtg          <  1:1560799/68‑1 (MQ=255)

                                                       gtgCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCa   <  1:774101/69‑1 (MQ=255)

                                                       gtgCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:3120707/36‑1 (MQ=37)

                                                       gtgCAGCGACATCTGCCGGGACTTCGCCCGTTTTCg                                    <  1:863150/36‑1 (MQ=37)

                                                          cAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCat  <  1:1395000/67‑1 (MQ=255)

TTCAGGCCAGCGGCTTCGACGATCGCTTTTACGTTATCCAGCGACTGACGTGCCTGTGCAGCGACGTCTGCCGGTACTTCGCCCGTTTTCGGATTTACCGGGATCTGACCGGAGGTGATGATCAT  >  minE/2912098‑2912222

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