BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,042	G→A	16.6%	D404D (GAC→GAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,042	0	G	A	16.6%	308.1 / 44.4	133	D404D (GAC→GAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (43/67); new base A (12/10); total (55/78)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 135 total aligned reads displayed.

ACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTC  >  minE/608973‑609108

aCCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAgg                                                                    <  1:3253912/70‑1 (MQ=255)

 ccAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGt                                                                   <  1:1329495/70‑1 (MQ=255)

     cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                              >  1:822877/1‑71 (MQ=255)

     cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                              <  1:4051844/71‑1 (MQ=255)

      gATCATCGGTCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt                                                             >  1:3346541/1‑71 (MQ=255)

       accattGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtt                                                            <  1:2958484/69‑1 (MQ=255)

        tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                                  >  1:947131/1‑64 (MQ=255)

        tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                                  >  1:1007848/1‑64 (MQ=255)

        tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                                  >  1:4006540/1‑64 (MQ=255)

        tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                                  >  1:1557215/1‑64 (MQ=255)

        tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                                                           <  1:2113837/71‑1 (MQ=255)

         catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                          >  1:1482523/1‑71 (MQ=255)

          atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                          <  1:439156/70‑1 (MQ=255)

           ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                        >  1:3290511/3‑71 (MQ=255)

                 ttACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                                                   <  1:2953341/70‑1 (MQ=255)

                  tACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                                                   >  1:683461/1‑69 (MQ=255)

                   aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                                                   >  1:4535093/1‑68 (MQ=255)

                    cGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                               <  1:2738835/71‑1 (MQ=255)

                     gcagcaGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCaagg                                              <  1:4350307/71‑4 (MQ=255)

                     gcagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                                                  <  1:373776/51‑1 (MQ=255)

                      cagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAgaga                                             <  1:2118396/71‑1 (MQ=255)

                       agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGATCTCATACTGCAGagaa                                            <  1:1570117/71‑1 (MQ=255)

                        gcaggcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAggaag                                           <  1:2366258/71‑1 (MQ=255)

                        gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                        <  1:1710426/58‑1 (MQ=255)

                        gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                        <  1:553177/58‑1 (MQ=255)

                        gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                        <  1:564540/58‑1 (MQ=255)

                        gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                           <  1:2007039/71‑1 (MQ=255)

                         caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                               <  1:154232/66‑1 (MQ=255)

                         caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                               <  1:312338/66‑1 (MQ=255)

                          aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                            >  1:1216229/1‑68 (MQ=255)

                          aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                            >  1:3710429/1‑68 (MQ=255)

                          aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                         <  1:316184/71‑1 (MQ=255)

                          aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                         >  1:3078264/1‑71 (MQ=255)

                           ggcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                        >  1:905889/2‑55 (MQ=255)

                           agcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                         <  1:3330683/70‑1 (MQ=255)

                            gcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgt                                       <  1:4835971/71‑1 (MQ=255)

                            gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                       <  1:419717/71‑1 (MQ=255)

                            gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                       <  1:1461097/71‑1 (MQ=255)

                            gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                       >  1:4337394/1‑71 (MQ=255)

                            gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                       <  1:1435471/71‑1 (MQ=255)

                             caaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATc                                                                  >  1:41501/1‑43 (MQ=12)

                              aaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgtga                                     <  1:4618897/71‑2 (MQ=255)

                               aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                    <  1:1551045/71‑1 (MQ=255)

                               aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                    >  1:2302400/1‑71 (MQ=255)

                               aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                    <  1:5588921/71‑1 (MQ=255)

                               aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                    <  1:1082422/71‑1 (MQ=255)

                                tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                               <  1:1803017/59‑1 (MQ=255)

                                tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                               <  1:1820379/59‑1 (MQ=255)

                                tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCg                                   <  1:5383566/71‑1 (MQ=255)

                                  gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                                                >  1:5055751/1‑40 (MQ=255)

                                  gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                                                >  1:4247255/1‑40 (MQ=255)

                                  gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                                                >  1:1150438/1‑40 (MQ=255)

                                  gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                                                           >  1:2581567/1‑45 (MQ=255)

                                  gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGAt                                 <  1:2751295/71‑1 (MQ=255)

                                  gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGa                                  <  1:940828/70‑1 (MQ=255)

                                      aGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATAc                                                     <  1:5126632/47‑1 (MQ=255)

                                       gTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCAt                                                       <  1:2293534/44‑1 (MQ=255)

                                        tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                              >  1:3663707/4‑68 (MQ=255)

                                        tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                              >  1:1970741/4‑68 (MQ=255)

                                        tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                              >  1:2152665/4‑68 (MQ=255)

                                        tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                              <  1:4903047/36‑1 (MQ=255)

                                        tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                                              <  1:3235701/36‑1 (MQ=255)

                                        tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                          <  1:235236/40‑1 (MQ=255)

                                        tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                              <  1:5667610/52‑1 (MQ=255)

                                          aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                         >  1:3379225/1‑71 (MQ=255)

                                          aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                         >  1:3179307/1‑71 (MQ=255)

                                          aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                         >  1:445539/1‑71 (MQ=255)

                                          aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                         >  1:705335/1‑71 (MQ=255)

                                           tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                           <  1:1334793/52‑1 (MQ=255)

                                           tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                           <  1:2456467/52‑1 (MQ=255)

                                           tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                           <  1:5640415/52‑1 (MQ=255)

                                           tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                         >  1:3110437/1‑70 (MQ=255)

                                            cGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                           <  1:2280545/51‑1 (MQ=255)

                                             gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa                       <  1:845802/70‑1 (MQ=255)

                                             gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAg                      >  1:5548385/1‑71 (MQ=255)

                                               tgttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                    >  1:456444/1‑71 (MQ=255)

                                                 ttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaaga                                          <  1:5177378/47‑1 (MQ=255)

                                                  tgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                              <  1:2901027/42‑1 (MQ=255)

                                                   gAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAgg                                        >  1:2996215/1‑47 (MQ=255)

                                                        ccGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAAt              <  1:1150458/68‑1 (MQ=255)

                                                         cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc           <  1:5337889/70‑1 (MQ=255)

                                                         cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc           <  1:2198931/70‑1 (MQ=255)

                                                           gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg        <  1:3077275/71‑1 (MQ=255)

                                                           gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg        <  1:5231204/71‑1 (MQ=255)

                                                            aCGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa       <  1:3909367/71‑1 (MQ=255)

                                                             cGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg        >  1:1036772/1‑69 (MQ=255)

                                                             cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTg        >  1:5175063/1‑69 (MQ=255)

                                                                 tCAGTTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc  <  1:3847743/71‑1 (MQ=255)

                                                                 tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                    <  1:611434/37‑1 (MQ=255)

                                                                 tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                    <  1:16782/37‑1 (MQ=255)

                                                                 tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg        <  1:2152486/65‑1 (MQ=255)

                                                                 tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc  <  1:4737692/71‑1 (MQ=255)

                                                                 tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc  <  1:2073074/71‑1 (MQ=255)

                                                                 tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc  <  1:4503201/71‑1 (MQ=255)

                                                                  cAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGtt   >  1:4818782/1‑69 (MQ=255)

                                                                   aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                    <  1:3710638/35‑1 (MQ=255)

                                                                   aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGa                                  >  1:5544640/1‑37 (MQ=255)

                                                                   aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa      >  1:2111126/1‑65 (MQ=255)

                                                                   aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa      >  1:1385118/1‑65 (MQ=255)

                                                                   aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa       >  1:2873580/1‑64 (MQ=255)

ACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTC  >  minE/608973‑609108

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