BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,777	A→G	22.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	22.1%	193.4 / 51.6	102	R159R (CGT→CGC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (45/34); new base G (16/6); total (62/40)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

Only 100 of 107 total aligned reads displayed.

GGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATT  >  minE/609707‑609846

ggTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCaa                                                                       >  1:277157/1‑71 (MQ=255)

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

    gTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCg                                                                     >  1:1254477/1‑69 (MQ=255)

     tCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAg                                                                  <  1:2118094/71‑1 (MQ=255)

        ccGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                                                               <  1:1036060/71‑1 (MQ=255)

         cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                                  >  1:2603756/1‑67 (MQ=255)

         cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAg                                                                  >  1:762822/1‑67 (MQ=255)

         cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATg                                                              >  1:1591625/1‑71 (MQ=255)

          gatgatACGGCCACGGCCATATGGATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGt                                                             >  1:3618208/1‑71 (MQ=255)

          gataatGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAtat                                                             <  1:3795941/71‑3 (MQ=255)

           atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc                                                                    >  1:563614/1‑63 (MQ=255)

           atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc                                                                    >  1:1713932/1‑63 (MQ=255)

           atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc                                                            >  1:8738/1‑71 (MQ=255)

           atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAAGAtag                                                             >  1:1777735/1‑68 (MQ=255)

              atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCgg                                                         <  1:2110042/71‑1 (MQ=255)

              atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                                          <  1:264934/70‑1 (MQ=255)

               tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                                          >  1:403082/1‑69 (MQ=255)

               tACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                                          >  1:1819363/1‑69 (MQ=255)

                  ggccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAg                                                       <  1:3154099/69‑1 (MQ=255)

                   gccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCa                                                                   <  1:681140/56‑1 (MQ=255)

                   cccacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAgc                                                                 >  1:2368655/2‑57 (MQ=255)

                    ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAgtgt                                                    >  1:2180042/1‑70 (MQ=255)

                    ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAgt                                                      >  1:2812195/1‑68 (MQ=255)

                     cacggccaTATGCATCTGTCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                                  <  1:5535/71‑1 (MQ=255)

                     cacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaa                                                  <  1:3006292/71‑1 (MQ=255)

                      acggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                 >  1:1205796/1‑71 (MQ=255)

                       cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                   <  1:1759126/68‑1 (MQ=255)

                       cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                   >  1:2006983/1‑68 (MQ=255)

                       cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                 >  1:2683014/1‑70 (MQ=255)

                       cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                 >  1:2393813/1‑70 (MQ=255)

                       cggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                 >  1:2389073/1‑70 (MQ=255)

                        tgccaGATGCATATGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                               >  1:2894211/2‑71 (MQ=255)

                        ggccaTATGCATCTGGCAGACAGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                               <  1:1154638/71‑1 (MQ=255)

                          ccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAAcat                                              >  1:3254109/1‑68 (MQ=255)

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

                              atGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGa                                                        >  1:2681569/1‑56 (MQ=255)

                              atGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGa                                                        >  1:3667396/1‑56 (MQ=255)

                               tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                 >  1:89121/1‑62 (MQ=255)

                               tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                        <  1:1557305/71‑1 (MQ=255)

                               tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                        <  1:3713503/71‑1 (MQ=255)

                               tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                        <  1:3863764/71‑1 (MQ=255)

                                gCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGa                                                        >  1:2082483/1‑54 (MQ=255)

                                ccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCggt                                                        >  1:1971649/2‑53 (MQ=16)

                                 cATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACg                                      <  1:3294142/71‑1 (MQ=255)

                                 cATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                                  >  1:3389602/1‑59 (MQ=255)

                                 cATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                                  >  1:2722572/1‑59 (MQ=255)

                                  aTCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATcaa                                           <  1:2573065/65‑3 (MQ=255)

                                  aTCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGc                                     >  1:934271/1‑71 (MQ=255)

                                  aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAg                                                       >  1:84688/1‑53 (MQ=255)

                                  aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                           >  1:622401/1‑65 (MQ=255)

                                  aTCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                           >  1:3472631/1‑65 (MQ=255)

                                   tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAg                                                                  >  1:2598261/1‑41 (MQ=19)

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

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

                                    cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                   >  1:3202735/1‑71 (MQ=255)

                                    cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                   >  1:683712/1‑71 (MQ=255)

                                    cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                   >  1:3124160/1‑71 (MQ=255)

                                    cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                   >  1:1200076/1‑71 (MQ=255)

                                      ggCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                           <  1:3232434/61‑1 (MQ=255)

                                       gCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGa                                           <  1:2694772/60‑1 (MQ=255)

                                         aGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                   <  1:1952353/50‑1 (MQ=255)

                                           aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                          <  1:540759/57‑1 (MQ=255)

                                           aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                            <  1:794064/71‑1 (MQ=255)

                                            ccGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAAc                                                >  1:3332771/1‑50 (MQ=255)

                                            ccGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgtg                           <  1:1080095/71‑1 (MQ=255)

                                             cGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                          <  1:472274/71‑1 (MQ=255)

                                              ggTCAGAACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCACTGGTtatggg                         >  1:2884392/1‑69 (MQ=255)

                                               gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                          >  1:2894459/1‑69 (MQ=255)

                                                tCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGGCGAACACGCCCTGGTTGTGAGtt                       <  1:1754112/71‑1 (MQ=255)

                                                 cAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                      >  1:10591/1‑71 (MQ=255)

                                                  aGAACTCCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt                     <  1:2585645/71‑1 (MQ=255)

                                                  aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTTTAAACGTCGAACACGCCCTGGGTGTGAGtttt                     <  1:3339867/71‑1 (MQ=255)

                                                  aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                                            <  1:3287562/48‑1 (MQ=255)

                                                  aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                      <  1:186373/70‑1 (MQ=255)

                                                  aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                      <  1:2490538/70‑1 (MQ=255)

                                                  aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                      <  1:858471/70‑1 (MQ=255)

                                                    aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                    <  1:816445/70‑1 (MQ=255)

                                                     acgccAGATTTACGGCATCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                  >  1:1815777/4‑71 (MQ=255)

                                                      cacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACgg                 >  1:437384/1‑71 (MQ=255)

                                                      cacCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACgg                 >  1:1129451/1‑71 (MQ=255)

                                                       gccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGttt                      <  1:1631026/64‑1 (MQ=255)

                                                       gccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                  >  1:2061190/2‑69 (MQ=255)

                                                          aGATTTACGGCAGCGCAGCATATCCGGTGAGTAAAc                                                >  1:3579797/1‑36 (MQ=12)

                                                          aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGtt                              >  1:1452631/1‑54 (MQ=255)

                                                           gATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                  >  1:1326924/1‑65 (MQ=255)

                                                           gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACg                                      >  1:1028313/1‑45 (MQ=255)

                                                           gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGt                >  1:1099139/1‑67 (MQ=255)

                                                           gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTc            <  1:1668607/71‑1 (MQ=255)

                                                              ttACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTg                                 >  1:2056657/1‑47 (MQ=255)

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

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

                                                               tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                                            <  1:627472/35‑1 (MQ=255)

                                                               tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGt         <  1:2080340/70‑1 (MQ=255)

                                                                aCGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg        <  1:3194530/70‑1 (MQ=255)

                                                                aCGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg        <  1:2314118/70‑1 (MQ=255)

                                                                 cGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGa       <  1:1607697/70‑1 (MQ=255)

                                                                     aaCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtt  >  1:1050678/1‑71 (MQ=255)

                                                                      aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                   >  1:1502410/1‑37 (MQ=255)

                                                                      aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                   >  1:649819/1‑37 (MQ=255)

                                                                      aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGa    >  1:3117557/1‑68 (MQ=255)

GGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATT  >  minE/609707‑609846

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`