BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,829,222	A→G	31.0%	V256V (GTT→GTC)	aspA ←	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,829,222	0	A	G	31.0%	36.4 / 19.4	45	V256V (GTT→GTC)	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand): ref base A (8/23); new base G (8/6); total (16/29)
Fisher's exact test for biased strand distribution p-value = 5.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGC  >  minE/2829155‑2829290

taacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTaa                                                                     >  1:318196/1‑69 (MQ=255)

taacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTaa                                                                     >  1:13/1‑69 (MQ=255)

taacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTaa                                                                     >  1:707573/1‑69 (MQ=255)

taacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTaa                                                                     >  1:1710500/1‑69 (MQ=255)

                gTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGtttttt                                                     <  1:1979563/69‑1 (MQ=255)

                gTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGtttttt                                                     <  1:1207608/69‑1 (MQ=255)

                gTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGtttttt                                                     <  1:2450402/69‑1 (MQ=255)

                gTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGtttttt                                                     <  1:1738069/69‑1 (MQ=255)

                gTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGtttttt                                                     <  1:1711864/69‑1 (MQ=255)

                 tCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGtttttt                                                     <  1:1385825/68‑1 (MQ=255)

                         cacTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAg                                            >  1:2819218/3‑69 (MQ=255)

                         cGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAg                                            >  1:1396266/1‑69 (MQ=255)

                             tCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTGACTTCAGCCAGTTTTTTCACTGCCAGCg                                          >  1:382526/1‑67 (MQ=255)

                             tCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTGACTTCAGCCAGTTTTTTCACTGCCAGCg                                          >  1:214277/1‑67 (MQ=255)

                             tCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTGACTTCAGCCAGTTTTTTCACTGCCAGCg                                          >  1:380147/1‑67 (MQ=255)

                             tCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTGACTTCAGCCAGTTTTTTCACTGCCAGCg                                          >  1:482710/1‑67 (MQ=255)

                             tCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTGACTTCAGCCAGTTTTTTCACTGCCAGCg                                          >  1:3211678/1‑67 (MQ=255)

                             tCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTGACTTCAGCCAGTTTTTTCACTGCCAGCg                                          >  1:912114/1‑67 (MQ=255)

                             tCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTGACTTCAGCCAGTTTTTTCACTGCCAGCg                                          >  1:732353/1‑67 (MQ=255)

                             tCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTGACTTCAGCCAGTTTTTTCACTGCCAGCg                                          >  1:166603/1‑67 (MQ=255)

                                aTCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGgag                                       <  1:2402126/67‑1 (MQ=255)

                                aTCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGgag                                       <  1:2289120/67‑1 (MQ=255)

                                 tCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGagag                                     <  1:2257128/68‑1 (MQ=255)

                                 tCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGagag                                     <  1:224542/68‑1 (MQ=255)

                                 tCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGagag                                     <  1:364250/68‑1 (MQ=255)

                                  cAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTa                                   <  1:1964735/69‑1 (MQ=255)

                                  cAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTa                                   <  1:76039/69‑1 (MQ=255)

                                  cAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTa                                   <  1:3293070/69‑1 (MQ=255)

                                  cAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGCTTTTTCACTGCCAGCGGagag                                     <  1:2543772/67‑1 (MQ=255)

                                          cAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGtttttt                                                     <  1:1714992/43‑1 (MQ=255)

                                                  aCGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCgtg                      <  1:1483949/66‑1 (MQ=255)

                                                  aCGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCgtg                      <  1:334558/66‑1 (MQ=255)

                                                  aCGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCgtg                      <  1:2797381/66‑1 (MQ=255)

                                                  aCGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCgtg                      <  1:1195114/66‑1 (MQ=255)

                                                   cGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCgtg                      <  1:309883/65‑1 (MQ=255)

                                                       tGGGAAGCCAGTAACTTCAGCCAGTTTTTTCATTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGAcc              <  1:1181656/69‑1 (MQ=255)

                                                         ggAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGga                                        >  1:2847416/1‑41 (MQ=255)

                                                         ggAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGga                                        >  1:3314338/1‑41 (MQ=255)

                                                         ggAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGga                                        >  1:2180476/1‑41 (MQ=255)

                                                         ggAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGga                                        >  1:74182/1‑41 (MQ=255)

                                                         ggAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGga                                        >  1:1146947/1‑41 (MQ=255)

                                                         ggAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGga                                        >  1:999205/1‑41 (MQ=255)

                                                              ccAGTGACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCg       <  1:2348471/69‑1 (MQ=255)

                                                              ccAGTGACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCg       <  1:2327689/69‑1 (MQ=255)

                                                              ccAGTGACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCg       <  1:1997762/69‑1 (MQ=255)

                                                              ccAGTGACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCg       <  1:575507/69‑1 (MQ=255)

                                                              ccAGTGACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCg       <  1:935107/69‑1 (MQ=255)

                                                              ccAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCgtg                      <  1:476119/54‑1 (MQ=255)

                                                               cAGTGACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCg       <  1:2665846/68‑1 (MQ=255)

                                                                   aaCTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGc  >  1:182437/1‑69 (MQ=255)

                                                                   aaCTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGc  >  1:1355259/1‑69 (MQ=255)

                                                                   aaCTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGc  >  1:1236057/1‑69 (MQ=255)

TAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGC  >  minE/2829155‑2829290

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