BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,727	G→A	25.2%	A509A (GCC→GCT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,727	0	G	A	25.2%	175.9 / 58.1	99	A509A (GCC→GCT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (31/43); new base A (13/12); total (44/55)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

CGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGC  >  minE/608660‑608796

ccgACAGGCCCGCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCt                                                                    <  1:1898037/69‑1 (MQ=255)

cGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCt                                                                    >  1:3421578/1‑71 (MQ=255)

 cgACAGGCCCGCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCtt                                                                   >  1:150386/2‑71 (MQ=255)

   aCAGGCCCGCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCg                                                                 >  1:943171/1‑71 (MQ=255)

   aCAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                                                 >  1:1107471/1‑71 (MQ=255)

   aCAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                                                 >  1:3372418/1‑71 (MQ=255)

    cAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                                                 >  1:3149032/1‑70 (MQ=255)

    cAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGt                                                                <  1:3694058/71‑1 (MQ=255)

      gACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCGGAGAGGCTTCGTa                                                               <  1:215811/70‑1 (MQ=255)

      gACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTa                                                               >  1:3193277/1‑70 (MQ=255)

       aCCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGc                                                             <  1:1685785/71‑1 (MQ=255)

           gCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGt                                                          >  1:2710663/1‑70 (MQ=255)

           gCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTa                                                         >  1:1630226/1‑71 (MQ=255)

            cGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTAc                                                        <  1:2629928/71‑1 (MQ=255)

              aTACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                                                           <  1:1217645/66‑1 (MQ=255)

              aTACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACtt                                                      <  1:1694886/71‑1 (MQ=255)

               tACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTg                                                     >  1:1711142/1‑71 (MQ=255)

                aCCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt                                                    >  1:1712925/1‑71 (MQ=255)

                aCCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt                                                    >  1:411857/1‑71 (MQ=255)

                 ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTc                                                   <  1:3614236/71‑1 (MQ=255)

                  cacaCGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCg                                                  <  1:256898/71‑1 (MQ=255)

                    catgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGt                                                 <  1:1649123/67‑1 (MQ=255)

                    cacGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt                                                    >  1:1554954/1‑67 (MQ=255)

                      tgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCa                                              >  1:2448498/2‑71 (MQ=255)

                      tgCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCa                                              >  1:514557/2‑71 (MQ=255)

                      cGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTa                                                         <  1:1576949/60‑1 (MQ=255)

                          aTGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATg                                            >  1:3881228/1‑69 (MQ=255)

                          aTGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTa                                          <  1:448813/71‑1 (MQ=255)

                          aTGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTa                                          <  1:1979514/71‑1 (MQ=255)

                             gTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGt                                        <  1:1766988/70‑1 (MQ=255)

                             gTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTg                                       <  1:2598374/71‑1 (MQ=255)

                              tGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTgg                                      <  1:3700027/71‑1 (MQ=255)

                              tGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTgg                                      <  1:693798/71‑1 (MQ=255)

                                cGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTc                                                                  >  1:617567/1‑41 (MQ=255)

                                cGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCTTGCATGTAGTGga                                     <  1:180471/70‑1 (MQ=255)

                                cGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                     <  1:115108/70‑1 (MQ=255)

                                cGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                     <  1:3045120/70‑1 (MQ=255)

                                 ggATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatg                                   <  1:1051507/71‑1 (MQ=255)

                                 ggATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgat                                    >  1:1712621/1‑70 (MQ=255)

                                 ggATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgat                                    >  1:2498510/1‑70 (MQ=255)

                                  gATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                  <  1:1510895/71‑1 (MQ=255)

                                    tA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGatgatg                                <  1:2111386/71‑1 (MQ=255)

                                     a‑aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                               >  1:2848019/1‑71 (MQ=255)

                                     a‑aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                               >  1:1109247/1‑71 (MQ=255)

                                       acatCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt                              <  1:2451649/68‑1 (MQ=255)

                                       aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTc                                                                  <  1:2784786/35‑1 (MQ=255)

                                       aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                               <  1:284503/70‑1 (MQ=255)

                                         gTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTca                            <  1:3765323/71‑1 (MQ=255)

                                         atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt                              >  1:1375221/2‑69 (MQ=255)

                                            aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                                           <  1:2153764/53‑1 (MQ=255)

                                            aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                                           <  1:3940187/53‑1 (MQ=255)

                                            aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                         >  1:2448796/1‑71 (MQ=255)

                                            aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                         >  1:3690028/1‑71 (MQ=255)

                                            aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                         <  1:3659170/71‑1 (MQ=255)

                                            aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                         <  1:583284/71‑1 (MQ=255)

                                            aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagt                          >  1:1746528/1‑70 (MQ=255)

                                             cGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCg                                                  <  1:2605998/45‑1 (MQ=255)

                                             cGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgc                        <  1:1547791/71‑1 (MQ=255)

                                              ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                              <  1:923735/64‑1 (MQ=255)

                                              ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                              <  1:1271516/64‑1 (MQ=255)

                                              ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                         >  1:3745091/1‑69 (MQ=255)

                                              ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTc                             >  1:2678580/1‑65 (MQ=255)

                                               gTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGc                                               <  1:2052728/46‑1 (MQ=255)

                                               gTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTca                            <  1:827209/65‑1 (MQ=255)

                                               atCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgc                      <  1:1986661/70‑3 (MQ=255)

                                                tCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagt                     >  1:1361555/1‑71 (MQ=255)

                                                 cGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgct                     >  1:3892463/1‑67 (MQ=255)

                                                 cGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgct                     >  1:3818077/1‑67 (MQ=255)

                                                  gTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcag                           <  1:3051711/63‑1 (MQ=255)

                                                  gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                    <  1:3364346/70‑1 (MQ=255)

                                                  gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                   >  1:3172857/1‑71 (MQ=255)

                                                  gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                   >  1:524081/1‑71 (MQ=255)

                                                   tGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                    <  1:1989274/69‑1 (MQ=255)

                                                    gCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGAt                  <  1:1474046/70‑1 (MQ=255)

                                                     cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgat                                 >  1:3949500/1‑54 (MQ=255)

                                                     cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                 >  1:1270476/1‑70 (MQ=255)

                                                     cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGt                <  1:1276183/71‑1 (MQ=255)

                                                      aGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGt                                        >  1:599214/1‑46 (MQ=255)

                                                      aGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTa               >  1:1871441/1‑71 (MQ=255)

                                                        cgcCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACt             >  1:1974860/1‑71 (MQ=255)

                                                          ccATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGt           <  1:2932575/71‑1 (MQ=255)

                                                           cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt          <  1:3283633/71‑1 (MQ=255)

                                                           cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt          <  1:1634670/71‑1 (MQ=255)

                                                            aTCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg                         <  1:2296146/55‑3 (MQ=255)

                                                            aTCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTgcac         >  1:1489148/1‑68 (MQ=255)

                                                             tCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTAc              <  1:1963599/65‑1 (MQ=255)

                                                             tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTgg                                      <  1:7523/41‑1 (MQ=255)

                                                             tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt          >  1:233172/1‑69 (MQ=255)

                                                              cAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                     <  1:3129243/41‑1 (MQ=255)

                                                              cAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTag       >  1:258021/1‑71 (MQ=255)

                                                               agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                           >  1:3525398/1‑50 (MQ=255)

                                                               agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                           >  1:925830/1‑50 (MQ=255)

                                                               agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc      >  1:638018/1‑71 (MQ=255)

                                                                gagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                  <  1:2817946/42‑1 (MQ=255)

                                                                 agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                               >  1:213058/1‑44 (MQ=255)

                                                                 agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                               >  1:3871573/1‑44 (MQ=255)

                                                                 agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                               >  1:3791007/1‑44 (MQ=255)

                                                                  gagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGatgatg                                <  1:326538/42‑1 (MQ=255)

                                                                  gagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccag   >  1:3099829/1‑71 (MQ=255)

                                                                   agGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc  >  1:3524790/1‑71 (MQ=255)

CGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGC  >  minE/608660‑608796

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