BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,759	A→C	10.9%	V165V (GTT→GTG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,759	0	A	C	10.9%	230.4 / 30.4	101	V165V (GTT→GTG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (48/42); new base C (9/2); total (57/44)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

Only 100 of 109 total aligned reads displayed.

GTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTAC  >  minE/609690‑609829

gTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                                                                       >  1:949343/1‑71 (MQ=255)

   cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                                                                       >  1:244799/1‑68 (MQ=255)

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

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

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

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

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

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

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

            gcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTAc                                                           >  1:2444628/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:1260641/71‑1 (MQ=255)

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

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

                   tAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                                                                       >  1:3271569/1‑52 (MQ=255)

                   tAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                                                                       >  1:2902626/1‑52 (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)

                      tCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGa                                                                        >  1:2817868/1‑48 (MQ=255)

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

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

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

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

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

                        aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTAc                                                           >  1:2257151/1‑59 (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)

                           gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAac                                                                      >  1:3220009/1‑45 (MQ=255)

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

                           gatgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                            >  1:426803/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)

                                         ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                              <  1:2958091/71‑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:410682/45‑1 (MQ=255)

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

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

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

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

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

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

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

                                                gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                         >  1:721075/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:4108083/1‑56 (MQ=255)

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

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

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

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

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

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

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

                                                              cGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                            <  1:3565293/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)

GTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTAC  >  minE/609690‑609829

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