BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,760	G→A	32.2%	I498I (ATC→ATT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,760	0	G	A	32.2%	279.4 / 176.4	217	I498I (ATC→ATT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (74/73); new base A (28/42); total (102/115)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 242 total aligned reads displayed.

ATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACT  >  minE/608695‑608829

aTAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                                                   >  1:1649975/1‑70 (MQ=255)

  aaCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                                                                >  1:5158606/1‑71 (MQ=255)

     atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgata                                                                 >  1:1286047/2‑66 (MQ=255)

      tCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                               <  1:5486016/68‑1 (MQ=255)

       cACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                                                      >  1:5182136/1‑60 (MQ=255)

       cACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTc                                                              <  1:2432898/68‑1 (MQ=255)

       cACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTc                                                              <  1:378611/68‑1 (MQ=255)

        aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          <  1:2905936/71‑1 (MQ=255)

        aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          >  1:3543678/1‑71 (MQ=255)

        aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          >  1:2256811/1‑71 (MQ=255)

        aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          >  1:2192726/1‑71 (MQ=255)

        aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg                                                          >  1:323897/1‑69 (MQ=255)

         cgatCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcag                                                            >  1:1028952/4‑68 (MQ=255)

         cgatCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcag                                                            >  1:3621482/4‑68 (MQ=255)

           gTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcag                                                       >  1:3317045/1‑71 (MQ=255)

           gTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcag                                                       >  1:3504231/1‑71 (MQ=255)

            tCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgct                                                      <  1:3292929/71‑4 (MQ=255)

            tCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgct                                                      <  1:5602762/71‑4 (MQ=255)

              gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                                                     <  1:5414927/70‑1 (MQ=255)

              gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                    >  1:1885736/1‑71 (MQ=255)

                gCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                                                  >  1:172179/1‑71 (MQ=255)

                gCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                                                  >  1:5679924/1‑71 (MQ=255)

                gCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                                                   >  1:4880499/1‑54 (MQ=255)

                gCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                                                  >  1:5244114/1‑71 (MQ=255)

                gCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                                                  >  1:755507/1‑71 (MQ=255)

                 cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                                                  <  1:3837388/70‑1 (MQ=255)

                 cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGt                                                 <  1:470009/71‑1 (MQ=255)

                    cgcCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACt                                              <  1:1975400/71‑1 (MQ=255)

                    cgcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                                                   >  1:1561869/1‑50 (MQ=255)

                    cgcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                                                   >  1:2777947/1‑50 (MQ=255)

                    cgcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                                                   >  1:714355/1‑50 (MQ=255)

                     gcCATCAGAGAGTCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg                                             <  1:3744199/71‑1 (MQ=255)

                     gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg                                             <  1:2095007/71‑1 (MQ=255)

                       cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg                                             >  1:2107867/1‑69 (MQ=255)

                       cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt                                           <  1:609636/71‑1 (MQ=255)

                        aTCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTgcac                                          <  1:4290162/71‑4 (MQ=255)

                         tCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcag                                                            >  1:3070682/1‑52 (MQ=255)

                         tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgca                                                        <  1:2762348/56‑1 (MQ=255)

                           agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGttt                                          >  1:2964557/1‑68 (MQ=255)

                           agagagGCTTCGTAGCTGTACTTGTCGAGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                       >  1:3157256/1‑71 (MQ=255)

                              gagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccag                                    >  1:974728/1‑71 (MQ=255)

                               agGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          <  1:928533/48‑1 (MQ=255)

                                ggCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                                                                <  1:4176974/41‑1 (MQ=255)

                                    tCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTc                                                              >  1:536302/1‑39 (MQ=38)

                                    tCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTc                                                              >  1:2325123/1‑39 (MQ=38)

                                     cGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                                                            >  1:3665452/1‑40 (MQ=255)

                                        agatGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGGCAGCCAGTCCAt                           >  1:2373097/4‑70 (MQ=255)

                                        aGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCAt                           >  1:2913572/1‑70 (MQ=255)

                                         gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                    >  1:2600751/1‑44 (MQ=255)

                                         gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                    >  1:1453534/1‑44 (MQ=255)

                                         gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                    >  1:3949930/1‑44 (MQ=255)

                                         gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                   <  1:824590/61‑1 (MQ=255)

                                         gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                          >  1:5101739/1‑70 (MQ=255)

                                         gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                          >  1:4338760/1‑70 (MQ=255)

                                         gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                          <  1:1396075/70‑1 (MQ=255)

                                           tGTACTTGTCGTGCATGTAGTGAATGATATTcagcgc                                                         >  1:4978199/1‑37 (MQ=25)

                                           tGTACTTGTCGTGCATGTAGTGAATGATATTcagcgc                                                         >  1:119248/1‑37 (MQ=25)

                                            gTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGt                                            <  1:988713/49‑1 (MQ=255)

                                             tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                       <  1:2415906/53‑1 (MQ=255)

                                             tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                       <  1:1791742/53‑1 (MQ=255)

                                             tACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTg                     <  1:808091/71‑1 (MQ=255)

                                              aCTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGa                    >  1:252003/1‑71 (MQ=255)

                                               cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                       <  1:2815589/51‑1 (MQ=255)

                                               cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                       <  1:4455116/51‑1 (MQ=255)

                                               cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                   <  1:4667816/71‑1 (MQ=255)

                                                ttGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcca                                     <  1:531762/52‑1 (MQ=255)

                                                ttGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                  >  1:5420974/1‑71 (MQ=255)

                                                ttGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                                                  <  1:2321909/39‑1 (MQ=16)

                                                 tGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTTATGTACTGTTTagcca                                     >  1:4660148/1‑51 (MQ=255)

                                                 tGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcca                                     >  1:5567874/1‑51 (MQ=255)

                                                 tGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcca                                     >  1:4194749/1‑51 (MQ=255)

                                                 tGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcca                                     >  1:317845/1‑51 (MQ=255)

                                                 tGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                   >  1:3428816/1‑69 (MQ=255)

                                                 tGTCGTGCATGTAGTGGATGATGTTCACTGCAGTGATGTACTGTTTagcca                                     >  1:395611/1‑51 (MQ=255)

                                                 tGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATc                  >  1:1219189/1‑70 (MQ=255)

                                                  gTCTTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATc                  <  1:1123840/69‑1 (MQ=255)

                                                  gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATcc                 <  1:2236535/70‑1 (MQ=255)

                                                  gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                  <  1:5453829/69‑1 (MQ=255)

                                                  gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                  <  1:4102748/69‑1 (MQ=255)

                                                  gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                  <  1:2985682/69‑1 (MQ=255)

                                                  gTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTcc                             <  1:481885/58‑1 (MQ=255)

                                                   tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                  <  1:3352880/68‑1 (MQ=255)

                                                   tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCAt               >  1:4930906/1‑71 (MQ=255)

                                                   tCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATcga                <  1:2636195/70‑3 (MQ=255)

                                                    cGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATg              >  1:5019723/1‑71 (MQ=255)

                                                       gCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagcc                                  <  1:895014/48‑1 (MQ=255)

                                                       gCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagcc                                  <  1:1966901/48‑1 (MQ=255)

                                                        cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCt          <  1:3455623/71‑1 (MQ=255)

                                                          tGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                          <  1:3589835/53‑1 (MQ=255)

                                                          tGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc         >  1:4003165/1‑70 (MQ=255)

                                                            tAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa       <  1:4580983/70‑1 (MQ=255)

                                                            tAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCAt      <  1:4527083/71‑1 (MQ=255)

                                                            tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGaa                        >  1:1852218/1‑53 (MQ=12)

                                                            tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGaa                        >  1:5577554/1‑53 (MQ=12)

                                                            tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATc                  <  1:1008395/59‑1 (MQ=255)

                                                              gTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCa    <  1:2696829/71‑1 (MQ=255)

                                                                 gatgatGTTCAGTGCAGTGATGTACTGTTTagccag                                    <  1:921637/36‑1 (MQ=255)

                                                                 gatgatGTTCAGTGCAGTGATGTACTGTTTagccag                                    <  1:3014893/36‑1 (MQ=255)

                                                                 gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGa                         >  1:5534629/1‑47 (MQ=255)

                                                                 gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                   >  1:1952591/1‑53 (MQ=255)

ATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACT  >  minE/608695‑608829

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