BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	329,175	A→T	14.3%	S194S (TCA→TCT)	gsk →	inosine/guanosine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	329,175	0	A	T	14.3%	98.0 / 9.1	49	S194S (TCA→TCT)	gsk	inosine/guanosine kinase
Reads supporting (aligned to +/- strand): ref base A (16/26); new base T (4/3); total (20/29)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGA  >  minE/329114‑329242

tGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTgg                                                              >  1:2623770/1‑69 (MQ=255)

tGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTgg                                                              >  1:984755/1‑69 (MQ=255)

tGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTgg                                                              >  1:640762/1‑69 (MQ=255)

tGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTgg                                                              >  1:3128739/1‑69 (MQ=255)

tGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTgg                                                              >  1:3057566/1‑69 (MQ=255)

tGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTgg                                                              >  1:2014629/1‑69 (MQ=255)

tGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTgg                                                              >  1:1885208/1‑69 (MQ=255)

              atgccGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCg                                                <  1:837877/66‑1 (MQ=255)

              aTTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCg                                                <  1:1976869/69‑1 (MQ=255)

                      aGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAAcc                                        <  1:1607157/69‑1 (MQ=255)

                      aGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAAcc                                        <  1:2866783/69‑1 (MQ=255)

                      aGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAAcc                                        <  1:130458/69‑1 (MQ=255)

                      aGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATAACTGGTGCGTTGCAAGCCGGGTGAAcc                                        <  1:1663244/69‑1 (MQ=255)

                         tgtgATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCAt                                     >  1:2954205/1‑69 (MQ=255)

                         tgtgATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCAt                                     >  1:2703246/1‑69 (MQ=255)

                               tGCCGGAGCCTCGGCTCTGGTTCTCACCTCTTATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGa                               >  1:1556324/1‑69 (MQ=255)

                               tGCCGGAGCCTCGGCTCTGGTTCTCACCTCTTATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGa                               >  1:1352509/1‑69 (MQ=255)

                               tGCCGGAGCCTCGGCTCTGGTTCTCACCTCTTATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGa                               >  1:1293043/1‑69 (MQ=255)

                               tGCCGGAGCCTCGGCTCTGGTTCTCACCTCTTATCTGGTGCGTTGCAAGACGGGTGAACCCATGCCgg                                >  1:1191902/1‑68 (MQ=255)

                                gCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGaa                                          <  1:483971/57‑1 (MQ=255)

                                gCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGaa                                          <  1:1375283/57‑1 (MQ=255)

                                gCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGaa                                          <  1:2478129/57‑1 (MQ=255)

                                gCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGaa                                          <  1:1928918/57‑1 (MQ=255)

                                     aGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAAc                         <  1:3220471/69‑1 (MQ=255)

                                     aGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAAc                         <  1:2222659/69‑1 (MQ=255)

                                     aGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAAc                         <  1:2879509/69‑1 (MQ=255)

                                     aGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGGTGCAAGCCGGGTGAACCCATGCCGGAAGCAAc                         <  1:475744/69‑1 (MQ=255)

                                      gCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAAc                         <  1:3216543/68‑1 (MQ=255)

                                                tGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCa              <  1:3128875/69‑1 (MQ=255)

                                                tGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCa              <  1:587442/69‑1 (MQ=255)

                                                tGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCa              <  1:775805/69‑1 (MQ=255)

                                                 ggTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGcc               <  1:3158647/67‑1 (MQ=255)

                                                 ggTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGcc               <  1:2742421/67‑1 (MQ=255)

                                                 ggTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGcc               <  1:1829558/67‑1 (MQ=255)

                                                 ggTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGcc               <  1:1137674/67‑1 (MQ=255)

                                                 ggTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGcc               <  1:2525434/67‑1 (MQ=255)

                                                 ggTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGcc               <  1:1953806/67‑1 (MQ=255)

                                                 ggTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGcc               <  1:2016595/67‑1 (MQ=255)

                                                 ggTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCa              >  1:1496305/1‑68 (MQ=255)

                                                 ggTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCa              >  1:3052601/1‑68 (MQ=255)

                                                 ggTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCa              >  1:127758/1‑68 (MQ=255)

                                                 ggTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCa              >  1:257287/1‑68 (MQ=255)

                                                 ggTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCa              >  1:805391/1‑68 (MQ=255)

                                                 ggTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCa              >  1:671400/1‑68 (MQ=255)

                                                     ctcACCTCTTATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAg         <  1:2274518/69‑1 (MQ=255)

                                                     ctcACCTCTTATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAg         <  1:2722467/69‑1 (MQ=255)

                                                     ctcACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAAc                         <  1:1262650/53‑1 (MQ=255)

                                                     ctcACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAg         >  1:541289/1‑69 (MQ=255)

                                                      tcACCTCTTATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAg         <  1:229190/68‑1 (MQ=255)

                                                             atatCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCga  >  1:2276987/1‑68 (MQ=255)

                                                             atatCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCga  >  1:2634812/1‑68 (MQ=255)

                                                             atatCTGGTGCGTTGCAAGCAGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCga  >  1:3021229/1‑68 (MQ=255)

                                                             atatCTCGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCga  >  1:507/1‑68 (MQ=255)

TGCGGGCTGAAAGCATTCCGGAAGATGTGATTGCCGGAGCCTCGGCACTGGTTCTCACCTCATATCTGGTGCGTTGCAAGCCGGGTGAACCCATGCCGGAAGCAACCATGAAAGCCATTGAGTACGCGA  >  minE/329114‑329242

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`