BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,762	A→G	11.0%	G164G (GGT→GGC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,762	0	A	G	11.0%	265.0 / 22.0	100	G164G (GGT→GGC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (47/42); new base G (9/2); total (56/44)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

Only 100 of 104 total aligned reads displayed.

GCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACG  >  minE/609695‑609830

gcgcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCa                                                                   <  1:77437/71‑1 (MQ=255)

 cgcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAacac                                                                     >  1:1221230/1‑68 (MQ=255)

 cgcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCa                                                                   <  1:2297647/70‑1 (MQ=255)

     aCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAtt                                                               <  1:80781/70‑1 (MQ=255)

     aCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAtt                                                               <  1:2139000/70‑1 (MQ=255)

     aCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAtt                                                               <  1:3316350/70‑1 (MQ=255)

       gcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTAc                                                            >  1:2444628/1‑71 (MQ=255)

       gcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTAc                                                            <  1:2527331/71‑1 (MQ=255)

       gcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTAc                                                            >  1:2742498/1‑71 (MQ=255)

       gcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTAc                                                            >  1:3102579/1‑71 (MQ=255)

         gACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACgg                                                          >  1:900120/1‑71 (MQ=255)

           cGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGc                                                         >  1:1514253/1‑70 (MQ=255)

            ggTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCaa                                                       <  1:161219/71‑1 (MQ=255)

            ggTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCaa                                                       <  1:3126426/71‑1 (MQ=255)

            ggTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCaa                                                       <  1:1260641/71‑1 (MQ=255)

             gTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCa                                                                   <  1:2574088/58‑1 (MQ=255)

              tAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACg                                                     <  1:3531195/71‑1 (MQ=255)

              tAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCg                                                     >  1:1785669/1‑71 (MQ=255)

               aGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGc                                                         <  1:2059451/66‑1 (MQ=255)

                gTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc                                                    <  1:2588397/70‑1 (MQ=255)

                gTCACCGATGATACGGCCACGGCCATATGCATCTAGCAGACCGGTCAGAACACCAGATTTACGGCAACGc                                                    <  1:112071/70‑1 (MQ=255)

                 tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                  <  1:898035/71‑1 (MQ=255)

                 tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                  <  1:3672418/71‑1 (MQ=255)

                 tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                  <  1:671543/71‑1 (MQ=255)

                 tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                  <  1:3719681/71‑1 (MQ=255)

                   aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTAc                                                            >  1:2257151/1‑59 (MQ=255)

                   aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                >  1:1913836/1‑71 (MQ=255)

                   aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                >  1:2124549/1‑71 (MQ=255)

                     cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                >  1:589932/1‑69 (MQ=255)

                     cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATg                                              <  1:906833/71‑1 (MQ=255)

                      gatgatACGGCCACGGCCATATGCATCTGGCAGCCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                             >  1:1728849/1‑71 (MQ=255)

                      gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                                               >  1:1716252/1‑69 (MQ=255)

                      gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                             >  1:426803/1‑71 (MQ=255)

                      gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                             >  1:170203/1‑71 (MQ=255)

                      gatgatACGGCCACGGCCATATGCATCTGGCAGACCGCTCAGAacac                                                                     >  1:3215177/1‑47 (MQ=255)

                       atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                  <  1:2665645/65‑1 (MQ=255)

                         gatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTAc                                                            <  1:766987/53‑1 (MQ=255)

                           tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAtt                                                               <  1:373277/48‑1 (MQ=255)

                           tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAgg                                                 >  1:312090/1‑62 (MQ=255)

                               gccacggccaTATGCATCTGGCAGACCGGTCAGAAGACCAGATTTACGGCAACGCAGGATGTCCGGAgtg                                     >  1:530807/1‑70 (MQ=255)

                               gccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACg                                                     >  1:2742079/1‑54 (MQ=255)

                               gccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                  <  1:861462/57‑1 (MQ=255)

                                ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                   <  1:3791377/71‑1 (MQ=255)

                                 cacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCggtga                                      >  1:1923519/1‑64 (MQ=255)

                                 cacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                                   <  1:293023/70‑1 (MQ=255)

                                     gccaTATGCATCTGGCAGACCGGTCAGAACACCAGAttt                                                              <  1:3075651/39‑1 (MQ=255)

                                     gccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGc                                                         >  1:336327/1‑44 (MQ=255)

                                     gccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                 >  1:2523725/1‑68 (MQ=255)

                                        ttatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                   <  1:4005554/62‑1 (MQ=255)

                                        atatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                           <  1:1666750/71‑1 (MQ=255)

                                         tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                  <  1:3007836/47‑1 (MQ=255)

                                         tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                  <  1:677384/47‑1 (MQ=255)

                                         tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                          <  1:1281979/71‑1 (MQ=255)

                                          atGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCgg                                         <  1:3523259/55‑1 (MQ=255)

                                           tGCATCTGGCAGACCGGTCAGAACACCAGATTTAc                                                            >  1:395313/1‑35 (MQ=255)

                                           tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCa                                                   <  1:2695245/44‑1 (MQ=255)

                                           tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                  >  1:2674695/1‑45 (MQ=255)

                                           tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                  <  1:410682/45‑1 (MQ=255)

                                           tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAg                                       >  1:2575142/1‑56 (MQ=255)

                                           tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                        >  1:1758809/1‑71 (MQ=255)

                                           tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                        >  1:1967169/1‑71 (MQ=255)

                                           tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                        >  1:2272082/1‑71 (MQ=255)

                                           tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                        >  1:2561492/1‑71 (MQ=255)

                                           gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                          >  1:721075/3‑69 (MQ=255)

                                           gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                          >  1:1757432/3‑69 (MQ=255)

                                           gccATCTGGTAAACAGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                          >  1:1214781/3‑69 (MQ=255)

                                             cATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATca                            >  1:1686080/1‑64 (MQ=255)

                                             cATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAgg                                                 >  1:541739/1‑44 (MQ=255)

                                               tCTGGCAGAGCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCgg                                         <  1:2699775/50‑1 (MQ=255)

                                               tCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGcc                    >  1:2871227/1‑71 (MQ=255)

                                               tCTGGCAGACC‑GTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                   >  1:3835619/1‑71 (MQ=255)

                                                cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                 >  1:3630091/1‑57 (MQ=255)

                                                 tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                                               >  1:674260/1‑42 (MQ=255)

                                                 tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                                               >  1:679048/1‑42 (MQ=255)

                                                 tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                 >  1:4131320/1‑56 (MQ=255)

                                                 tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                 >  1:4108083/1‑56 (MQ=255)

                                                 tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt                  >  1:2667979/1‑71 (MQ=255)

                                                 tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt                  >  1:1463254/1‑71 (MQ=255)

                                                    cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGt               >  1:3626158/1‑71 (MQ=255)

                                                     aaaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGt               >  1:2203540/3‑70 (MQ=255)

                                                     aaaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGt               >  1:588441/3‑70 (MQ=255)

                                                     aGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                  <  1:3335961/35‑1 (MQ=255)

                                                     aGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt              >  1:583371/1‑71 (MQ=255)

                                                     aGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt              >  1:2484232/1‑71 (MQ=255)

                                                      gACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                          <  1:482496/58‑1 (MQ=255)

                                                      gACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                          <  1:3222016/58‑1 (MQ=255)

                                                      gACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtg             >  1:2602522/1‑71 (MQ=255)

                                                        ccGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTtatg           >  1:916180/1‑71 (MQ=255)

                                                        ccGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt            >  1:635880/1‑70 (MQ=255)

                                                         cGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                             <  1:3565293/36‑1 (MQ=255)

                                                         cGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                             <  1:1073497/36‑1 (MQ=255)

                                                         cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt            >  1:1938058/1‑69 (MQ=255)

                                                         cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt            >  1:1610704/1‑69 (MQ=255)

                                                            tCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaa                                  <  1:1840715/44‑1 (MQ=12)

                                                               gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                        >  1:2404495/1‑51 (MQ=255)

                                                               gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt     <  1:32959/70‑1 (MQ=255)

                                                                aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                 >  1:3022398/1‑41 (MQ=255)

                                                                aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTAc   >  1:677873/1‑71 (MQ=255)

                                                                 acacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg  <  1:3980117/71‑1 (MQ=255)

                                                                 acacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACg  <  1:1709024/71‑1 (MQ=255)

GCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACG  >  minE/609695‑609830

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