BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,295	A→C	9.3%	A653A (GCT→GCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,295	0	A	C	9.3%	151.3 / 13.5	66	A653A (GCT→GCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (30/29); new base C (3/3); total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

CAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCC  >  minE/608229‑608363

cAGTGCGTTCGGAACGATGGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGc                                                                    <  1:2831040/69‑1 (MQ=255)

   tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                               <  1:542023/71‑1 (MQ=255)

   tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                               <  1:491343/71‑1 (MQ=255)

     cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                               >  1:2842881/1‑69 (MQ=255)

         cGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTc                                                          <  1:778552/70‑1 (MQ=255)

          ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                        <  1:787453/71‑1 (MQ=255)

          ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                        <  1:1559083/71‑1 (MQ=255)

          ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                                        <  1:3146303/71‑1 (MQ=255)

          ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                         <  1:1496844/70‑1 (MQ=255)

           gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                       <  1:2443807/71‑1 (MQ=255)

           gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                       <  1:3485917/71‑1 (MQ=255)

             aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                                       >  1:2310104/1‑67 (MQ=255)

             aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                                     <  1:3595411/71‑1 (MQ=255)

             aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                                     <  1:1848774/71‑1 (MQ=255)

              cGATAGAGAAGGTGTAGGAGATACCCTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagagag                                                    >  1:3732737/1‑71 (MQ=255)

              cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagagag                                                    >  1:2751276/1‑71 (MQ=255)

                aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                       >  1:3789780/1‑66 (MQ=255)

                aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                                  >  1:1160041/1‑71 (MQ=255)

                     gaAGGAGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                                >  1:116959/1‑68 (MQ=255)

                      aaGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGc                                            <  1:2397246/71‑1 (MQ=255)

                       aGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCa                                           >  1:66239/1‑71 (MQ=255)

                          tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                               <  1:1991040/48‑1 (MQ=255)

                            tAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCttt                                      >  1:1331584/1‑71 (MQ=255)

                              ggAGATACAATCTTTAGCGTAAGCAAACGGCAGTTTAGc                                                                    >  1:702882/1‑39 (MQ=255)

                               gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTTGCAACGGa                                                              <  1:2092376/44‑1 (MQ=255)

                                 gaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTc                                                          <  1:3062461/46‑1 (MQ=255)

                                  aTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGc                                                                    <  1:417379/35‑1 (MQ=255)

                                   tACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                                       <  1:3026541/47‑1 (MQ=255)

                                   tACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                 <  1:433874/69‑1 (MQ=255)

                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                         <  1:1565750/44‑1 (MQ=255)

                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                                         <  1:2312224/44‑1 (MQ=255)

                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                              >  1:1997276/1‑71 (MQ=255)

                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                              >  1:3711575/1‑71 (MQ=255)

                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                              >  1:2276328/1‑71 (MQ=255)

                                       aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagag                                                      <  1:3950012/44‑1 (MQ=255)

                                        tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCt                                                 <  1:612356/48‑1 (MQ=255)

                                        tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                    <  1:931244/61‑1 (MQ=255)

                                        tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                    <  1:57827/61‑1 (MQ=255)

                                        tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                  <  1:2276715/63‑1 (MQ=255)

                                        tCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                          >  1:2824860/1‑71 (MQ=255)

                                          tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                        >  1:86378/1‑71 (MQ=255)

                                            tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                       >  1:2274216/1‑70 (MQ=255)

                                            tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                       >  1:2871084/1‑70 (MQ=255)

                                            tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                       >  1:162429/1‑70 (MQ=255)

                                            tAGCGTAAGCAAACCGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                       >  1:3170732/1‑70 (MQ=255)

                                              gCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                             <  1:2989302/46‑1 (MQ=255)

                                                gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                  <  1:3083121/71‑1 (MQ=255)

                                                gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                  >  1:3206440/1‑71 (MQ=255)

                                                 tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                  <  1:1532937/54‑1 (MQ=255)

                                                   aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt               >  1:3070467/1‑71 (MQ=255)

                                                   aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt               >  1:1218867/1‑71 (MQ=255)

                                                   aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt               >  1:1485182/1‑71 (MQ=255)

                                                       aaCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt              <  1:51290/68‑1 (MQ=255)

                                                        aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                                >  1:1588358/1‑49 (MQ=255)

                                                         cGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAc         >  1:2351881/1‑71 (MQ=255)

                                                          ggCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAcc        <  1:1757453/71‑1 (MQ=255)

                                                            cAGTTTCGCCACCGAGGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGCCCg      >  1:244307/1‑71 (MQ=255)

                                                              gTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGCCCggc    >  1:442481/1‑70 (MQ=255)

                                                               tttCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAg            <  1:3423367/62‑1 (MQ=255)

                                                               tttCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAg            <  1:3418642/62‑1 (MQ=255)

                                                                 tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg                 >  1:2472862/1‑55 (MQ=255)

                                                                 tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc  >  1:249141/1‑70 (MQ=255)

                                                                 tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc  >  1:2686875/1‑70 (MQ=255)

                                                                 tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc  >  1:1308062/1‑70 (MQ=255)

                                                                 tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc  >  1:2777947/1‑70 (MQ=255)

                                                                 tAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTcc  >  1:816383/1‑70 (MQ=255)

                                                                  aGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt              <  1:3357585/57‑1 (MQ=255)

CAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCC  >  minE/608229‑608363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`