BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,724	A→C	25.2%	S510S (TCT→TCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,724	0	A	C	25.2%	127.8 / 70.0	99	S510S (TCT→TCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (31/43); new base C (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 109 total aligned reads displayed.

CAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCC  >  minE/608654‑608793

cagcaACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCaga                                                                       <  1:1218174/71‑1 (MQ=255)

 agcaACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagag                                                                      <  1:3678985/71‑1 (MQ=255)

      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:3799173/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:194067/1‑71 (MQ=255)

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

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

                   gATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTAc                                                     <  1:3079912/70‑1 (MQ=255)

                    aTACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                                                        <  1:1217645/66‑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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                           a‑aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                            >  1:2848019/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)

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

                                                 cACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagc                       >  1:1683013/1‑70 (MQ=255)

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

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

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

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

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

                                                  aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                      <  1:2616968/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:2285052/70‑3 (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)

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

                                                         tGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                 <  1:1989274/69‑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:3439651/71‑1 (MQ=255)

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

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

                                                                 cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt       <  1:3283633/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)

                                                                      gagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcc  <  1:1464048/71‑1 (MQ=255)

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

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

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

CAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCC  >  minE/608654‑608793

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