BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	383,076	C→G	32.1%	T162T (ACC→ACG)	entB →	isochorismatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	383,076	0	C	G	32.1%	120.8 / 79.6	109	T162T (ACC→ACG)	entB	isochorismatase
Reads supporting (aligned to +/- strand): ref base C (33/41); new base G (12/23); total (45/64)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 109 total aligned reads displayed.

GAGTGGACGTAACCAGCTGATTATTACCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAT  >  minE/383010‑383130

gagTGGACGTAACCAGCTGATTATTACCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGa                                                      <  1:916351/69‑1 (MQ=255)

       cGTAACCAGCTGATTATTACCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCAttt                                               >  1:2873065/1‑69 (MQ=255)

       cGTAACCAGCTGATTATTACCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCAtt                                                >  1:1897574/1‑68 (MQ=255)

          aaCCAGCTGATTATTACCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCAtt                                                <  1:2096471/65‑1 (MQ=255)

                     tattaCCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTaa                                 <  1:1816791/69‑1 (MQ=255)

                     tattaCCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTaa                                 <  1:1532339/69‑1 (MQ=255)

                     tattaCCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTaa                                 <  1:1771138/69‑1 (MQ=255)

                     tattaCCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTaa                                 <  1:3077753/69‑1 (MQ=255)

                     tattaCCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTaa                                 <  1:3218827/69‑1 (MQ=255)

                     tattaCCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTa                                  <  1:3001883/68‑1 (MQ=255)

                     tattaCCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTa                                  <  1:1250612/68‑1 (MQ=255)

                     tattaCCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTa                                  <  1:1768804/68‑1 (MQ=255)

                     tattaCCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTa                                  <  1:2439856/68‑1 (MQ=255)

                     tattaCCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTa                                  <  1:2435974/68‑1 (MQ=255)

                     tattaCCGGGGTATATGCCCACATTGGCTGTATGACCACCGAAACCGACGCATTTATGCGCGATATTaa                                 <  1:208856/69‑1 (MQ=255)

                     tattaCCGGCGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTaa                                 <  1:3222730/69‑1 (MQ=255)

                      attaCCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTa                                  <  1:896289/67‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:1951055/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:1951308/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:598935/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:2172705/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:2314886/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:2491107/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:921667/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:455174/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:3046432/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:3157251/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:3259136/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:3429841/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:1079664/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:1789076/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:1768806/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:1280384/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:1302524/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:869938/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:1536212/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:1245238/69‑1 (MQ=255)

                            ggggTATATGCCCACATTGGCTGTATGACCACCGCAACGGACCCATTTATGCGCGATATTAAACCGttt                          <  1:1937201/69‑1 (MQ=255)

                                  tatGCCCACATTGGCTGTATGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:724633/63‑1 (MQ=255)

                                  tatGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAtggtg                    >  1:869879/1‑69 (MQ=255)

                                    tGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAtggtgg                   >  1:3397804/1‑68 (MQ=255)

                                    tGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAtggtgg                   >  1:3340493/1‑68 (MQ=255)

                                    tGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAtggtgg                   >  1:2343069/1‑68 (MQ=255)

                                    tGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAtggtgg                   >  1:1127331/1‑68 (MQ=255)

                                    tGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAtggtg                    >  1:1128339/1‑67 (MQ=255)

                                    tGCCCACATTGGCTGTATGACCACCGCAACCGACGAATTTATGCGCGATATTAAACCGTTTAtggtgg                   >  1:1306481/1‑68 (MQ=255)

                                           aTTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcg           <  1:538695/69‑1 (MQ=255)

                                           aTTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcg           <  1:486355/69‑1 (MQ=255)

                                           aTTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcg           <  1:590273/69‑1 (MQ=255)

                                           aTTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcg           <  1:2821574/69‑1 (MQ=255)

                                           aTTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcg           <  1:388398/69‑1 (MQ=255)

                                           aTTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcg           <  1:252679/69‑1 (MQ=255)

                                           aTTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcg           <  1:1639301/69‑1 (MQ=255)

                                           aTTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcg           <  1:2094612/69‑1 (MQ=255)

                                            ttGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcgc          >  1:1138605/1‑69 (MQ=255)

                                            ttGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcgc          >  1:1151598/1‑69 (MQ=255)

                                            ttGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcgc          >  1:2975624/1‑69 (MQ=255)

                                            ttGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcgc          >  1:2166357/1‑69 (MQ=255)

                                            ttGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcgc          >  1:2974727/1‑69 (MQ=255)

                                              ggCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCg                 <  1:617267/60‑1 (MQ=255)

                                              ggCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCg                 <  1:207282/60‑1 (MQ=255)

                                              ggCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCg                 <  1:3200878/60‑1 (MQ=255)

                                              ggCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCg                 <  1:2920961/60‑1 (MQ=255)

                                              ggCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCg                 <  1:536552/60‑1 (MQ=255)

                                              ggCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCg                 <  1:2288008/60‑1 (MQ=255)

                                              ggCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCg                 <  1:3312402/60‑1 (MQ=255)

                                              ggCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCg                 <  1:1079206/60‑1 (MQ=255)

                                              ggCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCg                 <  1:2576777/60‑1 (MQ=255)

                                              ggCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTg        >  1:1134659/1‑69 (MQ=255)

                                              ggCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGAGGCg                 <  1:1514550/60‑1 (MQ=255)

                                               gCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCg                 <  1:414850/59‑1 (MQ=255)

                                                   tATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGa   >  1:227555/1‑69 (MQ=255)

                                                     tGACCACCTCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:1399690/1‑68 (MQ=255)

                                                     tGACCACCGCGACCGACGCATTTATGCGCGATATTAAACCGTTTAt                        >  1:1798662/1‑46 (MQ=255)

                                                     tGACCACCGCAACGGACGCATTTATGCGCGATATTAAAGCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:2299285/1‑68 (MQ=255)

                                                     tGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:78366/1‑68 (MQ=255)

                                                     tGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:983411/1‑68 (MQ=255)

                                                     tGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:831617/1‑68 (MQ=255)

                                                     tGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:291443/1‑68 (MQ=255)

                                                     tGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:105656/1‑68 (MQ=255)

                                                     tGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:270554/1‑68 (MQ=255)

                                                     tGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:2638209/1‑68 (MQ=255)

                                                     tGACCACCGCAACGGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:2439406/1‑68 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAt                        >  1:3038178/1‑46 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAt                        >  1:2958946/1‑46 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAt                        >  1:2907641/1‑46 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAt                        <  1:559475/46‑1 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAt                        >  1:2652445/1‑46 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAt                        >  1:68981/1‑46 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAt                        >  1:707303/1‑46 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAt                        >  1:746459/1‑46 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGa   >  1:1157275/1‑67 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:474710/1‑68 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:1298584/1‑68 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:1815257/1‑68 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAt  >  1:2839573/1‑68 (MQ=255)

                                                     tGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTCGCGGATGCGCTGGCCGAt  >  1:636248/1‑68 (MQ=255)

                                                      gACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTAt                        <  1:1133989/45‑1 (MQ=255)

                                                         caccGCAACGGACGCATTTATGCGCGATATTAAACCGttt                          <  1:1030157/40‑1 (MQ=255)

                                                                aCCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATgcg           <  1:1025226/48‑1 (MQ=255)

GAGTGGACGTAACCAGCTGATTATTACCGGGGTATATGCCCACATTGGCTGTATGACCACCGCAACCGACGCATTTATGCGCGATATTAAACCGTTTATGGTGGCGGATGCGCTGGCCGAT  >  minE/383010‑383130

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