BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,727	G→A	10.6%	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	10.6%	309.4 / 20.6	104	A509A (GCC→GCT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (40/53); new base A (7/4); total (47/57)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

Only 100 of 131 total aligned reads displayed.

ACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCC  >  minE/608659‑608797

aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagagag                                                                        >  1:3012353/1‑69 (MQ=255)

aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGc                                                                      >  1:1631695/1‑71 (MQ=255)

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

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

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

      aGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGt                                                                 <  1:3144915/70‑1 (MQ=255)

         ccAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCt                                                             >  1:2195975/1‑71 (MQ=255)

          ccgCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTg                                                            >  1:3833725/3‑71 (MQ=255)

          cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                                                            >  1:2141826/1‑71 (MQ=255)

            gCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTa                                                          >  1:3468631/1‑71 (MQ=255)

            gCGATACCACACGCCATGTTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTa                                                          <  1:3979132/71‑1 (MQ=255)

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

            gCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTa                                                          <  1:1252080/71‑1 (MQ=255)

             cGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTAc                                                         <  1:1141314/71‑1 (MQ=255)

               aTACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGt                                                                 <  1:509686/61‑1 (MQ=255)

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

                tACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTg                                                      <  1:3114234/71‑1 (MQ=255)

                tACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTg                                                      <  1:1738228/71‑1 (MQ=255)

                tACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTg                                                      <  1:3749532/71‑1 (MQ=255)

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

                  ccACCCGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTc                                                    >  1:1835232/1‑71 (MQ=255)

                  ccACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTc                                                    >  1:805204/1‑71 (MQ=255)

                    acacGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGt                                                  >  1:2724808/1‑71 (MQ=255)

                    acacGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGt                                                  >  1:2146397/1‑71 (MQ=255)

                     cacGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTg                                                 <  1:2036007/71‑1 (MQ=255)

                         ccATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATg                                             <  1:2932686/71‑1 (MQ=255)

                         ccATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATg                                             <  1:3656900/71‑1 (MQ=255)

                          cATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                                            <  1:2955700/71‑1 (MQ=255)

                          cATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                                            <  1:885448/71‑1 (MQ=255)

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

                            tGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAg                                          <  1:3583059/71‑1 (MQ=255)

                              gTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt                                                     <  1:3623476/58‑1 (MQ=255)

                              gTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTa                                           >  1:1338850/1‑68 (MQ=255)

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

                              gTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTg                                        >  1:2903820/1‑71 (MQ=255)

                              gTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTg                                        >  1:3553017/1‑71 (MQ=255)

                                  ggATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTg                                                      <  1:2604745/53‑1 (MQ=255)

                                    aTA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatgat                                  <  1:3531184/71‑1 (MQ=255)

                                     tA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCg                                                   >  1:1791894/1‑53 (MQ=255)

                                     tA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTa                                           >  1:3804094/1‑61 (MQ=255)

                                      a‑aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCa                                               >  1:1081561/1‑56 (MQ=255)

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

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

                                        acatCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatg                                    <  1:2497046/63‑1 (MQ=255)

                                         cGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTc                              <  1:1147769/71‑1 (MQ=255)

                                         cGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTc                              <  1:837458/71‑1 (MQ=255)

                                          gTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTgg                                       >  1:1133633/1‑61 (MQ=255)

                                          gTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                   >  1:3446339/1‑65 (MQ=255)

                                           tCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                            >  1:2991471/1‑71 (MQ=255)

                                           tCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                            <  1:4123714/71‑1 (MQ=255)

                                           tCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGc                                                <  1:2406634/51‑1 (MQ=255)

                                            cACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                            <  1:811859/70‑1 (MQ=255)

                                             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatg                                    <  1:308759/61‑1 (MQ=255)

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

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

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

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

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

                                             aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgata                                 >  1:766778/1‑63 (MQ=255)

                                             aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgata                                 >  1:1691388/1‑63 (MQ=255)

                                             aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg                          >  1:1016592/1‑69 (MQ=255)

                                                 tCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgct                      <  1:4139864/71‑4 (MQ=255)

                                                 tCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagt                      <  1:3079021/71‑1 (MQ=255)

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

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

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

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

                                                    tGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                                            <  1:2890810/46‑1 (MQ=255)

                                                      cAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                  <  1:218100/70‑1 (MQ=255)

                                                         cgcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACt              <  1:2627648/71‑1 (MQ=255)

                                                          gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGatgatg                                 >  1:1242730/1‑51 (MQ=255)

                                                          gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgca                        >  1:320981/1‑60 (MQ=255)

                                                          gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgca                        >  1:1868695/1‑60 (MQ=255)

                                                           ccATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAg                                          >  1:772579/1‑41 (MQ=255)

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

                                                             aTCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                                            <  1:2714742/37‑1 (MQ=255)

                                                             aTCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt           >  1:1522671/1‑70 (MQ=255)

                                                              tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                               <  1:4145358/49‑1 (MQ=255)

                                                              tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTa         <  1:2708923/71‑1 (MQ=255)

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

                                                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTa                                           >  1:802421/1‑35 (MQ=255)

                                                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAg                                          >  1:530679/1‑36 (MQ=255)

                                                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                     <  1:3079110/57‑1 (MQ=255)

                                                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                    >  1:673111/1‑58 (MQ=255)

                                                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGAt                   >  1:1827846/1‑59 (MQ=255)

                                                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                  <  1:1085544/60‑1 (MQ=255)

                                                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGttt          >  1:3886515/1‑68 (MQ=255)

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

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

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

                                                                agagagGCTTCGTAGCTGTACTTGTCGTGAATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc       >  1:2619530/1‑71 (MQ=255)

                                                                 gagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                            <  1:3372781/49‑1 (MQ=255)

                                                                  agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTgg                                       <  1:85475/37‑1 (MQ=255)

                                                                  agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                               <  1:2126694/45‑1 (MQ=255)

                                                                  agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                               <  1:1384714/45‑1 (MQ=255)

                                                                  agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                               <  1:861542/45‑1 (MQ=255)

                                                                  agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcag                       <  1:184330/53‑1 (MQ=255)

                                                                  agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTa                >  1:2418589/1‑60 (MQ=255)

                                                                   gagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcca     >  1:4136021/1‑70 (MQ=255)

                                                                    agGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTca                             >  1:4096061/1‑45 (MQ=255)

ACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCC  >  minE/608659‑608797

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