BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,148	A→G	30.2%	G702G (GGT→GGC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,148	0	A	G	30.2%	195.7 / 96.9	129	G702G (GGT→GGC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (55/35); new base G (28/11); total (83/46)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

Only 100 of 138 total aligned reads displayed.

ATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTT  >  Exported/608078‑608215

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

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

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

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

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

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

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

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

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

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

      ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                                                <  1:1375712/70‑1 (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)

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

                   cgcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                                              >  1:916362/1‑59 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                         cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCt                                                                  >  1:3524660/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)

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

                           gcagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAg                                          >  1:5649424/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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                             ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGccg                        <  1:3477946/69‑2 (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: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)

                                                    aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATcagcccagcc                     >  1:1906388/1‑67 (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:2440528/1‑71 (MQ=255)

                                                      cGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGATCAGCCAGGTTg               >  1:1063503/1‑71 (MQ=255)

                                                         tgaGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGtt            >  1:3347066/3‑70 (MQ=255)

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

                                                            ggTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa         >  1:2654879/1‑71 (MQ=255)

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

                                                             gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa         >  1:2324261/1‑70 (MQ=255)

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

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

                                                              tGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcca                                    >  1:4328242/1‑42 (MQ=255)

                                                               gCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc        <  1:4439662/69‑1 (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)

                                                                 tGACCGCCTTCGACATACGCTTCGTGGTGGAAATACACATCCAGCAGGCCGACAAGGTTGGTTTTACGTAc    >  1:1881386/1‑71 (MQ=255)

                                                                 tGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                          >  1:1664053/1‑49 (MQ=255)

                                                                 tGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATAAGAc                          >  1:3797227/1‑49 (MQ=255)

                                                                  gACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg       <  1:4104449/67‑1 (MQ=255)

                                                                   accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg       >  1:935576/1‑66 (MQ=255)

                                                                   accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg       >  1:3098408/1‑66 (MQ=255)

                                                                   accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt   <  1:1572301/70‑1 (MQ=255)

                                                                   accaccTACGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACtt  >  1:4563970/1‑71 (MQ=255)

ATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTT  >  Exported/608078‑608215

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`