BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,130	A→G	34.1%	V708V (GTT→GTC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,130	0	A	G	34.1%	158.6 / 111.3	126	V708V (GTT→GTC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (48/35); new base G (33/10); total (81/45)
Fisher's exact test for biased strand distribution p-value = 4.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 143 total aligned reads displayed.

GATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGT  >  Exported/608061‑608200

    gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa                                                                       <  1:388756/67‑1 (MQ=255)

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

      cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                                   >  1:672858/1‑69 (MQ=255)

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

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

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

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

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

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

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

                 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)

                  tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccac                                                     >  1:4866991/1‑71 (MQ=255)

                  tttttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGacc                                                       <  1:1020302/69‑1 (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)

                    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)

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

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

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

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

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

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

                            ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                           >  1:4004942/1‑71 (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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                          cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                                   >  1:5378634/2‑49 (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:2978923/1‑70 (MQ=12)

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

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

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

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

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

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

                                                 aTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATacc                       >  1:1072629/1‑68 (MQ=9)

                                                  tttCACGGTTCATCACGTTAACGTTCAGGTGCTGac                                                        >  1:3132951/1‑36 (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)

                                                    tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatc                   >  1:701406/1‑71 (MQ=255)

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

                                                       cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                <  1:3082681/71‑1 (MQ=255)

                                                         attCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa              <  1:4609794/70‑1 (MQ=12)

                                                            cattaCGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCcagcaggc           >  1:5406626/1‑69 (MQ=255)

                                                            cattaCGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCcagcaggc           >  1:2194135/1‑69 (MQ=255)

                                                            catcaCGTTAACGTTCAGGTGCTGACCACCTTCga                                               >  1:5383557/1‑35 (MQ=255)

                                                            catcaCGTTAACGTTCAGGTGCTGACCACCTTCga                                               >  1:387498/1‑35 (MQ=255)

                                                            catcaCGTTAACGTTCAGGTGCTGACCACCTTCga                                               >  1:2808846/1‑35 (MQ=255)

                                                            catcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc                       <  1:1334738/59‑1 (MQ=255)

                                                            catcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc           >  1:5157747/1‑71 (MQ=255)

                                                            catcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc           >  1:4958019/1‑71 (MQ=255)

                                                              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa              <  1:2189330/64‑1 (MQ=9)

                                                              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa              <  1:5129726/64‑1 (MQ=9)

                                                              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGccg         <  1:410356/69‑2 (MQ=255)

                                                              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGccg         <  1:5048742/69‑2 (MQ=255)

                                                              ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGccg         <  1:3477946/69‑2 (MQ=255)

                                                              tcaCGTTAACGTTCAGGTCCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAg             <  1:244548/67‑1 (MQ=255)

                                                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGaa                            >  1:3474597/1‑50 (MQ=255)

                                                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGaa                            >  1:764585/1‑50 (MQ=255)

                                                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                 >  1:3737946/1‑61 (MQ=255)

                                                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                 >  1:1947077/1‑61 (MQ=255)

                                                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcca                 >  1:1697490/1‑61 (MQ=255)

                                                                aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag        <  1:1231279/70‑1 (MQ=255)

                                                                 cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                              <  1:4196643/47‑1 (MQ=255)

                                                                 cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagc       >  1:1497237/1‑70 (MQ=255)

                                                                 cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccaga       >  1:3046544/1‑69 (MQ=255)

                                                                   ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagccag    >  1:5154923/1‑71 (MQ=255)

                                                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATcagcccagcc      >  1:1906388/1‑67 (MQ=255)

                                                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAcc          >  1:1661064/1‑63 (MQ=255)

                                                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAgg   <  1:2450950/70‑1 (MQ=255)

                                                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt  >  1:5228483/1‑71 (MQ=255)

                                                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt  >  1:3487062/1‑71 (MQ=255)

                                                                     aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt  >  1:4374389/1‑71 (MQ=255)

GATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGT  >  Exported/608061‑608200

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