BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,258	C→T	12.5%	P332P (CCG→CCA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,258	0	C	T	12.5%	248.4 / 19.3	96	P332P (CCG→CCA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (50/34); new base T (5/7); total (55/41)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

AGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACC  >  minE/609191‑609327

aGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCgg                                                                     <  1:2476638/70‑1 (MQ=255)

 ggagTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGt                                                                   <  1:2964678/68‑1 (MQ=255)

     tttttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTggg                                                                    >  1:4588029/1‑66 (MQ=255)

     tttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                                                <  1:151642/70‑1 (MQ=255)

      ttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                                                >  1:1515688/1‑69 (MQ=255)

      ttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                                                >  1:2138134/1‑69 (MQ=255)

       tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCg                                                                      >  1:2000489/1‑62 (MQ=255)

       tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGcc                                                               >  1:2712682/1‑69 (MQ=255)

       tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGcc                                                               >  1:2474570/1‑69 (MQ=255)

        ttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGccgg                                                             <  1:1269378/70‑3 (MQ=255)

          ggTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCagag                                                           >  1:1969865/1‑70 (MQ=255)

           gtcaCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGaa                                                         <  1:4104932/68‑1 (MQ=255)

           gTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAgaga                                                          <  1:2512797/70‑1 (MQ=255)

              aCCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAg                                                      >  1:845085/1‑71 (MQ=255)

              aCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGcc                                                               <  1:765651/62‑1 (MQ=255)

              aCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAg                                                      <  1:3508966/71‑1 (MQ=255)

                cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGtt                                                    >  1:4221212/1‑71 (MQ=255)

                cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                                                     >  1:3660695/1‑70 (MQ=255)

                cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                                                     >  1:2647042/1‑70 (MQ=255)

                cAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGccgg                                                             >  1:2796722/1‑60 (MQ=255)

                   ggTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCAt                                                 <  1:1293650/71‑1 (MQ=255)

                     tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                                               <  1:2079290/71‑1 (MQ=255)

                     tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                                               <  1:3889893/71‑1 (MQ=255)

                         aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAtt                                           >  1:1698968/1‑71 (MQ=255)

                         aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAt                                            >  1:624818/1‑70 (MQ=255)

                          ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                          >  1:4124624/1‑71 (MQ=255)

                          ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                          >  1:572776/1‑71 (MQ=255)

                            gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCgg                                                                     <  1:2518763/42‑1 (MQ=255)

                            gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                        <  1:1691527/71‑1 (MQ=255)

                             tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                                       <  1:3899602/55‑1 (MQ=255)

                             tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                                       <  1:4247526/55‑1 (MQ=255)

                             tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                                                       <  1:3307660/55‑1 (MQ=255)

                             tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                        >  1:4205076/1‑70 (MQ=255)

                                aGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAg                                     >  1:4514107/1‑70 (MQ=255)

                                    ccATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTGCGGAGTACg                                 >  1:3484677/1‑70 (MQ=255)

                                     cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                                        <  1:4291754/46‑1 (MQ=255)

                                     cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                                        <  1:2868013/46‑1 (MQ=255)

                                     cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                                        <  1:3358088/46‑1 (MQ=255)

                                      aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                              >  1:4629968/1‑71 (MQ=255)

                                       tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                                                        >  1:2804630/1‑44 (MQ=255)

                                       tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAg                                                      >  1:797551/1‑46 (MQ=255)

                                       tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                             <  1:3820930/71‑1 (MQ=255)

                                       tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                             <  1:1712763/71‑1 (MQ=255)

                                         ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa                            <  1:821195/70‑1 (MQ=255)

                                         ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa                            <  1:3361070/70‑1 (MQ=255)

                                           accGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg                                        >  1:378396/1‑56 (MQ=255)

                                            ccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                          >  1:5616054/1‑53 (MQ=255)

                                             cGATAGATTCGGTTGCCCAGATCGGGTCGCCAgaga                                                          >  1:2719212/1‑36 (MQ=255)

                                             cGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGaa                                                         >  1:4196245/1‑37 (MQ=255)

                                             cGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                        <  1:1657633/70‑1 (MQ=255)

                                             cGATAGATTCGGTTGCCCAGATCGGGTAGCCAGAGaa                                                         >  1:2396641/1‑37 (MQ=255)

                                              gATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATc                                                >  1:503729/1‑45 (MQ=255)

                                               aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTc                                                   <  1:5490046/41‑1 (MQ=255)

                                               aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc                     >  1:4952182/1‑71 (MQ=255)

                                                tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc                                          <  1:4114301/49‑1 (MQ=255)

                                                tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa                       <  1:2988351/68‑1 (MQ=255)

                                                 aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGc                         >  1:1460356/1‑65 (MQ=255)

                                                 aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                    >  1:5282517/1‑70 (MQ=255)

                                                 aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc                     >  1:4756260/1‑69 (MQ=255)

                                                 aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                   >  1:399609/1‑71 (MQ=255)

                                                  gATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                   <  1:4702818/70‑1 (MQ=255)

                                                  acttCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCa                   >  1:5649674/3‑70 (MQ=255)

                                                      cGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACg               >  1:4007674/1‑70 (MQ=255)

                                                           gCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGtt          >  1:1803539/1‑70 (MQ=255)

                                                            cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc        <  1:5186223/71‑1 (MQ=255)

                                                            cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc        <  1:3545860/71‑1 (MQ=255)

                                                             ccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa       >  1:3512965/1‑71 (MQ=255)

                                                              cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt      <  1:3130415/71‑1 (MQ=255)

                                                              cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt      <  1:2266684/71‑1 (MQ=255)

                                                              cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                        <  1:4897229/53‑1 (MQ=255)

                                                              cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                        <  1:3613448/53‑1 (MQ=255)

                                                              cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                        <  1:42153/53‑1 (MQ=255)

                                                               aGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt      >  1:3455217/1‑70 (MQ=255)

                                                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt      >  1:4891136/1‑70 (MQ=255)

                                                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt      >  1:123839/1‑70 (MQ=255)

                                                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt      >  1:3794498/1‑70 (MQ=255)

                                                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt      >  1:2193704/1‑70 (MQ=255)

                                                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt      >  1:3589419/1‑70 (MQ=255)

                                                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTACGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt      >  1:4897799/1‑70 (MQ=255)

                                                                gATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa    >  1:3407185/1‑71 (MQ=255)

                                                                gATCGGGTCGCCAGAGAACAGTTCATCGTATTTCGGAGTACGCAGGAAGCGa                       >  1:190192/1‑52 (MQ=255)

                                                                gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGc                                <  1:2008814/43‑1 (MQ=255)

                                                                gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                      >  1:4145257/1‑53 (MQ=255)

                                                                gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                      >  1:2798119/1‑53 (MQ=255)

                                                                gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                      >  1:5165173/1‑53 (MQ=255)

                                                                 aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGa                                      >  1:5447180/1‑36 (MQ=255)

                                                                 aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGa                                      >  1:440374/1‑36 (MQ=255)

                                                                 aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc   <  1:3607582/71‑1 (MQ=255)

                                                                  tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCa                               >  1:352507/1‑42 (MQ=255)

                                                                  tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa                       >  1:1816350/1‑50 (MQ=255)

                                                                  tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATa                 >  1:3601350/1‑56 (MQ=255)

                                                                  tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc  <  1:3888221/71‑1 (MQ=255)

                                                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                    >  1:5177806/1‑52 (MQ=255)

                                                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                    >  1:5637492/1‑52 (MQ=255)

                                                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                    >  1:2447193/1‑52 (MQ=255)

                                                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt         >  1:722803/1‑63 (MQ=255)

                                                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa       >  1:5626028/1‑65 (MQ=255)

                                                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc  <  1:5527717/70‑1 (MQ=255)

                                                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc  <  1:2468755/70‑1 (MQ=255)

                                                                   cGGGTCGCCAGAGAACAGTTCATCGTATTACGGAGTACGCAGGAAGCg                        >  1:5453611/1‑48 (MQ=255)

AGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACC  >  minE/609191‑609327

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