BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,777	A→G	15.1%	R159R (CGT→CGC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,777	0	A	G	15.1%	246.9 / 33.5	99	R159R (CGT→CGC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (47/37); new base G (10/5); total (57/42)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 102 total aligned reads displayed.

GGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAT  >  minE/609707‑609845

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

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

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

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

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

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

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

     tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGGAACGc                                                                   <  1:4016960/69‑1 (MQ=255)

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

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

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

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

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

       aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                               >  1:1913836/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:170203/1‑71 (MQ=255)

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

           atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                                 <  1:2665645/65‑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)

                         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)

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

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

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

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

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

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

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

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

                               tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                       >  1:2272082/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:679048/1‑42 (MQ=255)

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

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

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

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

                                     tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt                                 >  1:2667979/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:3222016/58‑1 (MQ=255)

                                          gACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                         <  1:482496/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:1073497/36‑1 (MQ=255)

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

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

                                             cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                           >  1:1938058/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)

                                                     acgccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                 <  1:241043/68‑1 (MQ=255)

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

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

                                                      cgccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                 >  1:952727/3‑70 (MQ=255)

                                                      cgccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                 <  1:1354525/68‑1 (MQ=255)

                                                        ccAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                         <  1:2034711/44‑1 (MQ=255)

                                                        ccAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTa              >  1:2956282/1‑71 (MQ=255)

                                                         cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAGACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt             >  1:3123276/1‑71 (MQ=255)

                                                         cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt             >  1:1292000/1‑71 (MQ=255)

                                                         cAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAt             >  1:3127472/1‑71 (MQ=255)

                                                          aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACg                                     >  1:3724182/1‑46 (MQ=255)

                                                          aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt            >  1:2109699/1‑71 (MQ=255)

                                                           gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTgg                               <  1:746411/51‑1 (MQ=255)

                                                            aTTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGttt                     <  1:3061216/60‑1 (MQ=12)

                                                            aTTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCa          <  1:2775210/71‑1 (MQ=255)

                                                              ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt        >  1:838025/1‑71 (MQ=255)

                                                              ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt        >  1:3992948/1‑71 (MQ=255)

                                                                 cGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAg         <  1:2792215/67‑1 (MQ=255)

                                                                 cGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAg         <  1:1061055/67‑1 (MQ=255)

                                                                   gCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                         <  1:577511/49‑1 (MQ=255)

                                                                    cAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAt  <  1:787331/71‑1 (MQ=255)

GGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAT  >  minE/609707‑609845

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