BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,557,643	A→C	14.3%	G99G (GGA→GGC)	secB →	protein export chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,557,643	0	A	C	14.3%	222.2 / 28.7	112	G99G (GGA→GGC)	secB	protein export chaperone
Reads supporting (aligned to +/- strand): ref base A (66/30); new base C (5/11); total (71/41)
Fisher's exact test for biased strand distribution p-value = 9.52e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 115 total aligned reads displayed.

AAGTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATCCCGCGGTA  >  minE/2557576‑2557710

aaGTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAg                                                                    <  1:1394194/69‑1 (MQ=255)

  gTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                                  <  1:3278901/69‑1 (MQ=255)

  gTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                                  <  1:1817181/69‑1 (MQ=255)

  gTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                                  <  1:133077/69‑1 (MQ=255)

  gTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                                  <  1:487776/69‑1 (MQ=255)

   ttCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                                  <  1:2270888/68‑1 (MQ=255)

                  ttttCTCCATCGCGGGTATCGAAGGCACCGAGATGGCGCATTGCCTGGGAGCATACTg                                                             >  1:2966808/1‑58 (MQ=255)

                  ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                             >  1:442836/1‑58 (MQ=255)

                  ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                             >  1:1116192/1‑58 (MQ=255)

                  ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                             >  1:2657825/1‑58 (MQ=255)

                  ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                             >  1:1559660/1‑58 (MQ=255)

                  ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                             >  1:152720/1‑58 (MQ=255)

                  ttttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTg                                                             >  1:1521224/1‑58 (MQ=255)

                    ttCTCCATCGCGGGTATCTAAGGCACCCAGATGGCGCATTGCCTGGGAg                                                                    <  1:2499124/49‑1 (MQ=255)

                    ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTTGGAg                                                                    <  1:3070856/49‑1 (MQ=255)

                    ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                                  >  1:2697246/1‑51 (MQ=255)

                    ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                                  >  1:927479/1‑51 (MQ=255)

                    ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCa                                                                  >  1:1661213/1‑51 (MQ=255)

                    ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTg                                                >  1:3211663/1‑69 (MQ=255)

                    ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTg                                                >  1:2217764/1‑69 (MQ=255)

                    ttCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTg                                                >  1:2288861/1‑69 (MQ=255)

                            cgcgGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTa                                        >  1:3275528/1‑69 (MQ=255)

                            cgcgGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTa                                        >  1:3301975/1‑69 (MQ=255)

                            cgcgGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTa                                        >  1:1950354/1‑69 (MQ=255)

                            cgcgGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTa                                        >  1:1608834/1‑69 (MQ=255)

                            cgcgGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTa                                        >  1:872905/1‑69 (MQ=255)

                               gggTATCGAAGGCACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGc                                     >  1:1608694/1‑69 (MQ=255)

                               gggTATCGAAGGCACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGc                                     >  1:1065191/1‑69 (MQ=255)

                               gggTATCGAAGGCACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGc                                     >  1:2768343/1‑69 (MQ=255)

                               gggTATCGAAGGCACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGc                                     >  1:2733894/1‑69 (MQ=255)

                                  tATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTc                                   >  1:2725410/1‑68 (MQ=255)

                                    tCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTg                                <  1:2038342/69‑1 (MQ=255)

                                    tCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTg                                <  1:2057310/69‑1 (MQ=255)

                                    tCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTg                                <  1:771044/69‑1 (MQ=255)

                                    tCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCCGTTCCCGTATGCTCGTg                                <  1:1057837/69‑1 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACa                                                     <  1:726917/46‑1 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTc                                             >  1:1656154/1‑54 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                                          <  1:254193/57‑1 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                                          <  1:1972026/57‑1 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                                          <  1:2066819/57‑1 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                                          <  1:896017/57‑1 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                                          <  1:1333543/57‑1 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                                          <  1:862784/57‑1 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCg                                          <  1:1649430/57‑1 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                              <  1:2170587/69‑1 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                              <  1:1365716/69‑1 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                              <  1:668640/69‑1 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                              >  1:1652118/1‑69 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                              >  1:2236176/1‑69 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                              <  1:222134/69‑1 (MQ=255)

                                      gAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAg                              <  1:2020557/69‑1 (MQ=255)

                                           cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt                         <  1:1048471/69‑1 (MQ=255)

                                           cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt                         <  1:2976618/69‑1 (MQ=255)

                                           cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt                         <  1:3009568/69‑1 (MQ=255)

                                           cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt                         <  1:826880/69‑1 (MQ=255)

                                           cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt                         <  1:3049050/69‑1 (MQ=255)

                                           cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt                         <  1:478169/69‑1 (MQ=255)

                                           cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt                         <  1:97286/69‑1 (MQ=255)

                                           cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt                         <  1:2050772/69‑1 (MQ=255)

                                           cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt                         <  1:2461186/69‑1 (MQ=255)

                                           cACCCAGATGGCACATTGCCTGGGCGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCAt                         <  1:1516002/69‑1 (MQ=255)

                                             cccAGATGGCGGATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATc                        <  1:1990563/68‑1 (MQ=255)

                                                    ggCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCAc                      <  1:2427240/63‑1 (MQ=255)

                                                    ggCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCAc                      <  1:424132/63‑1 (MQ=255)

                                                    ggCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCAt                >  1:3324936/1‑69 (MQ=255)

                                                    ggCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCAt                >  1:1361268/1‑69 (MQ=255)

                                                    ggCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCAt                >  1:2494988/1‑69 (MQ=255)

                                                    ggCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCCCGTGAGTGCATCACCAGCAt                >  1:2985941/1‑69 (MQ=255)

                                                       gcATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGt             >  1:2135788/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATg                                      >  1:2487907/1‑39 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGc                           >  1:874999/1‑50 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATgg              >  1:352482/1‑63 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATgg              >  1:768743/1‑63 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATgg              >  1:2074026/1‑63 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATgg              >  1:1557187/1‑63 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATgg              >  1:1491596/1‑63 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATgg              >  1:2239907/1‑63 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATgg              >  1:1453348/1‑63 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATgg              >  1:2292757/1‑63 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATgg              >  1:2862712/1‑63 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATgg              >  1:837081/1‑63 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATgg              >  1:2337119/1‑63 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATccc        >  1:771217/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATccc        >  1:782464/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATccc        >  1:1365977/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATccc        >  1:3377685/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATccc        >  1:1569907/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATccc        >  1:422688/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATccc        >  1:2096983/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATccc        >  1:221990/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATccc        >  1:2291866/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATccc        >  1:2335100/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATccc        >  1:2385981/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATccc        >  1:2594608/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTAGTGAGTGCATCACCAGCATGGTATccc        >  1:1228137/1‑69 (MQ=255)

                                                            gCCTGGGAGCATACTGCACGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATgg              >  1:2202806/1‑63 (MQ=255)

                                                            gCCTGGGAGCATACTCCCCGAACATTCTGTTCCCGTATGCt                                    >  1:1081629/1‑41 (MQ=255)

                                                                   aGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATCCCGCGGTa  <  1:1981733/68‑1 (MQ=255)

                                                                   aGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATCCCGCGGTa  <  1:2481231/68‑1 (MQ=255)

                                                                   aGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATCCCGCGGTa  <  1:1004220/68‑1 (MQ=255)

AAGTTCAGCAGGGCGGTATTTTCTCCATCGCGGGTATCGAAGGCACCCAGATGGCGCATTGCCTGGGAGCATACTGCCCGAACATTCTGTTCCCGTATGCTCGTGAGTGCATCACCAGCATGGTATCCCGCGGTA  >  minE/2557576‑2557710

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