BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,006	A→C	7.3%	A416A (GCT→GCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,006	0	A	C	7.3%	315.6 / 18.7	123	A416A (GCT→GCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (42/72); new base C (5/4); total (47/76)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 125 total aligned reads displayed.

CCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGG  >  minE/608936‑609070

ccAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                                                  <  1:2004439/71‑1 (MQ=255)

  aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                                                >  1:746943/1‑71 (MQ=255)

  aGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGc                                                                >  1:1523287/1‑71 (MQ=255)

   ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCATCAAGCAATAGCg                                                               >  1:3409611/1‑71 (MQ=255)

   ggTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgcg                                                               >  1:584689/1‑71 (MQ=255)

       tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                                           >  1:1450037/1‑71 (MQ=255)

         cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCa                                                         >  1:1887477/1‑71 (MQ=255)

         cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCa                                                         >  1:1403558/1‑71 (MQ=255)

           cgtgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTc                                                       <  1:3689070/68‑1 (MQ=255)

           cgcgCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc                                                       <  1:2905324/71‑1 (MQ=255)

           cgcgCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc                                                       <  1:3810445/71‑1 (MQ=255)

            gtgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTaa                                                            >  1:3461263/3‑64 (MQ=255)

              gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                                    <  1:2275509/71‑1 (MQ=255)

                aCCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctgt                                                  <  1:1484372/71‑1 (MQ=255)

                  cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCTTCgt                                                     <  1:203859/66‑1 (MQ=255)

                  cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                <  1:404777/71‑1 (MQ=255)

                  cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGAAAGCAATAGCGTAGTCATCGttgttg                                                <  1:3363081/71‑1 (MQ=255)

                  cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGAAAGCAATAGCGTAGTCATCGttgttg                                                <  1:76445/71‑1 (MQ=255)

                   gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                >  1:2659473/1‑70 (MQ=255)

                   gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGa                                               <  1:2835979/71‑1 (MQ=255)

                    aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                              <  1:2336403/71‑1 (MQ=255)

                    aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                                                >  1:2354153/1‑69 (MQ=255)

                     agaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                              >  1:3119929/1‑70 (MQ=255)

                      gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                  <  1:909458/65‑1 (MQ=255)

                      gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                  <  1:1567467/65‑1 (MQ=255)

                      gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                            <  1:3471068/71‑1 (MQ=255)

                      gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                            <  1:2244575/71‑1 (MQ=255)

                      gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                            <  1:2700263/71‑1 (MQ=255)

                      gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                            >  1:2804005/1‑71 (MQ=255)

                      gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                            <  1:3833286/71‑1 (MQ=255)

                      gaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt                                            <  1:3451581/71‑1 (MQ=255)

                       aaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                                               <  1:264776/51‑1 (MQ=255)

                       aaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                                               <  1:2014626/51‑1 (MQ=255)

                        aCTGCATTTGTTTAGCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                  <  1:1304360/63‑1 (MQ=255)

                        aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                  <  1:2532501/63‑1 (MQ=255)

                        aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                  <  1:1416918/63‑1 (MQ=255)

                        aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                  <  1:7253/63‑1 (MQ=255)

                        aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                  <  1:303174/63‑1 (MQ=255)

                         cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                         >  1:1524361/1‑71 (MQ=255)

                            cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCgg                                        <  1:785074/69‑1 (MQ=255)

                            cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCgg                                        <  1:908338/69‑1 (MQ=255)

                            cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGAc                                      >  1:1919654/1‑71 (MQ=255)

                              tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                                               <  1:1631342/44‑1 (MQ=255)

                              tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                <  1:1793417/59‑1 (MQ=255)

                              tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACg                                     >  1:1357/1‑70 (MQ=255)

                              tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACg                                     >  1:642692/1‑70 (MQ=255)

                               ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                       >  1:3136309/1‑67 (MQ=255)

                               ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                       >  1:2575094/1‑67 (MQ=255)

                                tGTTTACCAACGATCATCGGGCTTACGCAGCGAGCAATAGCGTa                                                             >  1:592778/1‑44 (MQ=255)

                                tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCa                                   <  1:3882253/70‑1 (MQ=255)

                                 gTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                                  >  1:1297535/1‑70 (MQ=255)

                                  tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                                                   <  1:1340539/52‑1 (MQ=255)

                                  tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                                                   <  1:144171/52‑1 (MQ=255)

                                  tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATc                                 <  1:2813371/70‑1 (MQ=255)

                                  tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATc                                 <  1:1992040/70‑1 (MQ=255)

                                    tACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                              >  1:1184445/1‑55 (MQ=255)

                                     aCCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGt                             >  1:1691041/1‑71 (MQ=255)

                                      ccAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGt                             <  1:2234179/70‑1 (MQ=255)

                                      ccAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                            <  1:2214462/71‑1 (MQ=255)

                                      ccAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                            <  1:534087/71‑1 (MQ=255)

                                       cAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                         >  1:1200391/1‑57 (MQ=255)

                                         aCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                       >  1:2558266/1‑57 (MQ=255)

                                         aCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                       >  1:2496742/1‑57 (MQ=255)

                                         aCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                          >  1:1839732/1‑70 (MQ=255)

                                          cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAgg                              <  1:596952/65‑1 (MQ=255)

                                          cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                        <  1:1121057/71‑1 (MQ=255)

                                          cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                        <  1:3198513/71‑1 (MQ=255)

                                           gATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                                           <  1:3776659/35‑1 (MQ=255)

                                           gATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt                       <  1:1449555/71‑1 (MQ=255)

                                           gATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt                       <  1:3827410/71‑1 (MQ=255)

                                             tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCa                                                         >  1:469559/1‑35 (MQ=255)

                                               atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt                                                 <  1:3935854/41‑1 (MQ=255)

                                               atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                    >  1:39150/1‑70 (MQ=255)

                                               atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                    >  1:2943661/1‑70 (MQ=255)

                                                tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                                    <  1:2213414/37‑1 (MQ=255)

                                                tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                         <  1:1884846/48‑1 (MQ=255)

                                                tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGGCATCGTTCTCa                  <  1:1346632/71‑1 (MQ=255)

                                                 cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                  >  1:850367/1‑70 (MQ=255)

                                                 cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                  >  1:3156580/1‑70 (MQ=255)

                                                  gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                         <  1:2121318/46‑1 (MQ=255)

                                                     cTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa                <  1:884056/68‑1 (MQ=255)

                                                     cTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg             <  1:3749711/71‑1 (MQ=255)

                                                     cTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTg             <  1:2620261/71‑1 (MQ=255)

                                                        aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                       >  1:1964762/1‑42 (MQ=255)

                                                         cGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCa           >  1:1262507/1‑69 (MQ=255)

                                                         cGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCa           >  1:477454/1‑69 (MQ=255)

                                                            agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag        >  1:2136143/1‑69 (MQ=255)

                                                            agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa      <  1:3280242/71‑1 (MQ=255)

                                                            agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa      <  1:2359009/71‑1 (MQ=255)

                                                             gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATAGTTCTCATACTGCAGagaag     <  1:1125767/71‑1 (MQ=255)

                                                              caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAgaga       <  1:3022563/68‑1 (MQ=255)

                                                               aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAg                               >  1:59411/1‑43 (MQ=255)

                                                               aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                    <  1:1411631/54‑1 (MQ=255)

                                                               aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                    <  1:3266642/54‑1 (MQ=255)

                                                               aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag   <  1:2791412/71‑1 (MQ=255)

                                                               aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag   <  1:325690/71‑1 (MQ=255)

                                                               aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAACGTTCTCATACTGCAGAGAagag   <  1:1999685/71‑1 (MQ=255)

                                                                agcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt                       <  1:3252434/50‑1 (MQ=255)

                                                                agcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag   <  1:438393/70‑1 (MQ=255)

                                                                  caaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga         >  1:2145448/1‑62 (MQ=255)

CCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGG  >  minE/608936‑609070

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