BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,000	A→G	6.8%	A418A (GCT→GCC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,000	0	A	G	6.8%	293.7 / 10.7	103	A418A (GCT→GCC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (63/33); new base G (4/3); total (67/36)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 135 total aligned reads displayed.

TCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGA  >  minE/608933‑609068

tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                                    <  1:3395187/70‑1 (MQ=255)

 cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                                                   <  1:1199145/70‑1 (MQ=255)

 cGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                                                   <  1:1670833/70‑1 (MQ=255)

 cGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCagcaggca                                                                   <  1:1900468/70‑1 (MQ=255)

    cAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg                                                               <  1:1472220/71‑1 (MQ=255)

      ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAg                                                               >  1:4078426/1‑69 (MQ=255)

      ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                                             <  1:2076321/71‑1 (MQ=255)

      ggTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgcg                                                             >  1:1221440/1‑71 (MQ=255)

       gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGt                                                            <  1:1555201/71‑1 (MQ=255)

        tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTa                                                           <  1:232203/71‑1 (MQ=255)

         ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                                  >  1:1117979/1‑63 (MQ=255)

         ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                                                >  1:2022989/1‑65 (MQ=255)

          tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                                         >  1:2004154/1‑71 (MQ=255)

          tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                                         >  1:2525102/1‑71 (MQ=255)

             aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                                    <  1:3363694/57‑1 (MQ=255)

             aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCAt                                                      <  1:3381657/71‑1 (MQ=255)

               gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                                  >  1:282008/1‑57 (MQ=255)

               gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCg                                                    <  1:2744558/71‑1 (MQ=255)

                tgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTa                                                           <  1:456201/62‑1 (MQ=255)

                tgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATAGCGTa                                                           <  1:802870/62‑1 (MQ=255)

                 gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                                  >  1:3364338/1‑71 (MQ=255)

                 gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                                  >  1:1211478/1‑71 (MQ=255)

                 gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgct                                                  >  1:1385780/1‑69 (MQ=255)

                 gcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgct                                                  >  1:808507/1‑69 (MQ=255)

                    ccGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                                    <  1:2944900/50‑1 (MQ=255)

                    ccGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                >  1:123820/1‑70 (MQ=255)

                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAAt                                                                 <  1:1006285/52‑1 (MQ=255)

                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc                                                     <  1:965716/64‑1 (MQ=255)

                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt                                               >  1:2056076/1‑70 (MQ=255)

                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt                                               >  1:2061036/1‑70 (MQ=255)

                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                              <  1:151874/71‑1 (MQ=255)

                     cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                              <  1:1725687/71‑1 (MQ=255)

                      gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt                                               >  1:1840749/1‑69 (MQ=255)

                      gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                                                 >  1:528020/1‑67 (MQ=255)

                      gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGAAGCAAGCAATAGCGTAGTCATCgttg                                                 >  1:2618494/1‑67 (MQ=255)

                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                              <  1:618985/69‑1 (MQ=255)

                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                              >  1:2744738/1‑69 (MQ=255)

                       aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                            <  1:2852718/71‑1 (MQ=255)

                       aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGaa                                            <  1:1272267/71‑1 (MQ=255)

                        agaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAg                                           >  1:3157389/1‑71 (MQ=255)

                           aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTcc                                        <  1:1934522/71‑1 (MQ=255)

                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc                                                                    <  1:4002379/42‑1 (MQ=255)

                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                                  >  1:1850420/1‑44 (MQ=255)

                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                       <  1:3254358/71‑1 (MQ=255)

                            cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                       <  1:421580/71‑1 (MQ=255)

                              gCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                >  1:1060756/1‑60 (MQ=255)

                               cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                                  <  1:244663/41‑1 (MQ=255)

                               cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                                  <  1:1817824/41‑1 (MQ=255)

                               cATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                                  <  1:3988942/41‑1 (MQ=255)

                                 tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAg                                           <  1:643970/62‑1 (MQ=255)

                                 tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGc                                  <  1:3672883/71‑1 (MQ=255)

                                   tGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                                  >  1:2006482/1‑37 (MQ=255)

                                   tGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                                  >  1:4081085/1‑37 (MQ=255)

                                   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTc                                         >  1:56076/1‑62 (MQ=255)

                                   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                                >  1:3799987/1‑71 (MQ=255)

                                   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                                >  1:2830485/1‑71 (MQ=255)

                                   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                                >  1:2190831/1‑71 (MQ=255)

                                   tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                                >  1:259203/1‑71 (MQ=255)

                                    gTTTACCAACGATCTTCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                                <  1:2443336/70‑1 (MQ=255)

                                           aaCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCa                                                       <  1:2976428/40‑1 (MQ=255)

                                             cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc                                                     >  1:1723181/1‑40 (MQ=255)

                                             cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc                                                     >  1:2747646/1‑40 (MQ=255)

                                             cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                     >  1:410887/1‑56 (MQ=255)

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

                                                 catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                                                 >  1:708533/1‑40 (MQ=255)

                                                  atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCa                         <  1:1326218/63‑1 (MQ=255)

                                                  atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGtt                    <  1:1624796/68‑1 (MQ=255)

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

                                                  atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                  >  1:3265333/1‑70 (MQ=255)

                                                  atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                  >  1:3056759/1‑70 (MQ=255)

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

                                                   tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                              <  1:2833879/41‑1 (MQ=255)

                                                   tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                            >  1:1508039/1‑43 (MQ=255)

                                                   tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                            >  1:1983774/1‑43 (MQ=255)

                                                   tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                <  1:3865728/71‑1 (MQ=255)

                                                    cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                >  1:2351553/1‑70 (MQ=255)

                                                    cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                >  1:3344593/1‑70 (MQ=255)

                                                    cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                >  1:49787/1‑70 (MQ=255)

                                                    cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                >  1:789443/1‑70 (MQ=255)

                                                    cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                >  1:3877880/1‑70 (MQ=255)

                                                    cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                >  1:3647067/1‑70 (MQ=255)

                                                     gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                            >  1:3452161/1‑41 (MQ=255)

                                                     gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                     >  1:468891/1‑48 (MQ=255)

                                                     gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa              >  1:1970506/1‑71 (MQ=255)

                                                      ggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATAc             <  1:1611743/71‑1 (MQ=255)

                                                       gCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACt            >  1:2839531/1‑71 (MQ=255)

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

                                                            cGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCag        <  1:492426/70‑1 (MQ=255)

                                                             gcagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCa         >  1:2162635/1‑68 (MQ=255)

                                                             gcagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCa         >  1:581631/1‑68 (MQ=255)

                                                             gcagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCa         >  1:2606828/1‑68 (MQ=255)

                                                               agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATc                               >  1:1150044/1‑44 (MQ=255)

                                                               agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATc                               >  1:1031630/1‑44 (MQ=255)

                                                               agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag      >  1:3108699/1‑69 (MQ=255)

                                                               agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa    <  1:182696/71‑1 (MQ=255)

                                                               agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa    <  1:112456/71‑1 (MQ=255)

                                                               agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa    <  1:2087585/71‑1 (MQ=255)

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

                                                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTGTCATa              >  1:2334013/1‑58 (MQ=255)

                                                                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa              >  1:503913/1‑58 (MQ=255)

TCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGA  >  minE/608933‑609068

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