BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,012	G→A	21.5%	D414D (GAC→GAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,012	0	G	A	21.5%	258.3 / 61.9	130	D414D (GAC→GAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (40/62); new base A (16/12); total (56/74)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

GCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGA  >  Exported/608948‑609080

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)

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

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

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

      caaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                                                          <  1:5360228/69‑1 (MQ=255)

      caaagaaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGctgttg                                                          <  1:3782391/69‑1 (MQ=255)

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

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

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

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

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

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

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

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

           aaCTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                                                     <  1:1569837/71‑1 (MQ=255)

             cTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGttgttg                                                          <  1:2567882/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)

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

                    tGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCAt                                            >  1:268613/1‑71 (MQ=255)

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

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

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

                        tACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAg                                         <  1:1626106/70‑1 (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:822877/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

                                   atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAg                                         <  1:2042020/59‑1 (MQ=255)

                                   atcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                              <  1:439156/70‑1 (MQ=255)

                                    ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAg                                         <  1:5126323/56‑1 (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:642556/1‑44 (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:3345073/1‑37 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                   aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                 >  1:4328136/1‑35 (MQ=255)

                                                   aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                 >  1:3000644/1‑35 (MQ=255)

                                                   aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGa                                                 >  1:3979835/1‑35 (MQ=255)

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

                                                   aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag             >  1:3078264/1‑71 (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:1461097/71‑1 (MQ=255)

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

                                                      caaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATc                                      >  1:41501/1‑43 (MQ=12)

                                                       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)

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

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

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

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

                                                           gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAg                                         >  1:3662086/1‑35 (MQ=255)

                                                           gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCAt                                    >  1:1150438/1‑40 (MQ=255)

                                                           gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGtt                                >  1:938943/1‑44 (MQ=255)

                                                           gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                               >  1:2581567/1‑45 (MQ=255)

                                                           gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGAt     <  1:2751295/71‑1 (MQ=255)

                                                           gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGa      <  1:940828/70‑1 (MQ=255)

                                                           gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGAt     >  1:2137280/1‑71 (MQ=255)

                                                              tAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGa  <  1:3301044/71‑1 (MQ=255)

                                                               aGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATAc                         <  1:2223316/47‑1 (MQ=255)

                                                               aGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATAc                         <  1:5126632/47‑1 (MQ=255)

                                                                gTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCAt                           <  1:2293534/44‑1 (MQ=255)

GCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGA  >  Exported/608948‑609080

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