BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,777	A→G	33.3%	R159R (CGT→CGC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,777	0	A	G	33.3%	129.2 / 92.7	108	R159R (CGT→CGC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (49/23); new base G (19/17); total (68/40)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

Only 100 of 111 total aligned reads displayed.

TAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATT  >  Exported/609709‑609846

tAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACg                                                                     >  1:3230450/1‑71 (MQ=255)

  gTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACg                                                                     >  1:1429881/1‑69 (MQ=255)

   tCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAg                                                                  >  1:5050522/1‑71 (MQ=255)

     aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                >  1:830147/1‑71 (MQ=255)

     aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                >  1:2137121/1‑71 (MQ=255)

     aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                >  1:1910296/1‑71 (MQ=255)

     aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                >  1:3452840/1‑71 (MQ=255)

     aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                >  1:5506485/1‑71 (MQ=255)

     aCCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                >  1:5385352/1‑71 (MQ=255)

       cGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATg                                                              <  1:269486/71‑1 (MQ=255)

        gataatGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAtat                                                             >  1:840490/1‑69 (MQ=255)

         atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc                                                            >  1:417548/1‑71 (MQ=255)

         atgatACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTc                                                            >  1:3351202/1‑71 (MQ=255)

            atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                                          <  1:1230781/70‑1 (MQ=255)

            atACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCg                                                          >  1:5603563/1‑70 (MQ=255)

               cgcccacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAt                                                               >  1:291098/4‑62 (MQ=255)

               cgcccacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCAt                                                               >  1:2696991/4‑62 (MQ=255)

                  ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGTCAACGCAGGATGTCCGGAGTGTa                                                   <  1:4426013/71‑1 (MQ=255)

                  ccacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                   <  1:119365/71‑1 (MQ=255)

                  ccacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                                                   <  1:4890959/71‑1 (MQ=255)

                   cacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaa                                                  <  1:609488/71‑1 (MQ=255)

                   cacggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                   <  1:472761/70‑1 (MQ=255)

                   cacggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaa                                                  <  1:5234997/71‑1 (MQ=255)

                    acggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                 >  1:1256356/1‑71 (MQ=255)

                    acggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                                                   <  1:826443/69‑1 (MQ=255)

                     cggccaTATGCATCTGGCAGACCTTTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                 >  1:2340349/1‑70 (MQ=255)

                     cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                 >  1:1699264/1‑70 (MQ=255)

                     cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAAc                                                >  1:5312540/1‑71 (MQ=255)

                     cggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAAc                                                >  1:2896784/1‑71 (MQ=255)

                      ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGc                                                                    >  1:709788/1‑50 (MQ=255)

                      ggccaTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACg                                               <  1:114108/71‑1 (MQ=255)

                      ggccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAAca                                               <  1:5066423/71‑2 (MQ=255)

                        ccaTAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATc                                             <  1:4805222/71‑1 (MQ=255)

                         caTATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                                            <  1:1392433/71‑1 (MQ=255)

                           tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGAt                                                               <  1:1389149/50‑1 (MQ=255)

                           tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                   <  1:2781513/62‑1 (MQ=255)

                           tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                          <  1:4442888/71‑1 (MQ=255)

                           tatGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                          <  1:1112119/71‑1 (MQ=255)

                            agccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                          >  1:4424970/4‑70 (MQ=255)

                            agccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                          >  1:4223194/4‑70 (MQ=255)

                            agccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaa                                          >  1:2467907/4‑70 (MQ=255)

                             tGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGa                                                                >  1:1136120/1‑47 (MQ=255)

                             tGCATCTGGCAGACCGGGCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                        >  1:3822322/1‑71 (MQ=255)

                             gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTaaa                                                 >  1:35780/3‑62 (MQ=255)

                             gccATCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCaaat                                         <  1:4228793/68‑2 (MQ=255)

                                 tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATcc                                                           <  1:2465098/48‑1 (MQ=16)

                                 tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATcc                                                           <  1:4288415/48‑1 (MQ=16)

                                 tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCggtg                                                       <  1:3678607/52‑3 (MQ=16)

                                 tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                                                   <  1:145809/56‑1 (MQ=255)

                                 tCTGGTAAACCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTa                                                   <  1:685436/56‑1 (MQ=255)

                                 tCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                   <  1:5121823/56‑1 (MQ=255)

                                  cTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGccc                                   >  1:4247352/1‑71 (MQ=255)

                                   tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCa                                                                   >  1:2553216/1‑38 (MQ=255)

                                   tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                   <  1:1846544/54‑1 (MQ=255)

                                   tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCg                                            >  1:3039087/1‑61 (MQ=255)

                                   tGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCt                                  >  1:5075076/1‑71 (MQ=255)

                                   tGGCAGACCGGTCAGAACACCAGATTTACGG‑AACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTg                                 >  1:2475998/1‑71 (MQ=255)

                                      cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                          >  1:3868989/1‑60 (MQ=255)

                                      cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                          >  1:1376146/1‑60 (MQ=255)

                                      cAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGT‑CGGAGTGTAAACGTCGaa                                          >  1:4552518/1‑59 (MQ=255)

                                       aaaCCGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGt                               >  1:1024750/3‑70 (MQ=255)

                                        gACCGGTCAGAACACCAGATTTACGGCAACGCAgg                                                                 <  1:1134571/35‑1 (MQ=255)

                                         aCCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTtgt                            <  1:183171/71‑1 (MQ=255)

                                           cGGTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTtatgg                          >  1:3710304/1‑70 (MQ=255)

                                            ggTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTa                                                   <  1:3069347/45‑1 (MQ=255)

                                             gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                                   >  1:4106300/1‑60 (MQ=12)

                                             gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                                   >  1:4581481/1‑60 (MQ=12)

                                             gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                                   >  1:4575174/1‑60 (MQ=12)

                                             gTCAGCACGCCAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGccc                                   >  1:5571530/1‑60 (MQ=12)

                                             gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                        >  1:1149188/1‑55 (MQ=255)

                                             gTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAac                                         >  1:1333273/1‑54 (MQ=255)

                                              tCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                 >  1:3094220/1‑45 (MQ=255)

                                              tCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTaaa                                                 >  1:2025101/1‑45 (MQ=255)

                                               cAGAACACCAGATTTACGGCAACGCATGATGTCCGGAGTGTAAACGTCg                                            <  1:417847/49‑1 (MQ=255)

                                                aGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGttt                      >  1:4540961/1‑70 (MQ=255)

                                                 gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTTAACGTCGaa                                          >  1:4909086/1‑49 (MQ=255)

                                                 gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                          >  1:1354427/1‑49 (MQ=255)

                                                 gAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                          >  1:2810451/1‑49 (MQ=255)

                                                  aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                        >  1:3661589/1‑50 (MQ=255)

                                                  aaCACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTa                    >  1:1366394/1‑70 (MQ=255)

                                                      ccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGtt                       <  1:1475341/63‑1 (MQ=16)

                                                      ccAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACgcag               <  1:5526446/71‑4 (MQ=255)

                                                       cAGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                  >  1:4448415/1‑67 (MQ=255)

                                                        aGATTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                  >  1:5020873/1‑66 (MQ=255)

                                                        aGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt             >  1:4220959/1‑71 (MQ=255)

                                                         gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                        >  1:425310/1‑43 (MQ=255)

                                                         gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAaca                                        >  1:3553148/1‑43 (MQ=255)

                                                         gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTAtt             >  1:5335648/1‑70 (MQ=255)

                                                         gATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTc            <  1:2987674/71‑1 (MQ=255)

                                                          aTTTACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGt                        <  1:1816085/58‑1 (MQ=12)

                                                           tttAAGGCAACGCAGGATGTCCGGAGTGTAAACGTCGaa                                          >  1:5282951/1‑39 (MQ=255)

                                                             tACGGCAGCGCAGCATATCCGGTGAGTAAACATCAAATACGCCCTGGTTATGGGTTTTACg                  >  1:624521/1‑61 (MQ=16)

                                                             tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                          <  1:5670167/53‑1 (MQ=255)

                                                             tACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                          <  1:718323/53‑1 (MQ=255)

                                                              aCGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg        >  1:4481659/1‑70 (MQ=255)

                                                                 gCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGa                          <  1:257130/49‑1 (MQ=255)

                                                                  cAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGtttt                     >  1:1238475/1‑53 (MQ=255)

                                                                   aaCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGAtt  >  1:3581036/1‑71 (MQ=255)

                                                                    aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg        >  1:2324190/1‑64 (MQ=255)

                                                                    aCGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTg        <  1:2607159/64‑1 (MQ=255)

TAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGCAACGCAGGATGTCCGGAGTGTAAACGTCGAACACGCCCTGGTTGTGAGTTTTACGGTATTCAGTGAAGATT  >  Exported/609709‑609846

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