BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,246,353	T→C	33.3%	T49T (ACA→ACG)	znuA ←	zinc transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,246,353	0	T	C	33.3%	29.8 / 21.5	45	T49T (ACA→ACG)	znuA	zinc transporter subunit
Reads supporting (aligned to +/- strand): ref base T (10/20); new base C (2/13); total (12/33)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAACCC  >  minE/1246290‑1246394

ttttACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCtctgtttc                                      <  1:1880761/69‑1 (MQ=255)

    aCATCCGATGGGCGCAGTGAATAATCATGTTTTGAAGCGCCGTCAGGAAGTAAAACCtctgtttctgt                                   <  1:2668689/68‑1 (MQ=255)

    aCATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCtctgtttctgt                                   <  1:2064587/68‑1 (MQ=255)

        ccGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTAccc                              <  1:2307036/69‑1 (MQ=255)

        ccGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTAccc                              <  1:1220211/69‑1 (MQ=255)

         cGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTAcccc                             >  1:1446423/1‑69 (MQ=255)

         cGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTAcccc                             >  1:2717327/1‑69 (MQ=255)

                cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCGATCAGc                       <  1:439987/68‑1 (MQ=255)

                cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:328388/68‑1 (MQ=255)

                cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:2391271/68‑1 (MQ=255)

                cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:2558948/68‑1 (MQ=255)

                cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:2212824/68‑1 (MQ=255)

                cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:2194651/68‑1 (MQ=255)

                cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:1325376/68‑1 (MQ=255)

                cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:1160504/68‑1 (MQ=255)

                cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:888729/68‑1 (MQ=255)

                cgcAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:1038728/68‑1 (MQ=255)

                      gAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:19539/62‑1 (MQ=255)

                      gAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:627760/62‑1 (MQ=255)

                       aataatCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:3394776/61‑1 (MQ=255)

                         taatCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCa                      >  1:188292/1‑60 (MQ=255)

                         taatCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCa                      >  1:646836/1‑60 (MQ=255)

                         taatCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCa                      >  1:1019831/1‑60 (MQ=255)

                         taatCAAGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCa                      >  1:1022/1‑60 (MQ=255)

                           atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2472727/69‑1 (MQ=255)

                           atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2441270/69‑1 (MQ=255)

                           atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:842386/69‑1 (MQ=255)

                           atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:838085/69‑1 (MQ=255)

                           atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2045452/69‑1 (MQ=255)

                           atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:272318/69‑1 (MQ=255)

                           atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2875961/69‑1 (MQ=255)

                           atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:721927/69‑1 (MQ=255)

                           atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:691073/69‑1 (MQ=255)

                           atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:1153876/69‑1 (MQ=255)

                           atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:1254785/69‑1 (MQ=255)

                           atCATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:1301558/69‑1 (MQ=255)

                             cATGTTCTGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2931890/67‑1 (MQ=255)

                                   cTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGc                       <  1:1491809/49‑1 (MQ=255)

                                    tGAAGCCCCGTCAGGTAGTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAAccc  >  1:2812376/1‑69 (MQ=255)

                                    tGAAGCCCCGTCAGGTAGTAAAACCTACGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAAccc  >  1:1064795/1‑69 (MQ=255)

                                                aGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           >  1:451971/1‑48 (MQ=255)

                                                aGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           >  1:1815814/1‑48 (MQ=255)

                                                aGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           >  1:2141935/1‑48 (MQ=255)

                                                aGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           <  1:2282096/48‑1 (MQ=255)

                                                aGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGc           >  1:2550879/1‑48 (MQ=255)

                                                  gAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAAccc  <  1:29824/55‑1 (MQ=255)

                                                  gAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAAccc  <  1:2880252/55‑1 (MQ=255)

                                                   aaGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAAccc  <  1:3100925/54‑1 (MQ=255)

TTTTACATCCGATGGGCGCAGTGAATAATCATGTTCTGAAGCGCCGTCAGGAAGTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAACCC  >  minE/1246290‑1246394

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