BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,015	A→G	5.4%	D413D (GAT→GAC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,015	0	A	G	5.4%	343.5 / 6.7	117	D413D (GAT→GAC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (37/73); new base G (2/4); total (39/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 124 total aligned reads displayed.

CACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTT  >  minE/608945‑609085

cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCa                                                                        >  1:1887477/1‑71 (MQ=255)

cACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCa                                                                        >  1:1403558/1‑71 (MQ=255)

  cgtgcACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTc                                                                      <  1:3689070/68‑1 (MQ=255)

  cgcgCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc                                                                      <  1:2905324/71‑1 (MQ=255)

  cgcgCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATc                                                                      <  1:3810445/71‑1 (MQ=255)

   gcgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCg                                                                     <  1:2264441/71‑1 (MQ=255)

      cACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctg                                                                  <  1:444757/71‑1 (MQ=255)

       aCCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctgt                                                                 <  1:1484372/71‑1 (MQ=255)

         caaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctgtt                                                                >  1:2861963/3‑70 (MQ=255)

         cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCTTCgt                                                                    <  1:203859/66‑1 (MQ=255)

         cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                               <  1:404777/71‑1 (MQ=255)

         cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGAAAGCAATAGCGTAGTCATCGttgttg                                                               <  1:3363081/71‑1 (MQ=255)

         cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGAAAGCAATAGCGTAGTCATCGttgttg                                                               <  1:76445/71‑1 (MQ=255)

          gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                               >  1:2659473/1‑70 (MQ=255)

           aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                                             <  1:2336403/71‑1 (MQ=255)

           aagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                                                               >  1:2354153/1‑69 (MQ=255)

            agaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                                             >  1:3119929/1‑70 (MQ=255)

             gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                                 <  1:1567467/65‑1 (MQ=255)

             gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                                 <  1:909458/65‑1 (MQ=255)

             gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                                           <  1:3471068/71‑1 (MQ=255)

             gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                                           <  1:2700263/71‑1 (MQ=255)

             gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                                           <  1:3833286/71‑1 (MQ=255)

             gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                                           <  1:2244575/71‑1 (MQ=255)

             gaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                                           >  1:2804005/1‑71 (MQ=255)

             gaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt                                                           <  1:3451581/71‑1 (MQ=255)

              aaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTc                                                          <  1:1842783/71‑1 (MQ=255)

               aCTGCATTTGTTTAGCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                                 <  1:1304360/63‑1 (MQ=255)

               aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                                 <  1:2532501/63‑1 (MQ=255)

               aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                                 <  1:1412339/63‑1 (MQ=255)

               aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                                 <  1:7253/63‑1 (MQ=255)

                cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                                        >  1:1524361/1‑71 (MQ=255)

                   cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCgg                                                       <  1:785074/69‑1 (MQ=255)

                   cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCgg                                                       <  1:908338/69‑1 (MQ=255)

                     tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACg                                                    >  1:642692/1‑70 (MQ=255)

                      ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                      >  1:3136309/1‑67 (MQ=255)

                      ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCa                                                  >  1:1980019/1‑71 (MQ=255)

                        gTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                                                 >  1:1297535/1‑70 (MQ=255)

                         tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                                                                  <  1:1340539/52‑1 (MQ=255)

                         tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                                                                  <  1:144171/52‑1 (MQ=255)

                         tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATc                                                <  1:1992040/70‑1 (MQ=255)

                         tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATc                                                <  1:2813371/70‑1 (MQ=255)

                           tACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                                             >  1:1184445/1‑55 (MQ=255)

                            aCCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGt                                            >  1:1691041/1‑71 (MQ=255)

                             ccAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGt                                            <  1:2234179/70‑1 (MQ=255)

                              cAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                                        >  1:1200391/1‑57 (MQ=255)

                                aCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                      >  1:2496742/1‑57 (MQ=255)

                                aCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                      >  1:2558266/1‑57 (MQ=255)

                                aCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                         >  1:1839732/1‑70 (MQ=255)

                                 cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAgg                                             <  1:596952/65‑1 (MQ=255)

                                 cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                       <  1:1121057/71‑1 (MQ=255)

                                 cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                       <  1:3198513/71‑1 (MQ=255)

                                 cGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCg                                       <  1:3905242/71‑1 (MQ=255)

                                  gATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt                                      <  1:3827410/71‑1 (MQ=255)

                                  gATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt                                      <  1:1449555/71‑1 (MQ=255)

                                    tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCa                                                                        >  1:469559/1‑35 (MQ=255)

                                      atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgtt                                                                <  1:3935854/41‑1 (MQ=255)

                                      atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                   >  1:39150/1‑70 (MQ=255)

                                      atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                   >  1:2943661/1‑70 (MQ=255)

                                       tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                                                   <  1:2213414/37‑1 (MQ=255)

                                       tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                                        <  1:1884846/48‑1 (MQ=255)

                                       tcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGGCATCGTTCTCa                                 <  1:1346632/71‑1 (MQ=255)

                                        cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                 >  1:850367/1‑70 (MQ=255)

                                        cGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                 >  1:3156580/1‑70 (MQ=255)

                                            cTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa                               <  1:884056/68‑1 (MQ=255)

                                            cTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                            >  1:1478820/1‑71 (MQ=255)

                                            cTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTg                            <  1:2620261/71‑1 (MQ=255)

                                               aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                      >  1:1964762/1‑42 (MQ=255)

                                                cGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCa                          >  1:477454/1‑69 (MQ=255)

                                                   agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                       >  1:2136143/1‑69 (MQ=255)

                                                   agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                     <  1:2359009/71‑1 (MQ=255)

                                                   agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                     <  1:3280242/71‑1 (MQ=255)

                                                    gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATAGTTCTCATACTGCAGagaag                    <  1:1125767/71‑1 (MQ=255)

                                                     caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAgaga                      <  1:3022563/68‑1 (MQ=255)

                                                     caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaaga                   <  1:3780678/71‑1 (MQ=255)

                                                     caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                     >  1:2598318/1‑69 (MQ=255)

                                                      aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAg                                              >  1:59411/1‑43 (MQ=255)

                                                      aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                   <  1:1411631/54‑1 (MQ=255)

                                                      aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                   <  1:2453730/54‑1 (MQ=255)

                                                      aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                   <  1:3266642/54‑1 (MQ=255)

                                                      aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                  <  1:325690/71‑1 (MQ=255)

                                                      aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                  <  1:2791412/71‑1 (MQ=255)

                                                      aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAACGTTCTCATACTGCAGAGAagag                  <  1:1999685/71‑1 (MQ=255)

                                                       agcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCa                                                  <  1:591344/38‑1 (MQ=255)

                                                       agcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGt                                      <  1:3252434/50‑1 (MQ=255)

                                                       agcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                  <  1:438393/70‑1 (MQ=255)

                                                       agcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAgg                 >  1:2176052/1‑71 (MQ=255)

                                                         caaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                        >  1:2145448/1‑62 (MQ=255)

                                                              gCGTAGTCATGGTTGTTGAAGTCCGGACGCATCAGGTCATCATTCTCATACTGCAGAGAAGATGTGTCGAt          <  1:1644853/71‑1 (MQ=255)

                                                              gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGa           <  1:2296389/70‑1 (MQ=255)

                                                              gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGa           <  1:1201388/70‑1 (MQ=255)

                                                              gCGTAGTCATAGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtctc                                  >  1:2093133/1‑47 (MQ=255)

                                                               cGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGAt          <  1:769049/70‑1 (MQ=255)

                                                               cGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATg         <  1:1365343/71‑1 (MQ=255)

                                                                   gTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGaca     <  1:2081968/71‑1 (MQ=255)

                                                                   gTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGac      <  1:2083711/70‑1 (MQ=255)

                                                                    tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                                    <  1:277149/39‑1 (MQ=255)

                                                                     cATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACt   >  1:2659549/1‑71 (MQ=255)

                                                                     cATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACt   >  1:1924629/1‑71 (MQ=255)

                                                                      aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                <  1:2249136/57‑1 (MQ=255)

                                                                      aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt  >  1:3858771/1‑71 (MQ=255)

CACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTT  >  minE/608945‑609085

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