BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,142	C→T	34.0%	Q704Q (CAG→CAA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,142	0	C	T	34.0%	194.8 / 120.4	142	Q704Q (CAG→CAA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (60/33); new base T (35/13); total (96/46)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

CGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAA  >  Exported/608073‑608212

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                                       >  1:4004942/1‑71 (MQ=255)

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

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

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

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

                            tCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                                             <  1:5622409/53‑1 (MQ=255)

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

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

                              gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                          <  1:2140849/70‑1 (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)

                              gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                          <  1:2369404/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: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)

                                 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)

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

                                      tttCACGGTTCATCACGTTAACTTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca                                 >  1:4915494/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)

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

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

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

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

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

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

                                            ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                            >  1:5186034/1‑70 (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:2808846/1‑35 (MQ=255)

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

                                                 attacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGcc                      <  1:4386120/68‑1 (MQ=9)

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

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

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

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

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

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

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

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

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

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

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

                                                         gaCGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                          >  1:847201/2‑59 (MQ=9)

                                                         gaCGTTGAGGTGTTGACAGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCa                          >  1:2651620/2‑59 (MQ=255)

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

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

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

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

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

                                                            gTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTgg           >  1:1924003/1‑71 (MQ=255)

                                                            gTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTgg           <  1:4275958/71‑1 (MQ=255)

                                                               cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt        >  1:4177110/1‑71 (MQ=255)

                                                                  gtgtTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACccatccag                            <  1:5696337/48‑2 (MQ=255)

                                                                  gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc     <  1:1212158/71‑1 (MQ=255)

                                                                  gTGATGACCACCTTCGATGGATGCTTCGTGGTGGAAGt                                      >  1:4922195/1‑38 (MQ=255)

                                                                    gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg    >  1:5183963/1‑70 (MQ=255)

                                                                    gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg    >  1:2530585/1‑70 (MQ=255)

                                                                    gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg    >  1:1493020/1‑70 (MQ=255)

                                                                     cTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGaa  <  1:12246/71‑1 (MQ=255)

CGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAA  >  Exported/608073‑608212

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