BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,323,511	G→A	53.3%	G281G (GGC→GGT)	pfkA ←	6‑phosphofructokinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,323,511	0	G	A	53.3%	3.0 / 64.5	75	G281G (GGC→GGT)	pfkA	6‑phosphofructokinase I
Reads supporting (aligned to +/- strand): ref base G (14/21); new base A (17/23); total (31/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

TCGATGATGTCGTGGTGAACCAGCTGTTCGTTCTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACG  >  minE/2323458‑2323570

tCGATGATGTCGTGGTGAACCAGCTGTTCGTTCTGGATACCTACACAACGACCGCCGTAACCTGCcagc                                              >  1:3057877/1‑69 (MQ=255)

     gatgTCGTGGTGAACCAGCTGTTCGTTCTGGATACCTACACAACGACCGCCGTAACCTGCcagcagcag                                         <  1:2846209/69‑1 (MQ=255)

     gatgTCGTGGTGAACCAGCTGTTCGTTCTGGATACCTACACAACGACCGCCGTAACCTGCcagcagcag                                         <  1:2970017/69‑1 (MQ=255)

     gatgTCGTGGTGAACCAGCTGTTCGTTCTGGATACCTACACAACGACCGCCGTAACCTGCcagcagcag                                         <  1:659957/69‑1 (MQ=255)

     gatgTCGTGGTGAACCAGCTGTTCGTTCTGGATACCTACACAACGACCGCCGTAACCTGCcagcagcag                                         <  1:787922/69‑1 (MQ=255)

             ggtgAACCAGCTGTTCGTTCTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAg                                 <  1:1858990/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:1241346/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:2615969/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:2254338/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:38808/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:2104254/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:648882/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:1992493/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:1827933/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:1708824/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:1400149/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:145281/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:1807854/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:1470022/69‑1 (MQ=255)

                gAACCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCAAGCAGCAGCTCGATAGCGt                              <  1:1526756/69‑1 (MQ=255)

                 aaCCAGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:2480164/68‑1 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:601687/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:2521671/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:1354856/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:1359702/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:2268084/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:3200115/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:362979/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:93357/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:2065230/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:1494352/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:66379/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:1543635/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:1894148/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:184856/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:1848141/1‑69 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAg                           >  1:879082/1‑67 (MQ=255)

                     aGCTGTTCGTTCTGGATACCTAC‑GAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                         >  1:416705/1‑68 (MQ=255)

                         gttcgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                     >  1:1245800/1‑69 (MQ=255)

                         gttcgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                     >  1:2036442/1‑69 (MQ=255)

                            cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATa                   <  1:2843954/68‑1 (MQ=255)

                            cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATa                   <  1:3045183/68‑1 (MQ=255)

                            cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATa                   <  1:3055639/68‑1 (MQ=255)

                            cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATa                   <  1:2018933/68‑1 (MQ=255)

                            cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATa                   <  1:2187042/68‑1 (MQ=255)

                            cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATAc                  >  1:1725362/1‑69 (MQ=255)

                            cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATAc                  >  1:2315724/1‑69 (MQ=255)

                            cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATAc                  >  1:837423/1‑69 (MQ=255)

                            cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATAc                  >  1:2225426/1‑69 (MQ=255)

                            cgttcTGGATACCTACACAACGACCGCAGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATAc                  >  1:2110762/1‑69 (MQ=255)

                                cTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGt                              <  1:2291226/53‑1 (MQ=38)

                                         tACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAAt     >  1:1897715/1‑69 (MQ=255)

                                          acacAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                     <  1:984062/52‑1 (MQ=255)

                                          acacAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                     <  1:2560595/52‑1 (MQ=255)

                                          acacAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                     <  1:589294/52‑1 (MQ=255)

                                          acacAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                     <  1:3297257/52‑1 (MQ=255)

                                          acacAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCa                     <  1:2420389/52‑1 (MQ=255)

                                           cacaACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCg                               >  1:726326/1‑41 (MQ=255)

                                           cacaACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTa                             >  1:219305/1‑43 (MQ=255)

                                           cacaACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTa                             >  1:3101123/1‑43 (MQ=255)

                                           cacaACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTa                             >  1:624244/1‑43 (MQ=255)

                                           cacaACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTa                             >  1:313893/1‑43 (MQ=255)

                                             caACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACg  <  1:433744/68‑1 (MQ=255)

                                             caACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACg  <  1:837087/68‑1 (MQ=255)

                                             caACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACg  <  1:3016590/68‑1 (MQ=255)

                                             caACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACg  <  1:267176/68‑1 (MQ=255)

                                             caACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACg  <  1:1482090/68‑1 (MQ=255)

                                             caACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACg  <  1:3258550/68‑1 (MQ=255)

                                               aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg               <  1:219928/53‑1 (MQ=39)

                                               aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg               <  1:555699/53‑1 (MQ=39)

                                               aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg               <  1:2047111/53‑1 (MQ=39)

                                               aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg               <  1:3135311/53‑1 (MQ=39)

                                               aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg               <  1:2181920/53‑1 (MQ=39)

                                               aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg               <  1:3080509/53‑1 (MQ=39)

                                               aCGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAGCGCCCATACggg               <  1:818720/53‑1 (MQ=39)

TCGATGATGTCGTGGTGAACCAGCTGTTCGTTCTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAATACG  >  minE/2323458‑2323570

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