BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,244,760	A→T	41.8%	A264A (GCT→GCA)	yebA ←	predicted peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,244,760	0	A	T	41.8%	74.6 / 156.2	153	A264A (GCT→GCA)	yebA	predicted peptidase
Reads supporting (aligned to +/- strand): ref base A (49/40); new base T (15/49); total (64/89)
Fisher's exact test for biased strand distribution p-value = 1.18e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 153 total aligned reads displayed.

GCAAGAATCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAC  >  minE/1244696‑1244824

gCAAGAATCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGt                                                              >  1:559941/1‑69 (MQ=255)

gCAAGAATCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGt                                                              >  1:730404/1‑69 (MQ=255)

      aTCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTaataat                                                        >  1:2860880/1‑69 (MQ=255)

      aTCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTaataat                                                        >  1:2093460/1‑69 (MQ=255)

      aTCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTaataat                                                        >  1:909918/1‑69 (MQ=255)

      aTCCTTTCGCCAGACAAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTaataat                                                        >  1:1034055/1‑69 (MQ=255)

        ccTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATc                                                       >  1:1867327/1‑68 (MQ=255)

        ccTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATc                                                       >  1:666401/1‑68 (MQ=255)

        ccTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATc                                                       >  1:2432685/1‑68 (MQ=255)

         cTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCtt                                                     <  1:583670/69‑1 (MQ=255)

         cTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCtt                                                     <  1:1042544/69‑1 (MQ=255)

         cTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCtt                                                     <  1:1038660/69‑1 (MQ=255)

         cTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCtt                                                     <  1:330239/69‑1 (MQ=255)

            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTa                                                   >  1:503704/1‑68 (MQ=255)

            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:1948166/1‑69 (MQ=255)

            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:2078260/1‑69 (MQ=255)

            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:1673656/1‑69 (MQ=255)

            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:1657130/1‑69 (MQ=255)

            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:1734489/1‑69 (MQ=255)

            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:584212/1‑69 (MQ=255)

            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:689093/1‑69 (MQ=255)

            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTAc                                                  >  1:753058/1‑69 (MQ=255)

            tCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGAAATAATCTTTAc                                                  >  1:1448665/1‑69 (MQ=255)

                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:2825418/67‑1 (MQ=255)

                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:2466833/67‑1 (MQ=255)

                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:3239043/67‑1 (MQ=255)

                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:2269874/67‑1 (MQ=255)

                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:3336014/67‑1 (MQ=255)

                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:3402083/67‑1 (MQ=255)

                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:1966162/67‑1 (MQ=255)

                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:940434/67‑1 (MQ=255)

                    ccAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:1636168/67‑1 (MQ=255)

                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:1241028/1‑69 (MQ=255)

                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:2436333/1‑69 (MQ=255)

                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:2602487/1‑69 (MQ=255)

                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:1633723/1‑69 (MQ=255)

                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:3264868/1‑69 (MQ=255)

                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:905006/1‑69 (MQ=255)

                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:1707965/1‑69 (MQ=255)

                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:2164984/1‑69 (MQ=255)

                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:2109491/1‑69 (MQ=255)

                         aCCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaa                                     >  1:605553/1‑69 (MQ=255)

                          ccGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaaa                                    >  1:321655/1‑69 (MQ=255)

                          ccGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCaaa                                    >  1:1949683/1‑69 (MQ=255)

                            gTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACg                                  <  1:3054522/69‑1 (MQ=255)

                               aCGATCGTAGAACTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTAc                               >  1:1169282/1‑69 (MQ=255)

                               aCGATCGTAGAACTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTAc                               >  1:1697294/1‑69 (MQ=255)

                               aCGATCGTAGAACTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTAc                               >  1:1482159/1‑69 (MQ=255)

                                cGATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTAATAATCTTTACCTTCTg                                            <  1:1484796/55‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:2948949/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:3098962/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:2925066/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:2559859/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:1540651/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:3233246/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:1166556/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:328234/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:434696/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:445238/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:127837/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:510836/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:562129/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:222989/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:173778/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:1853376/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:1110124/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:1650590/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:266073/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:2402183/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:2441525/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:819321/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:255404/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTGCGTaataat                                                        <  1:2651372/42‑1 (MQ=255)

                                 gATCGTAGAATTTGCCATCTTCAGCGCGGATTCCGTaataat                                                        <  1:2405485/42‑1 (MQ=38)

                                             tGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGc                 >  1:1703813/1‑69 (MQ=255)

                                                cATCTTGAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACg                                  <  1:2100809/49‑1 (MQ=255)

                                                    ttCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTc          >  1:1118077/1‑69 (MQ=255)

                                                    ttCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTc          >  1:131986/1‑69 (MQ=255)

                                                    ttCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTc          >  1:1217996/1‑69 (MQ=255)

                                                    ttCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTc          >  1:1878939/1‑69 (MQ=255)

                                                    ttCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTc          >  1:180400/1‑69 (MQ=255)

                                                         cgcgGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGtt     <  1:2425197/69‑1 (MQ=255)

                                                         cgcgGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGtt     <  1:1417341/69‑1 (MQ=255)

                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:113202/69‑1 (MQ=255)

                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:1740448/69‑1 (MQ=255)

                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:1135828/69‑1 (MQ=255)

                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:554350/69‑1 (MQ=255)

                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:539922/69‑1 (MQ=255)

                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:344529/69‑1 (MQ=255)

                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:1458076/69‑1 (MQ=255)

                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:3286485/69‑1 (MQ=255)

                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:3265404/69‑1 (MQ=255)

                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:2512465/69‑1 (MQ=255)

                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:2501010/69‑1 (MQ=255)

                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:1544470/69‑1 (MQ=255)

                                                            ggATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCCGGCTTTGCTCACGTTTAc  <  1:1914997/69‑1 (MQ=255)

                                                             gATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAc  <  1:2590691/68‑1 (MQ=255)

                                                               tAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGtt     >  1:523302/1‑63 (MQ=255)

                                                               tAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGtt     >  1:1733748/1‑63 (MQ=255)

                                                               tAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGtt     >  1:1010475/1‑63 (MQ=255)

GCAAGAATCCTTTCGCCAGACCAGTACCGTTACGATCGTAGAATTTGCCATCTTCAGCGCGGATAGCGTAATAATCTTTACCTTCTGAACGCAAACGTACGCCCAGCAGCTGGCTTTGCTCACGTTTAC  >  minE/1244696‑1244824

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