BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,649,709	A→G	43.7%	A310A (GCA→GCG)	yfgC →	predicted peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,649,709	0	A	G	43.7%	25.5 / 90.4	112	A310A (GCA→GCG)	yfgC	predicted peptidase
Reads supporting (aligned to +/- strand): ref base A (40/23); new base G (13/36); total (53/59)
Fisher's exact test for biased strand distribution p-value = 1.27e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

AGCTCACCAGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTAA  >  minE/1649651‑1649768

aGCTCACCAGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATg                                                   <  1:1791723/69‑1 (MQ=255)

aGCTCACCAGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATg                                                   <  1:2286038/69‑1 (MQ=255)

aGCTCACCAGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATg                                                   <  1:1306629/69‑1 (MQ=255)

aGCTCACCAGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATg                                                   <  1:877009/69‑1 (MQ=255)

aGCTCACCAGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATg                                                   <  1:2365256/69‑1 (MQ=255)

 gCTCACCAGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATg                                                   <  1:3406800/68‑1 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTg                                             >  1:708936/1‑67 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGc                                            >  1:1690189/1‑68 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGc                                            >  1:968699/1‑68 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGc                                            >  1:2308681/1‑68 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGc                                            >  1:2633107/1‑68 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGc                                            >  1:467652/1‑68 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGc                                            >  1:2869092/1‑68 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGc                                            >  1:3133461/1‑68 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGc                                            >  1:3275083/1‑68 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGc                                            >  1:3408152/1‑68 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGc                                            >  1:3311802/1‑68 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGc                                            >  1:1464410/1‑68 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGc                                            >  1:117198/1‑68 (MQ=255)

        aGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGc                                            >  1:946968/1‑68 (MQ=255)

            aTTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTa                                        >  1:2528403/1‑68 (MQ=255)

            aTTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTAc                                       >  1:1682529/1‑69 (MQ=255)

            aTTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTAc                                       >  1:793950/1‑69 (MQ=255)

            aTTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTAc                                       >  1:3407371/1‑69 (MQ=255)

            aTTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTAc                                       >  1:3034431/1‑69 (MQ=255)

             tttGCTGGATTAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:446250/52‑1 (MQ=255)

             tttGCTGGATGAATGGGCGAAAGGAAAGGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:2277300/52‑1 (MQ=38)

             tttGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:544176/52‑1 (MQ=39)

             tttGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:3109176/52‑1 (MQ=39)

             tttGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:42834/52‑1 (MQ=39)

             tttGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:1392642/52‑1 (MQ=39)

             tttGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:3000398/52‑1 (MQ=39)

             tttGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:449832/52‑1 (MQ=39)

             tttGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:1327358/52‑1 (MQ=39)

             tttGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:373486/52‑1 (MQ=39)

             tttGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:1216178/52‑1 (MQ=39)

             tttGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:194365/52‑1 (MQ=39)

             tttGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:718096/52‑1 (MQ=39)

             tttGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:1663148/52‑1 (MQ=39)

                        aaTGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACaa                                                       <  1:395444/41‑1 (MQ=255)

                        aaTGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACaa                                                       <  1:887192/41‑1 (MQ=255)

                        aaTGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACaa                                                       <  1:1629352/41‑1 (MQ=255)

                        aaTGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACaa                                                       <  1:1521010/41‑1 (MQ=255)

                          tGGGCGAAAGGAAACGTTCGTCAGCAACGAGCGGCGCaa                                                       <  1:533255/39‑1 (MQ=38)

                          tGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGc                          <  1:1554254/68‑1 (MQ=255)

                          tGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGc                          <  1:948840/68‑1 (MQ=255)

                          tGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGc                          <  1:3373816/68‑1 (MQ=255)

                          tGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGc                          <  1:1482523/68‑1 (MQ=255)

                          tGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGc                          <  1:2421348/68‑1 (MQ=255)

                          tGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGc                          <  1:185702/68‑1 (MQ=255)

                          tGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGc                          <  1:2352120/68‑1 (MQ=255)

                          tGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGc                          <  1:575348/68‑1 (MQ=255)

                          tGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCCTGCTTTACAGGCGATGGAAGc                          <  1:743304/68‑1 (MQ=255)

                           gggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGcc                         <  1:504734/68‑1 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGcc                         >  1:3032425/1‑67 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGcc                         >  1:3206953/1‑67 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGcc                         >  1:660747/1‑67 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCa                        >  1:2620385/1‑68 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCa                        >  1:1747165/1‑68 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCa                        >  1:713983/1‑68 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCa                        >  1:3283894/1‑68 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCa                        >  1:3147824/1‑68 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCa                        >  1:675829/1‑68 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCa                        >  1:2127483/1‑68 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCa                        >  1:3094390/1‑68 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCa                        >  1:1410285/1‑68 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCa                        >  1:3028253/1‑68 (MQ=255)

                            ggCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCa                        >  1:2480707/1‑68 (MQ=255)

                               gAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATg                                                   <  1:1361914/38‑1 (MQ=255)

                                caaGGAAACGTTCGTCAGCAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGa                                 >  1:2542593/2‑55 (MQ=255)

                                aaaGGAAACGTTCGTCAGCAACGAGCGGCGCAATATGGTCGTGCTTTACAGGc                                   >  1:1017806/1‑53 (MQ=255)

                                aaaGGAAACGTTCGTCAGCAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCg                                  >  1:3119058/1‑54 (MQ=255)

                                aaaGGAAACGTTCGTCAGCAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGa                                 >  1:851908/1‑55 (MQ=255)

                                aaaGGAAACGTTCGTCAGCAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGa                                 >  1:1836550/1‑55 (MQ=255)

                                aaaGGAAACGTTCGTCAGCAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGa                                 >  1:2098132/1‑55 (MQ=255)

                                aaaGGAAACGTTCGTCAGCAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGa                                 >  1:1438156/1‑55 (MQ=255)

                                aaaGGAAACGTTCGTCAGCAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGa                                 >  1:2159768/1‑55 (MQ=255)

                                aaaGGAAACGTTCGTCAGCAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGa                                 >  1:2276024/1‑55 (MQ=255)

                                aaaGGAAACGTTCGTCAGCAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGa                                 >  1:832001/1‑55 (MQ=255)

                                aaaGGAAACGTTCGTCAGCAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGa                                 >  1:2749832/1‑55 (MQ=255)

                                aaaGGAAACGTTCGTCAGCAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGa                                 >  1:329835/1‑55 (MQ=255)

                                aaaGGAAACGTTCGTCAGCAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGa                                 >  1:1015542/1‑55 (MQ=255)

                                  aGGAAACGTTCGTCAGCAACGAGCAGCACAATATGTTCGTGCTTTACAGGCGATGGAAGCCaataaat                  >  1:1213304/1‑68 (MQ=255)

                                  aGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCaataaat                  >  1:2865632/1‑68 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:588481/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:1312624/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:1420848/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:1118043/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:1726545/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:1927835/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:2042203/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:339969/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:2142919/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:2388215/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:2668791/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:3040052/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:3067843/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:3420235/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTaa  <  1:3357923/68‑1 (MQ=255)

                                                  cAACGAGCGGCGCAATATGGTCGTGCTTTACAGGCGATGCAAGCCAATAAATACGACGAGGCGCGTaa  <  1:1419343/68‑1 (MQ=255)

AGCTCACCAGTGATTTGCTGGATGAATGGGCGAAAGGAAACGTTCGTCAGCAACGAGCAGCACAATATGGTCGTGCTTTACAGGCGATGGAAGCCAATAAATACGACGAGGCGCGTAA  >  minE/1649651‑1649768

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