BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,000	A→G	11.4%	A418A (GCT→GCC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,000	0	A	G	11.4%	281.7 / 22.9	105	A418A (GCT→GCC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (35/58); new base G (6/6); total (41/64)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

TTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGG  >  minE/608932‑609070

ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag                                                                       >  1:527550/1‑70 (MQ=255)

    ccAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAtc                                                                  <  1:2743812/71‑2 (MQ=255)

     cAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAtcg                                                                 <  1:2406376/71‑3 (MQ=255)

      aGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGc                                                                >  1:1523287/1‑71 (MQ=255)

       ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCATCAAGCAATAGCg                                                               >  1:3409611/1‑71 (MQ=255)

       ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                                               >  1:1617089/1‑71 (MQ=255)

         tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                                               <  1:3466562/69‑1 (MQ=255)

          ttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg                                                            >  1:3045171/1‑71 (MQ=255)

           tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAg                                                            <  1:1345190/70‑1 (MQ=255)

           tGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                                           >  1:1450037/1‑71 (MQ=255)

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

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

              aCGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa                                                                    <  1:435288/59‑1 (MQ=255)

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

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

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

                  gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                                    <  1:2275509/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)

                       gaagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGa                                               <  1:2835979/71‑1 (MQ=255)

                        aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagca                                                                     <  1:2171720/48‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                           aaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                                               <  1:2014626/51‑1 (MQ=255)

                           aaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                                               <  1:264776/51‑1 (MQ=255)

                           aaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                                               <  1:2019029/51‑1 (MQ=255)

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

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

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

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

                            aCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttgt                                                  <  1:1412339/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)

                                cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGAc                                      >  1:1919654/1‑71 (MQ=255)

                                cATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGCTGTTGAAGTCCGGAc                                      <  1:796670/71‑1 (MQ=255)

                                  tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                                               <  1:1631342/44‑1 (MQ=255)

                                  tttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                <  1:1793417/59‑1 (MQ=255)

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

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

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

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

                                    tGTTTACCAACGATCATCGGGCTTACGCAGCGAGCAATAGCGTa                                                             >  1:592778/1‑44 (MQ=255)

                                    tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCa                                   <  1:3882253/70‑1 (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)

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

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

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

                                          ccAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                            <  1:2214462/71‑1 (MQ=255)

                                          ccAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                            <  1:534087/71‑1 (MQ=255)

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

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

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

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

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

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

                                               gATCATCGGGCTTACGCAGCAAGCAATAGCGTAGt                                                           <  1:3776659/35‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag        >  1:2136143/1‑69 (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)

                                                                  caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa      >  1:650405/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:2791412/71‑1 (MQ=255)

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

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

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

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

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

TTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGG  >  minE/608932‑609070

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