BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,006	A→C	21.9%	A416A (GCT→GCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,006	0	A	C	21.9%	203.8 / 75.2	128	A416A (GCT→GCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (38/62); new base C (15/13); total (53/75)
Fisher's exact test for biased strand distribution p-value = 1.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

AGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCG  >  minE/608938‑609075

aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                                                     >  1:1211034/1‑71 (MQ=255)

aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                                                     >  1:945354/1‑71 (MQ=255)

aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                                                     >  1:2480333/1‑71 (MQ=255)

aGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGc                                                                     <  1:4923163/71‑1 (MQ=255)

 ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                                                    <  1:224252/71‑1 (MQ=255)

  gTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt                                                                   <  1:2925223/71‑1 (MQ=255)

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

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

           cgcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgt                                                          <  1:1980355/71‑1 (MQ=255)

            gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                                         >  1:3664825/1‑71 (MQ=255)

            gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                                         >  1:4528867/1‑71 (MQ=255)

            gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgtt                                                         >  1:1677563/1‑71 (MQ=255)

            gcACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgt                                                          >  1:140475/1‑70 (MQ=255)

             cACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                                                     <  1:3759162/58‑1 (MQ=255)

             cACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                                                     <  1:4355125/58‑1 (MQ=255)

             cACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCgttg                                                        <  1:2154030/71‑1 (MQ=255)

             cACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCgctg                                                        >  1:3185884/1‑71 (MQ=255)

                cGAAGATCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                     <  1:2340665/71‑1 (MQ=255)

                cGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCg                                                                    <  1:192719/56‑1 (MQ=255)

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

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

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

                 aaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                                                     >  1:633249/2‑70 (MQ=255)

                 aaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                                                     >  1:1616873/2‑70 (MQ=255)

                  aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                                                       >  1:4163056/1‑51 (MQ=255)

                  aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATa                                                                       >  1:1469284/1‑51 (MQ=255)

                  aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc                                                                     <  1:3781532/53‑1 (MQ=255)

                  aagaaCTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                     >  1:4654861/1‑69 (MQ=255)

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

                    gaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt                                                 >  1:24494/1‑71 (MQ=255)

                     aaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcgc                                                          <  1:3580464/61‑2 (MQ=255)

                       cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                     <  1:2567882/64‑1 (MQ=255)

                       cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                     <  1:302389/64‑1 (MQ=255)

                       cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGaa                                                   >  1:5081927/1‑66 (MQ=255)

                        tGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTc                                                               <  1:762743/53‑1 (MQ=255)

                           aTTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACg                                          <  1:1878427/71‑1 (MQ=255)

                             ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTc                                                               <  1:3881398/48‑1 (MQ=255)

                             ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                              <  1:659474/65‑1 (MQ=255)

                             ttGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCg                                              <  1:5176578/65‑1 (MQ=255)

                               gTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                            >  1:5244188/1‑65 (MQ=255)

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

                                tttACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                                       >  1:4548540/1‑69 (MQ=255)

                                   aCCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAgg                                   <  1:3253912/70‑1 (MQ=255)

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

                                      aaCGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                               <  1:2096741/71‑1 (MQ=255)

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

                                          accattGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtt                           <  1:2958484/69‑1 (MQ=255)

                                           tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                 >  1:4006540/1‑64 (MQ=255)

                                           tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                 >  1:947131/1‑64 (MQ=255)

                                           tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                 >  1:1007848/1‑64 (MQ=255)

                                           tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                 >  1:1557215/1‑64 (MQ=255)

                                           tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                          <  1:3016969/71‑1 (MQ=255)

                                           tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                          <  1:2113837/71‑1 (MQ=255)

                                            catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                         >  1:5542797/1‑71 (MQ=255)

                                            catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                         >  1:4750025/1‑71 (MQ=255)

                                              ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAg                                    <  1:5126323/56‑1 (MQ=255)

                                              ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                       >  1:3290511/3‑71 (MQ=255)

                                              ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                       >  1:5684195/3‑71 (MQ=255)

                                              ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                       >  1:4556021/3‑71 (MQ=255)

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

                                                gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCAt                      <  1:4952891/70‑1 (MQ=255)

                                                gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa                     >  1:230515/1‑71 (MQ=255)

                                                gggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTca                                               >  1:1596959/1‑44 (MQ=255)

                                                    ttACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                  <  1:2953341/70‑1 (MQ=255)

                                                     tACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                  >  1:683461/1‑69 (MQ=255)

                                                      aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                                 <  1:3564746/37‑1 (MQ=255)

                                                      aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGt                                                 >  1:3345073/1‑37 (MQ=255)

                                                      aCGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                  >  1:4535093/1‑68 (MQ=255)

                                                       cGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga              <  1:2738835/71‑1 (MQ=255)

                                                        gcagcaGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCaagg             <  1:3895082/71‑4 (MQ=255)

                                                        gcagcaGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCaagg             <  1:4350307/71‑4 (MQ=255)

                                                        gcagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTc                                 <  1:373776/51‑1 (MQ=255)

                                                         cagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAgaga            <  1:2118396/71‑1 (MQ=255)

                                                          agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGa                                            >  1:1982035/1‑38 (MQ=255)

                                                          agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGATCTCATACTGCAGagaa           <  1:1570117/71‑1 (MQ=255)

                                                           gcaggcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAggaag          <  1:2366258/71‑1 (MQ=255)

                                                           gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                       <  1:564540/58‑1 (MQ=255)

                                                           gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                       <  1:1710426/58‑1 (MQ=255)

                                                           gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                       <  1:553177/58‑1 (MQ=255)

                                                           gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag          <  1:2007039/71‑1 (MQ=255)

                                                            caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga              <  1:154232/66‑1 (MQ=255)

                                                            caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga              <  1:312338/66‑1 (MQ=255)

                                                             aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa           >  1:515667/1‑68 (MQ=255)

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

                                                             aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag        >  1:3078264/1‑71 (MQ=255)

                                                              ggcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                       >  1:905889/2‑55 (MQ=255)

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

                                                               gcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgt      <  1:4835971/71‑1 (MQ=255)

                                                               gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt      <  1:1435471/71‑1 (MQ=255)

                                                               gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt      >  1:4337394/1‑71 (MQ=255)

                                                               gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt      <  1:419717/71‑1 (MQ=255)

                                                               gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt      <  1:1461097/71‑1 (MQ=255)

                                                                 aaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgtga    <  1:4618897/71‑2 (MQ=255)

                                                                  aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCa                                     <  1:2824134/37‑1 (MQ=255)

                                                                  aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCa                                <  1:4323708/42‑1 (MQ=255)

                                                                  aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc   <  1:1082422/71‑1 (MQ=255)

                                                                  aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc   <  1:1551045/71‑1 (MQ=255)

                                                                   tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                       <  1:1385924/50‑1 (MQ=255)

                                                                   tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga              <  1:1820379/59‑1 (MQ=255)

                                                                   tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCg  <  1:5383566/71‑1 (MQ=255)

AGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCG  >  minE/608938‑609075

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