BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,323,502	A→G	60.8%	C284C (TGT→TGC)	pfkA ←	6‑phosphofructokinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,323,502	0	A	G	60.8%	24.6 / 32.6	54	C284C (TGT→TGC)	pfkA	6‑phosphofructokinase I
Reads supporting (aligned to +/- strand): ref base A (11/10); new base G (16/16); total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CGCTTCATGTTTTCGATAGCGTCGATGATGTCGTGGTGAACCAGCTGTTCGTTCTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAAT  >  minE/2323437‑2323567

cGCTTCATGTTTTCGATAGCGTCGATGATGTCGTGGTGAACCAGCTGTTCGTTCTGGATACCTAcaca                                                                 >  1:1996240/1‑68 (MQ=255)

     cATGTTTTCGATAGCGTCGATGATGTCGTGGTGAACCAGCTGTTCGTTCTGGATACCTACACAACGAcc                                                           >  1:1781256/1‑69 (MQ=255)

     cATGTTTTCGATAGCGTCGATGATGTCGTGGTGAACCAGCTGTTCGTTCTGGATACCTACACAACGAcc                                                           >  1:501181/1‑69 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                          aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                      >  1:93357/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)

                                          aGCTGTTCGTTCTGGATACCTACGCAACGACCACCGTAACCTGCCAGCAGCAGCTCGATAGCGTAAgcg                      >  1:1494352/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:3045183/68‑1 (MQ=255)

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

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

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

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

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

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

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

                                                 cgttcTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATAc               >  1:1725362/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:589294/52‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

CGCTTCATGTTTTCGATAGCGTCGATGATGTCGTGGTGAACCAGCTGTTCGTTCTGGATACCTACACAACGACCGCCGTAACCTGCCAGCAGCAGATCGATAGCGTAAGCGCCCATACGGGAAGCCAGAAT  >  minE/2323437‑2323567

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