BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,705	A→G	15.6%	R183R (CGT→CGC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,705	0	A	G	15.6%	263.2 / 33.5	109	R183R (CGT→CGC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (53/39); new base G (8/9); total (61/48)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

AGCCTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGC  >  minE/609636‑609775

aGCCTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGAc                                                                       >  1:686992/1‑71 (MQ=255)

aGCCTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGAc                                                                       >  1:3186792/1‑71 (MQ=255)

   cTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGt                                                                    <  1:927197/71‑1 (MQ=255)

          gaAGTGAACTTTTCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCAcc                                                             <  1:161517/71‑1 (MQ=255)

           aaGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACgg                                                                     <  1:757716/62‑1 (MQ=255)

            aGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACgg                                                                     <  1:1958647/61‑1 (MQ=255)

            aGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCg                                                            <  1:5489663/70‑1 (MQ=255)

               gAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACagtagtaa                                                               >  1:4271934/1‑62 (MQ=255)

               gAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCga                                                           >  1:1597522/1‑68 (MQ=255)

               gAACTCTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACAGTAGTCa                                                               >  1:4278588/1‑64 (MQ=255)

                aaCTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAAGGCGACGGTAGTCACCgat                                                          >  1:4388833/1‑68 (MQ=255)

                aaCTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatga                                                        >  1:5592300/1‑70 (MQ=255)

                 aCTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatga                                                        <  1:4533258/69‑1 (MQ=255)

                  cTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatga                                                        <  1:1501583/68‑1 (MQ=255)

                   tgtgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacg                                                    >  1:227242/1‑71 (MQ=255)

                     tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatga                                                        >  1:5575671/1‑65 (MQ=255)

                     tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatga                                                        >  1:549199/1‑65 (MQ=255)

                     tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggc                                                  >  1:58150/1‑71 (MQ=255)

                     tgCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggc                                                  >  1:274207/1‑71 (MQ=255)

                        cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCAcc                                                             <  1:1285523/57‑1 (MQ=255)

                        cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccac                                               <  1:3335765/71‑1 (MQ=255)

                        cAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccac                                               >  1:4086692/1‑71 (MQ=255)

                            ttGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATTATacggccacggc                                            <  1:939193/70‑1 (MQ=255)

                             tGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCGCCgatgat                                                       >  1:1554669/1‑58 (MQ=255)

                             tGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatgat                                                       >  1:5210927/1‑58 (MQ=255)

                             tGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatgat                                                       >  1:1510659/1‑58 (MQ=255)

                             tGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccac                                               >  1:3453361/1‑66 (MQ=255)

                              gTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTa                                                                   >  1:2004540/1‑45 (MQ=255)

                              gTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATAcggccacggcca                                          <  1:4987662/70‑1 (MQ=255)

                              gTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATAcggccacggcca                                          >  1:4368058/1‑70 (MQ=255)

                              gTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATAcggccacggcca                                          <  1:4981263/70‑1 (MQ=255)

                                cTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatg                                                         >  1:1698368/1‑53 (MQ=255)

                                cTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatg                                                         >  1:4441436/1‑53 (MQ=255)

                                cTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatg                                                         >  1:1497996/1‑53 (MQ=255)

                                cTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCgatg                                                         >  1:1216460/1‑53 (MQ=255)

                                   tcatcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGc                                     <  1:2795615/70‑1 (MQ=255)

                                   tcatcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCa                                    >  1:5208177/1‑71 (MQ=255)

                                     atcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCAt                                   <  1:739018/70‑1 (MQ=255)

                                      tcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCt                                 <  1:2921582/71‑1 (MQ=255)

                                      tcaGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCt                                 <  1:5165410/71‑1 (MQ=255)

                                          gTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATAcggccacggcca                                          <  1:3335217/58‑1 (MQ=255)

                                             gTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                          <  1:1626945/71‑1 (MQ=255)

                                              tCGATACCGTACAGCGCACCGCGACGGTAGTCAcc                                                             <  1:2471051/35‑1 (MQ=255)

                                              tCGATACCGTACAGCGCAACGCGACGGTAGTCAcc                                                             <  1:2487276/35‑1 (MQ=255)

                                              tCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGAc                          <  1:5438625/70‑1 (MQ=255)

                                               tgATGCCATACAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCg                        >  1:4832129/2‑71 (MQ=255)

                                               cGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCg                        >  1:3647090/1‑71 (MQ=255)

                                                gATACCTTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                       >  1:2364034/1‑71 (MQ=255)

                                                gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATAcggccacggcca                                          <  1:2112933/52‑1 (MQ=255)

                                                gATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCgg                       >  1:4304771/1‑71 (MQ=255)

                                                  tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATacggccacggc                                            >  1:3762709/1‑48 (MQ=255)

                                                  tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCAt                                   >  1:3107127/1‑57 (MQ=255)

                                                  tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCAt                                   >  1:5246443/1‑57 (MQ=255)

                                                  tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCAt                                   >  1:2385760/1‑57 (MQ=255)

                                                  tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCAt                                   >  1:867204/1‑57 (MQ=255)

                                                  tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGt                      >  1:1882342/1‑70 (MQ=255)

                                                  tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGt                      >  1:652340/1‑70 (MQ=255)

                                                  tACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGACACGGCCATATGCAt                                   >  1:5125436/1‑57 (MQ=255)

                                                    ccatACAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATGTGGTAAACCGGTCAg                   <  1:3584754/68‑1 (MQ=255)

                                                    ccatACAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAg                   <  1:1652425/68‑1 (MQ=255)

                                                    ccGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCata                                        <  1:3301076/50‑1 (MQ=255)

                                                    ccGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAg                   <  1:819147/71‑1 (MQ=255)

                                                    ccGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAg                   <  1:4882749/71‑1 (MQ=255)

                                                    ccGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAg                   <  1:4061678/71‑1 (MQ=255)

                                                      gTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                 >  1:1245798/1‑71 (MQ=255)

                                                       tACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCa                                    <  1:4261029/51‑1 (MQ=255)

                                                       tACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGaa                 >  1:1790933/1‑70 (MQ=255)

                                                        aCAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCacg               >  1:2452247/1‑70 (MQ=255)

                                                        aCAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCa                                    <  1:3353916/50‑1 (MQ=255)

                                                         cAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTggt                              >  1:1897407/1‑54 (MQ=16)

                                                         cAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTggt                              >  1:2959397/1‑54 (MQ=16)

                                                         cAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTggt                              >  1:2969872/1‑54 (MQ=16)

                                                         cAGCGCTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTggt                              >  1:2963213/1‑54 (MQ=16)

                                                         cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCAt                                   >  1:2296876/1‑50 (MQ=255)

                                                         cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATc                                  >  1:3114443/1‑51 (MQ=255)

                                                         cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAacac              >  1:3290189/1‑71 (MQ=255)

                                                         cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAacac              <  1:2323848/71‑1 (MQ=255)

                                                         cAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAacac              >  1:2925148/1‑71 (MQ=255)

                                                           gcgcTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTggta                             <  1:111725/54‑3 (MQ=12)

                                                            cgcTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTggta                             <  1:4648563/53‑3 (MQ=12)

                                                            cgcTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTggta                             <  1:1037542/53‑3 (MQ=12)

                                                            cgcTACGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTggta                             <  1:1468219/53‑3 (MQ=12)

                                                            cgcAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAacac              >  1:5112888/1‑68 (MQ=255)

                                                             gctaCGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAg           <  1:2540851/67‑1 (MQ=255)

                                                             gctaCGCGGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGa          <  1:3725503/68‑1 (MQ=255)

                                                             gcAACGCGACGGTAGTCACCGATGATAcggccacggcca                                          <  1:4708354/39‑1 (MQ=255)

                                                             gcAACGCGACGGTAGTCACCGATGATAcggccacggcca                                          <  1:2937721/39‑1 (MQ=255)

                                                             gcAACGCGACGGTAGTCACCGATGATAcggccacggcca                                          <  1:4342510/39‑1 (MQ=255)

                                                               aaCGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGAtt        <  1:2322149/71‑1 (MQ=255)

                                                               aaCGCGACGGTAGTCACCGATGATACGGCCAAGGCCATATGCATCTGGAAGAc                          >  1:3947898/1‑53 (MQ=39)

                                                               aaCGCGACGGTAGTCACAGATGATACGGCCAAGGCCATATGCATCTGGCAGAc                          >  1:1878734/1‑53 (MQ=255)

                                                                 cgcgGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCa                    >  1:2939260/1‑57 (MQ=12)

                                                                 cgcgGCGATAGTCACCGATAATGCGCCCACGGCCATAGCCATCTGGTAAACCGGTCAGCACGCCAGATTTa      >  1:1043154/1‑71 (MQ=255)

                                                                 cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa      >  1:1977826/1‑71 (MQ=255)

                                                                 cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa      <  1:3262171/71‑1 (MQ=255)

                                                                 cgcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa      <  1:4904901/71‑1 (MQ=255)

                                                                  gcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTa      <  1:3497576/70‑1 (MQ=255)

                                                                  gcgACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGCTCAGAACACCAGATTTAc     <  1:2408572/71‑1 (MQ=255)

                                                                     aCGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGc  <  1:1593060/71‑1 (MQ=255)

                                                                     aCGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGc  <  1:1586519/71‑1 (MQ=255)

AGCCTGCAGAGAAGTGAACTGTGCCAGTTTGTCTTTCATCAGGTAGTCGATACCGTACAGCGCAACGCGACGGTAGTCACCGATGATACGGCCACGGCCATATGCATCTGGCAGACCGGTCAGAACACCAGATTTACGGC  >  minE/609636‑609775

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`