BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,724	A→C	10.6%	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	10.6%	241.7 / 27.2	104	S510S (TCT→TCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (40/53); new base C (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 = 9.80e-01

Only 100 of 133 total aligned reads displayed.

CAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGC  >  minE/608654‑608792

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

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

  gcaACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCaga                                                                      >  1:660784/1‑69 (MQ=255)

        gACAGGCCCGCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCATCagcgaa                                                                   >  1:2593470/1‑65 (MQ=255)

        gACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCaga                                                                      >  1:1344901/1‑63 (MQ=255)

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

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

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

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

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

                aGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGt                                                      <  1:1411613/71‑1 (MQ=255)

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

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

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

                 gCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTa                                                     >  1:2492386/1‑71 (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:3749532/71‑1 (MQ=255)

                     tACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTg                                                 <  1:3114234/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:2146397/1‑71 (MQ=255)

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

                          cacGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTg                                            >  1:360706/1‑71 (MQ=255)

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

                               cATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCGATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                                       <  1:3422701/71‑1 (MQ=255)

                               cATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTa                                                     >  1:4093604/1‑57 (MQ=255)

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

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

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

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

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

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

                                 tGGTGCGGATA‑ACGTCACGGTCGTGAAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAg                                     <  1:662727/71‑1 (MQ=255)

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

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

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

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

                                     gCGGATA‑ACGTCACGGTCGTGCAGCGCCATCagaga                                                                    >  1:795318/1‑36 (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‑aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatgat                             >  1:348213/1‑69 (MQ=255)

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

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

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

                                              cGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                          >  1:413721/1‑70 (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)

                                                tCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                                                       >  1:2795471/1‑39 (MQ=255)

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

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

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

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

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

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

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

                                                  aCGGTAGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                     <  1:3826744/71‑1 (MQ=255)

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

                                                  aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgata                            >  1:766778/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:1062797/70‑1 (MQ=255)

                                                        gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                <  1:50034/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)

                                                           cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg                     >  1:1562948/1‑60 (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)

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

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

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

                                                                   tCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGt            <  1:3051200/63‑1 (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)

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

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

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

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

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

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

                                                                     agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc  >  1:935191/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:861542/45‑1 (MQ=255)

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

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

CAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGC  >  minE/608654‑608792

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