BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,706	G→A	8.7%	D516D (GAC→GAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,706	0	G	A	8.7%	359.1 / 17.4	115	D516D (GAC→GAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (53/52); new base A (8/2); total (61/54)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 116 total aligned reads displayed.

GCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTG  >  minE/608638‑608773

gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGt                                                                    <  1:1035655/70‑1 (MQ=255)

 cAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCg                                                                  <  1:3133595/71‑1 (MQ=255)

    aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTg                                                                >  1:1797318/1‑70 (MQ=255)

     cAGGGAGTCATCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCa                                                              <  1:3074958/71‑1 (MQ=255)

      aGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAg                                                             >  1:1548224/1‑71 (MQ=255)

       gggAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAgc                                                            <  1:1725805/71‑1 (MQ=255)

          aGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGcc                                                         <  1:136066/71‑1 (MQ=255)

            tCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                                       >  1:2157333/1‑71 (MQ=255)

            tCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTAGTGCAGCGcc                                                         >  1:2489165/1‑69 (MQ=255)

             cagcagCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                                       >  1:2482461/1‑70 (MQ=255)

              agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                                       >  1:3023697/1‑69 (MQ=255)

              agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                                       >  1:1006879/1‑69 (MQ=255)

              agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                                       >  1:1709765/1‑69 (MQ=255)

              agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                                       >  1:295821/1‑69 (MQ=255)

              agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                                       >  1:1113619/1‑69 (MQ=255)

              agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                                       >  1:1030582/1‑69 (MQ=255)

              agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCAt                                                       >  1:2085778/1‑69 (MQ=255)

              agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATc                                                      <  1:1893571/70‑1 (MQ=255)

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

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

                cagcaACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCa                                                     >  1:287102/1‑69 (MQ=255)

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

                  gcaACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCag                                                    <  1:1617229/68‑1 (MQ=255)

                  gcaACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCag                                                    >  1:3579619/1‑68 (MQ=255)

                     aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGt                                                                    >  1:2999098/1‑49 (MQ=255)

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

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

                        gACAGGCCCGCGATGCCGCATGCCATAGTGC‑GATAGACATCACGATCGTGCAGCGCCAt                                                       >  1:3973355/1‑59 (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)

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

                             aCCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCag                                                    <  1:508078/57‑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)

                               cAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                                    <  1:868672/71‑1 (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‑ACGTCACGGTCGTGCAGCGCCATCag                                                    <  1:3197811/53‑1 (MQ=255)

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

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

                                  cGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTa                                  <  1:2145720/70‑1 (MQ=255)

                                    aTACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGt                                         <  1:509686/61‑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:2724808/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)

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

                                              ccATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATg                     <  1:3656900/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:3116043/71‑1 (MQ=255)

                                               cATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt                    <  1:182571/71‑1 (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:3553017/1‑71 (MQ=255)

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

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

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

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

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

                                                           a‑aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatgat          >  1:348213/1‑69 (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:837458/71‑1 (MQ=255)

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

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

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

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

                                                               atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTAc                                 >  1:122766/2‑43 (MQ=37)

                                                                tCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                                    >  1:2795471/1‑39 (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)

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

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

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

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

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

                                                                  aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagt   >  1:421674/1‑70 (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)

                                                                   cGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACt                                >  1:81826/1‑40 (MQ=255)

                                                                   cGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag    <  1:2037560/68‑1 (MQ=255)

                                                                    ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatgat          >  1:2564569/1‑61 (MQ=255)

GCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATA‑ACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTG  >  minE/608638‑608773

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