BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,820	C→G	10.0%	E478D (GAG→GAC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,820	0	C	G	10.0%	281.9 / 17.3	100	E478D (GAG→GAC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (49/41); new base G (4/6); total (53/47)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 109 total aligned reads displayed.

CATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTT  >  minE/608751‑608888

cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATgcgc                                                                      <  1:35490/70‑1 (MQ=255)

cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCt                                                                     >  1:5178894/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

      gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCa                                                               >  1:2323210/1‑71 (MQ=255)

       tgaatgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCAc                                                              <  1:159215/68‑1 (MQ=255)

       tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCAc                                                              <  1:954640/71‑1 (MQ=255)

       tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCAc                                                              <  1:2240940/71‑1 (MQ=255)

         gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACt                                                             <  1:2768294/70‑1 (MQ=255)

         aatgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtt                                                            >  1:1423752/2‑71 (MQ=255)

          atgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACttt                                                           >  1:4764583/1‑70 (MQ=255)

          atgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACttt                                                           <  1:4779975/71‑2 (MQ=255)

           tgatgTTCAGTGCAGTGATGTACTGTTTACCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtca                                                          <  1:2401335/71‑1 (MQ=255)

           tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg                                                          <  1:2684964/71‑3 (MQ=255)

           tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg                                                          <  1:1182116/71‑3 (MQ=255)

           tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg                                                          >  1:3939805/1‑69 (MQ=255)

              tgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca                                                       <  1:1855108/71‑1 (MQ=255)

                ttCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca                                                       >  1:5059907/1‑69 (MQ=255)

                ttCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca                                                       >  1:4852286/1‑69 (MQ=255)

                ttCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCa                                                               <  1:5063728/61‑1 (MQ=255)

                 tCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcat                                                         >  1:1206162/1‑66 (MQ=255)

                 tCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                                                    >  1:4597266/1‑71 (MQ=255)

                 tCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                                                    >  1:1383646/1‑71 (MQ=255)

                  cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                                                    >  1:528916/1‑70 (MQ=255)

                  cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                                                    >  1:2147734/1‑70 (MQ=255)

                  cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                                                    >  1:5662568/1‑70 (MQ=255)

                  cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGt                                                   >  1:87301/1‑71 (MQ=255)

                    gtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTc                                                 <  1:5513822/71‑1 (MQ=255)

                    gtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTc                                                 <  1:1142408/71‑1 (MQ=255)

                    gtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTc                                                 <  1:3086584/71‑1 (MQ=255)

                       cagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAg                                               <  1:3585521/70‑1 (MQ=255)

                       cagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAg                                               <  1:5444385/70‑1 (MQ=255)

                         gtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACtt                                                            <  1:4432361/55‑1 (MQ=255)

                         gtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAc                                            <  1:2126794/71‑1 (MQ=255)

                         gtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAc                                            <  1:5240645/71‑1 (MQ=255)

                          tgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCa                                                <  1:2086920/66‑1 (MQ=255)

                          tgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAc                                            <  1:1278528/70‑1 (MQ=255)

                          tgATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCa                                                               <  1:5043685/51‑1 (MQ=255)

                            aTGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAg                                               >  1:3079151/1‑65 (MQ=255)

                               tACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                                                  <  1:4835358/43‑1 (MQ=255)

                                 cTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCtt                                    >  1:295020/1‑71 (MQ=255)

                                 cTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCtt                                    >  1:1035317/1‑71 (MQ=255)

                                  tGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCttt                                   >  1:687099/1‑71 (MQ=255)

                                   gTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                                                  >  1:4808103/1‑39 (MQ=255)

                                   gTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                                                  >  1:524488/1‑39 (MQ=255)

                                   gTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                                                  >  1:1775835/1‑39 (MQ=255)

                                   gTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCttt                                   >  1:610017/1‑70 (MQ=255)

                                      tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCa                                                               >  1:1327243/1‑39 (MQ=255)

                                      tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCa                                                               >  1:2599519/1‑39 (MQ=255)

                                      tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATc                               >  1:2914802/1‑71 (MQ=255)

                                      tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATc                               >  1:5209155/1‑71 (MQ=255)

                                        gccagccaGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCgg                             <  1:772482/71‑1 (MQ=255)

                                        gccagccaGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCgg                             <  1:4783903/71‑1 (MQ=255)

                                         ccagccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGt                            >  1:3239341/1‑71 (MQ=255)

                                            gccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                           >  1:5283357/1‑69 (MQ=255)

                                            gccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                           >  1:3053173/1‑69 (MQ=255)

                                            gccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCa                         <  1:4025166/71‑1 (MQ=255)

                                             ccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                           >  1:658166/1‑68 (MQ=255)

                                             ccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                           >  1:1093147/1‑68 (MQ=255)

                                              cagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGa                       >  1:3063933/1‑71 (MQ=255)

                                                gTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTg                                  <  1:3393067/58‑1 (MQ=255)

                                                gTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGAc                      >  1:1429563/1‑70 (MQ=255)

                                                   cATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCg                  >  1:2315182/1‑71 (MQ=255)

                                                   cATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCg                  >  1:5658349/1‑71 (MQ=255)

                                                   cATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCg                  >  1:3291293/1‑71 (MQ=255)

                                                     tGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGa                >  1:3450673/1‑71 (MQ=255)

                                                     tGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGa                >  1:2971293/1‑71 (MQ=255)

                                                      gAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAc                                            >  1:4064017/1‑42 (MQ=255)

                                                      gAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCgg                             >  1:3950326/1‑57 (MQ=255)

                                                      gAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCg                  >  1:4070583/1‑68 (MQ=255)

                                                      gAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGAc               >  1:621988/1‑71 (MQ=255)

                                                        aGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTc                   >  1:258179/1‑65 (MQ=255)

                                                         gTGATCCATGGGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCaa            >  1:2444895/1‑71 (MQ=255)

                                                         gTGATCCATGCGCTCCATCACTTCATCATAGTTCAg                                               <  1:4032457/36‑1 (MQ=255)

                                                         gTGATCCATGCGCTCCATCACTTCATCATAGTTCAg                                               <  1:4025868/36‑1 (MQ=255)

                                                         gTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCaa            >  1:5097586/1‑71 (MQ=255)

                                                          tGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTg                                  <  1:1825032/48‑1 (MQ=255)

                                                          tGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGa                       >  1:1382234/1‑59 (MQ=255)

                                                          tGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAAc           >  1:30573/1‑71 (MQ=255)

                                                             tCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAAc           <  1:1831273/68‑1 (MQ=255)

                                                             tCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAAc           <  1:5052332/68‑1 (MQ=255)

                                                             tCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTg        <  1:1473586/71‑1 (MQ=255)

                                                              ccATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTc                   >  1:5469228/1‑59 (MQ=255)

                                                               cATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTc                   <  1:3597882/58‑1 (MQ=255)

                                                               cATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTc                   <  1:669781/58‑1 (MQ=255)

                                                                aTGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCa                         <  1:4017214/51‑1 (MQ=255)

                                                                 tGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAtt    <  1:2256830/71‑1 (MQ=255)

                                                                  gcgcTCCATCACTTCATCATAGTTCAGGACATCGCCt                                     >  1:361244/1‑37 (MQ=255)

                                                                  gcgcTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAttt   >  1:337572/1‑71 (MQ=255)

                                                                   cgcTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAtttt  >  1:204513/1‑71 (MQ=255)

                                                                   cgcTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAttt   >  1:5403660/1‑70 (MQ=255)

                                                                   cgcTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAttt   >  1:4786776/1‑70 (MQ=255)

                                                                   cgcTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAttt   >  1:3546414/1‑70 (MQ=255)

                                                                    gcTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAttt   >  1:2525482/1‑69 (MQ=255)

CATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTT  >  minE/608751‑608888

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