BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	875,974	A→G	17.0%	N397S (AAC→AGC)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	875,974	0	A	G	17.0%	177.0 / 23.0	111	N397S (AAC→AGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand): ref base A (63/29); new base G (3/16); total (66/45)
Fisher's exact test for biased strand distribution p-value = 4.24e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

Only 100 of 111 total aligned reads displayed.

AGGGCCTAAACGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGAC  >  minE/875909‑876036

aGGGCCTAAACGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACg                                                              <  1:2990536/68‑1 (MQ=255)

aGGGCCTAAACGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACg                                                              <  1:1852890/68‑1 (MQ=255)

   gCCTAAACGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAAc                                                           >  1:2331971/1‑68 (MQ=255)

   gCCTAAACGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAAc                                                           >  1:755105/1‑68 (MQ=255)

   gCCTAAACGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAAc                                                           >  1:204842/1‑68 (MQ=255)

   gCCTAAACGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAAc                                                           >  1:1414594/1‑68 (MQ=255)

    ccTAAACGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACt                                                          >  1:1808578/1‑68 (MQ=255)

    ccTAAACGCC‑‑CAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGa                                                       >  1:1092182/1‑69 (MQ=255)

        aaCGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGaa                                                      >  1:1050913/1‑68 (MQ=255)

        aaCGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGaa                                                      >  1:2601012/1‑68 (MQ=255)

        aaCGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGaa                                                      >  1:2378638/1‑68 (MQ=255)

        aaCGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGaa                                                      >  1:2224163/1‑68 (MQ=255)

        aaCGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGaa                                                      >  1:2192274/1‑68 (MQ=255)

        aaCGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGaa                                                      >  1:2938764/1‑68 (MQ=255)

        aaCGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGaa                                                      >  1:546305/1‑68 (MQ=255)

        aaCGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGaa                                                      >  1:557814/1‑68 (MQ=255)

        aaCGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGaa                                                      >  1:1468976/1‑68 (MQ=255)

        aaCGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGaa                                                      >  1:1363181/1‑68 (MQ=255)

        aaCGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGaa                                                      >  1:1436909/1‑68 (MQ=255)

          cGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTg                                                   <  1:665975/69‑1 (MQ=255)

          cGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTg                                                   <  1:2752723/69‑1 (MQ=255)

          cGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTg                                                   <  1:417014/69‑1 (MQ=255)

          cGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTg                                                   <  1:1135862/69‑1 (MQ=255)

          cGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTg                                                   <  1:3048256/69‑1 (MQ=255)

          cGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTg                                                   <  1:803003/69‑1 (MQ=255)

          cGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTg                                                   <  1:3013900/69‑1 (MQ=255)

          cGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTg                                                   <  1:2171387/69‑1 (MQ=255)

          cGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTg                                                   <  1:1384197/69‑1 (MQ=255)

          cGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTg                                                   <  1:2137720/69‑1 (MQ=255)

           gCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTg                                                   <  1:1418770/68‑1 (MQ=255)

                 aGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGCACTGGAAGTGAATGCCa                                            >  1:2515368/1‑69 (MQ=255)

                 aGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGcc                                             >  1:2849679/1‑68 (MQ=255)

                 aGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCa                                            >  1:2812281/1‑69 (MQ=255)

                 aGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCa                                            >  1:1517867/1‑69 (MQ=255)

                 aGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCa                                            >  1:1369436/1‑69 (MQ=255)

                 aGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCa                                            >  1:3008521/1‑69 (MQ=255)

                 aGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCa                                            >  1:1760576/1‑69 (MQ=255)

                 aGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCa                                            >  1:381726/1‑69 (MQ=255)

                 aGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCa                                            >  1:2904461/1‑69 (MQ=255)

                  ggATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCa                                            >  1:2215169/1‑68 (MQ=255)

                  ggATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCa                                            >  1:1289391/1‑68 (MQ=255)

                  ggATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCa                                            >  1:2481744/1‑68 (MQ=255)

                    aTCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACg                                                              >  1:89376/1‑48 (MQ=255)

                    aTCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACg                                                              >  1:687621/1‑48 (MQ=255)

                    aTCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACg                                                              >  1:2499328/1‑48 (MQ=255)

                    aTCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACg                                                              >  1:2495277/1‑48 (MQ=255)

                    aTCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACg                                                              >  1:3170496/1‑48 (MQ=255)

                    aTCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACg                                                              >  1:2404090/1‑48 (MQ=255)

                    aTCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACg                                                              >  1:2819976/1‑48 (MQ=255)

                    aTCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACg                                                              >  1:3168797/1‑48 (MQ=255)

                    aTCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACg                                                              >  1:2850042/1‑48 (MQ=255)

                    aTCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACg                                                              >  1:1415430/1‑48 (MQ=255)

                            gCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGa                                                             >  1:936075/1‑41 (MQ=255)

                            gCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGa                                                             >  1:1872864/1‑41 (MQ=255)

                               cTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATc                                >  1:40582/1‑67 (MQ=255)

                               cTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATc                                >  1:1526574/1‑67 (MQ=255)

                               cTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATc                                >  1:627394/1‑67 (MQ=255)

                               cTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATc                                >  1:1373059/1‑67 (MQ=255)

                                  cccGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAg                           >  1:1353630/1‑69 (MQ=255)

                                  cccGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAg                           >  1:3054232/1‑69 (MQ=255)

                                  cccGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAg                           >  1:2971068/1‑69 (MQ=255)

                                  cccGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAg                           >  1:799154/1‑69 (MQ=255)

                                  cccGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAg                           >  1:2847854/1‑69 (MQ=255)

                                  cccGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAg                           >  1:167936/1‑69 (MQ=255)

                                  cccGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAg                           >  1:1087011/1‑69 (MQ=255)

                                  accGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCa                            >  1:1114150/2‑68 (MQ=255)

                                     gATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTg                                                   <  1:1939187/42‑1 (MQ=255)

                                           ccAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGAt                  <  1:3249151/69‑1 (MQ=255)

                                           ccAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGAt                  <  1:2366926/69‑1 (MQ=255)

                                           ccAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGAt                  <  1:528006/69‑1 (MQ=255)

                                                aGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGAttatgt             <  1:3199148/69‑1 (MQ=255)

                                                    tttGCTGCCAGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGAttatgttatg         >  1:2997486/1‑69 (MQ=255)

                                                    tttGCTGCCAGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGAttatgttatg         >  1:2351196/1‑69 (MQ=255)

                                                    tttGCTGCCAGTAGAGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGAttatgttat          >  1:1767336/1‑68 (MQ=255)

                                                      tGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGAt                  <  1:675746/58‑1 (MQ=255)

                                                          gccAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGa   <  1:2059948/69‑1 (MQ=255)

                                                          gccAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGa   <  1:2065580/69‑1 (MQ=255)

                                                          gccAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGa   <  1:1606530/69‑1 (MQ=255)

                                                          gccAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGa   <  1:1607026/69‑1 (MQ=255)

                                                          gccAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGa   <  1:2391097/69‑1 (MQ=255)

                                                          gccAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGa   <  1:2466605/69‑1 (MQ=255)

                                                          gccAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGa   <  1:2557311/69‑1 (MQ=255)

                                                          gccAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGa   <  1:1650062/69‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:3110712/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:1433024/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:535533/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:3430555/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:3266315/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:1684701/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:262801/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:3079110/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:3010666/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:2936390/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:2091115/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:2867991/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:2850672/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:2720807/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATGAACGGAc  <  1:2632515/67‑1 (MQ=255)

                                                             aGTAGCGAACTGGAAGTGAATGCCACGCATAAAGATCGCCTGCCGGTCGATTATGTTATAAACGGAc  <  1:1724625/67‑1 (MQ=255)

                                                             aGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGAt                  <  1:895598/51‑1 (MQ=255)

AGGGCCTAAACGCCCACAGGATCGCGTTGCACTGCCCGATGTACCAAAAGCATTTGCCGCCAGTAACGAACTGGAAGTGAATGCCACGCATAAAGATCGCCAGCCGGTCGATTATGTTATGAACGGAC  >  minE/875909‑876036

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