BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,115	A→C	38.5%	R713R (CGT→CGG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,115	0	A	C	38.5%	109.1 / 155.0	130	R713R (CGT→CGG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (50/30); new base C (33/17); total (83/47)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 133 total aligned reads displayed.

CGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCC  >  Exported/608045‑608184

cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCa                                                                       >  1:686846/1‑71 (MQ=255)

   aGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCga                                                                    >  1:611247/1‑70 (MQ=12)

    gCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCgat                                                                   >  1:204410/1‑69 (MQ=12)

       agagACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtc                                                                 <  1:3476824/70‑1 (MQ=9)

       aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                                                >  1:4707475/1‑71 (MQ=255)

        gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                                                               <  1:1208217/68‑1 (MQ=9)

        gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                                                               <  1:3047823/68‑1 (MQ=9)

        gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                                                               <  1:1890052/68‑1 (MQ=9)

          tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                                              <  1:3986411/70‑1 (MQ=255)

          gacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcatta                                                             >  1:2583400/2‑69 (MQ=9)

            cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg                                                           <  1:5515929/71‑1 (MQ=12)

            cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg                                                           <  1:2878476/71‑1 (MQ=12)

            cacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACg                                                           <  1:3994210/71‑1 (MQ=12)

              cGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                                         <  1:724340/71‑1 (MQ=255)

                    gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGt                                                    <  1:623050/70‑1 (MQ=255)

                      cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                                              >  1:3268764/1‑58 (MQ=255)

                       aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                                                >  1:4238926/1‑71 (MQ=255)

                          ttaggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTg                                             <  1:430491/68‑1 (MQ=255)

                            cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                             <  1:3561424/69‑1 (MQ=255)

                            cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                           <  1:2347595/71‑1 (MQ=255)

                            cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                           >  1:4815554/1‑71 (MQ=255)

                             ggATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                 <  1:684597/64‑1 (MQ=255)

                             ggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgtt                                           <  1:1529062/70‑3 (MQ=255)

                              gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGa                                         >  1:5338729/1‑71 (MQ=12)

                               atatTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTg                                             <  1:5669114/66‑1 (MQ=255)

                                tatTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGAc                                                      >  1:28388/1‑56 (MQ=255)

                                tatTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGacc                                       >  1:3784295/1‑71 (MQ=12)

                                 atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                           >  1:223045/1‑66 (MQ=255)

                                 atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                                      >  1:762954/1‑71 (MQ=255)

                                 atTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccg                                      >  1:2994233/1‑70 (MQ=255)

                                 atTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccg                                      >  1:3157662/1‑70 (MQ=255)

                                 atTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccg                                      >  1:4047357/1‑70 (MQ=255)

                                 atTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGacc                                       >  1:4518422/1‑70 (MQ=255)

                                  tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                 >  1:1148539/1‑59 (MQ=255)

                                  tttttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGacc                                       <  1:1020302/69‑1 (MQ=255)

                                   ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCa                                                                       <  1:659991/36‑1 (MQ=255)

                                   ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                                                      >  1:4531187/1‑37 (MQ=255)

                                   ttttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                                    >  1:2690368/1‑71 (MQ=255)

                                    tttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTg                                          >  1:1028818/1‑64 (MQ=255)

                                    tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                                            >  1:4919348/1‑61 (MQ=255)

                                    tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                                            >  1:2252157/1‑61 (MQ=255)

                                    tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                                    <  1:41402/70‑1 (MQ=255)

                                    tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                                    <  1:1711094/70‑1 (MQ=255)

                                    tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                   >  1:4776960/1‑71 (MQ=255)

                                    tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                   >  1:4862636/1‑71 (MQ=255)

                                    tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                   >  1:5486463/1‑71 (MQ=255)

                                    tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                   >  1:3596185/1‑71 (MQ=255)

                                     ttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACg                                                                     >  1:3913732/1‑36 (MQ=255)

                                     ttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                                               >  1:2108252/1‑42 (MQ=255)

                                      tCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTc                                 >  1:4702231/1‑71 (MQ=12)

                                      tCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTc                                 >  1:894824/1‑71 (MQ=12)

                                      tCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTc                                 >  1:2933824/1‑71 (MQ=12)

                                       ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                                         <  1:3789738/46‑1 (MQ=255)

                                       ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                 <  1:1375712/70‑1 (MQ=255)

                                          ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                         >  1:4418540/1‑59 (MQ=255)

                                          ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                         >  1:3606674/1‑59 (MQ=255)

                                              ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                               >  1:1784602/1‑65 (MQ=255)

                                               tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgat                              >  1:1564142/1‑65 (MQ=255)

                                                ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                       >  1:2694238/1‑71 (MQ=255)

                                                ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                       >  1:1309634/1‑71 (MQ=255)

                                                ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                       >  1:1071379/1‑71 (MQ=255)

                                                     gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtg                  >  1:5125089/1‑71 (MQ=255)

                                                     gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgt                   <  1:4659684/70‑1 (MQ=255)

                                                       gTCGAGAAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                <  1:2142944/71‑1 (MQ=255)

                                                        tCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                              <  1:4177272/40‑1 (MQ=255)

                                                        tCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                                              <  1:4758452/40‑1 (MQ=255)

                                                         cGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggtg               >  1:2851578/1‑70 (MQ=255)

                                                         cGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggtg               >  1:1158104/1‑70 (MQ=255)

                                                         cGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggtg               >  1:229330/1‑70 (MQ=255)

                                                          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg              <  1:1019380/70‑1 (MQ=255)

                                                          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg              <  1:4259562/70‑1 (MQ=255)

                                                          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg              <  1:3282078/70‑1 (MQ=255)

                                                          gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg              <  1:2385207/70‑1 (MQ=255)

                                                          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                   >  1:5378634/2‑49 (MQ=12)

                                                          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                   >  1:3524660/2‑49 (MQ=12)

                                                          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCga                               >  1:449414/2‑53 (MQ=12)

                                                          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg              <  1:4463507/69‑1 (MQ=255)

                                                          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg              >  1:5280240/2‑70 (MQ=255)

                                                          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa             <  1:2301901/70‑1 (MQ=12)

                                                            gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa           >  1:2948280/1‑70 (MQ=12)

                                                            gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa           >  1:2978923/1‑70 (MQ=12)

                                                            gcagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                 <  1:2723896/49‑1 (MQ=255)

                                                            gcagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                 <  1:5141500/49‑1 (MQ=255)

                                                            gcagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAg           >  1:4247507/1‑71 (MQ=255)

                                                            gcagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAg           >  1:5649424/1‑71 (MQ=255)

                                                             cagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg                    <  1:1564587/61‑1 (MQ=255)

                                                              agcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATa         >  1:4622590/1‑71 (MQ=12)

                                                               gcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATac        >  1:5666752/1‑70 (MQ=12)

                                                                caTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATacc       <  1:3871174/71‑3 (MQ=9)

                                                                caTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc       <  1:3379660/71‑1 (MQ=255)

                                                                caTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc       >  1:1334170/1‑71 (MQ=255)

                                                                  tttCACGGTTCATCACGTTAACTTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca     >  1:4915494/1‑71 (MQ=255)

                                                                  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGac                                        >  1:3132951/1‑36 (MQ=255)

                                                                  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca     >  1:2504342/1‑71 (MQ=255)

                                                                  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc      >  1:4281342/1‑70 (MQ=255)

                                                                  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAaca     >  1:3366388/1‑71 (MQ=255)

                                                                   ttCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa           >  1:2831214/1‑63 (MQ=255)

                                                                    tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg              <  1:2085760/60‑1 (MQ=255)

                                                                    tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccat    <  1:700995/70‑1 (MQ=255)

                                                                      aCGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc  <  1:2767935/70‑1 (MQ=255)

CGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCC  >  Exported/608045‑608184

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