BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,124	C→T	37.3%	V710V (GTG→GTA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,124	0	C	T	37.3%	138.9 / 111.5	118	V710V (GTG→GTA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (41/33); new base T (30/14); total (71/47)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

Only 100 of 131 total aligned reads displayed.

TACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAC  >  Exported/608055‑608191

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

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

 acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                                   >  1:125761/1‑71 (MQ=255)

 acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                                   >  1:378263/1‑71 (MQ=255)

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

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

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

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

  cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACg                                                                  >  1:1900190/1‑71 (MQ=255)

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

     ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                                                >  1:2256425/1‑70 (MQ=255)

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

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

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

             aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa                                                              >  1:5346504/1‑64 (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:4815554/1‑71 (MQ=255)

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

                   ggATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTg                                                 <  1:4759502/71‑1 (MQ=255)

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

                    gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGa                                                >  1:352683/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:3157662/1‑70 (MQ=255)

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

                       atTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccg                                             >  1:4047357/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)

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

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

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

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

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

                             ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                        <  1:1375712/70‑1 (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:1071379/1‑71 (MQ=255)

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

                                         tCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                                                         <  1:574831/41‑1 (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:4758452/40‑1 (MQ=255)

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

                                              tCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                        <  1:5622409/53‑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:2369404/70‑1 (MQ=255)

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

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

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

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

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

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

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

                                                cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg                     <  1:4463507/69‑1 (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)

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

                                                   cagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg                           <  1:3371798/61‑1 (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:1334170/1‑71 (MQ=255)

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

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

                                                        tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc             >  1:4274676/1‑70 (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)

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

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

                                                              ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat       >  1:5186034/1‑70 (MQ=255)

                                                              ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat       >  1:5401046/1‑70 (MQ=255)

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

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

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

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

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

                                                                  catcaCGTTAACGTTCAGGTGCTGACCACCTTCga                                      >  1:387498/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)

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

TACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAC  >  Exported/608055‑608191

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