BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,142	0	C	T	16.8%	288.5 / 41.8	119	Q704Q (CAG→CAA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (46/53); new base T (3/17); total (49/70)
Fisher's exact test for biased strand distribution p-value = 1.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

GCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACG  >  minE/608072‑608210

gCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGc                                                                      >  1:3690568/1‑71 (MQ=255)

   gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                                                   <  1:1929460/71‑1 (MQ=255)

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

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

       tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccac                                                               <  1:2149654/71‑1 (MQ=255)

        ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                                              >  1:3251582/1‑71 (MQ=255)

        ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                                              >  1:855116/1‑71 (MQ=255)

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

             cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                                         <  1:868636/71‑1 (MQ=255)

             cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                                         <  1:1278827/71‑1 (MQ=255)

              gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgat                                                        <  1:1627687/71‑1 (MQ=255)

                gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                                       >  1:1553947/1‑70 (MQ=255)

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

                   ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat                                                    >  1:1977977/1‑70 (MQ=255)

                    tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc                                                  >  1:3689004/1‑71 (MQ=255)

                     ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTGAGGTGCTGACCACCTTCGATGGATGc                                                  <  1:3419300/70‑1 (MQ=255)

                     ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc                                                  <  1:1213134/70‑1 (MQ=255)

                      cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCtt                                                >  1:3171065/1‑71 (MQ=255)

                       aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCtt                                                            >  1:1577473/1‑58 (MQ=255)

                       aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTc                                               >  1:350886/1‑71 (MQ=255)

                            gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                                                   >  1:3597212/1‑46 (MQ=255)

                            gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                                          >  1:2158204/1‑71 (MQ=255)

                             tCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                                         <  1:2830291/55‑1 (MQ=255)

                               gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                                                     >  1:3094586/1‑41 (MQ=255)

                               gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                        >  1:1551415/1‑70 (MQ=255)

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

                               cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa                                       <  1:1307701/70‑1 (MQ=255)

                               cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGa                                       <  1:1099635/70‑1 (MQ=255)

                                agcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                                                              <  1:1502184/47‑1 (MQ=255)

                                agcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaa                                      <  1:2845961/71‑1 (MQ=12)

                                 gcagcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                                     <  1:2209198/71‑2 (MQ=255)

                                  cagcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGt                                    <  1:410139/71‑1 (MQ=255)

                                    gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa                                  >  1:686513/1‑71 (MQ=255)

                                    gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa                                  >  1:690240/1‑71 (MQ=255)

                                       tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                                          <  1:3859320/60‑1 (MQ=255)

                                        ttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccat                              >  1:1077683/1‑71 (MQ=255)

                                         tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat                                                    <  1:2052940/48‑1 (MQ=255)

                                          cACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                                          >  1:923323/1‑57 (MQ=255)

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

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

                                            cGGTTCATCACGTTAACGTTCAGGTGCTTACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                          <  1:3467065/71‑1 (MQ=255)

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

                                              gTTCATCACGTTAACGTTCAGGTTCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                        <  1:1865490/71‑1 (MQ=255)

                                              gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                          >  1:2550594/1‑69 (MQ=255)

                                              gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                        <  1:743175/71‑1 (MQ=255)

                                              gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                        <  1:1093349/71‑1 (MQ=255)

                                              gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                        <  1:3007658/71‑1 (MQ=255)

                                               ttCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                        <  1:400357/70‑1 (MQ=255)

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

                                                 catcaCGTTAACGATCAGGTGCTGACCACCTACgatgga                                                     >  1:2430576/1‑39 (MQ=37)

                                                  atcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                     <  1:1054824/70‑1 (MQ=255)

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

                                                    tacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaa                                      >  1:556602/2‑51 (MQ=255)

                                                     aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagc                 >  1:1065758/1‑71 (MQ=255)

                                                     aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagc                 >  1:2987300/1‑71 (MQ=255)

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

                                                     aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag                  >  1:462972/1‑70 (MQ=255)

                                                       gTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACccatccag                          <  1:3062689/60‑2 (MQ=255)

                                                       gTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACccatccag                          <  1:2263703/60‑2 (MQ=255)

                                                       gTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagcca               <  1:10364/71‑1 (MQ=255)

                                                       gTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagcca               <  1:1102502/71‑1 (MQ=255)

                                                        ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGt                                    <  1:2996603/49‑1 (MQ=255)

                                                        ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa                                  >  1:1356578/1‑51 (MQ=255)

                                                        ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                     >  1:3383687/1‑64 (MQ=255)

                                                         tAACGTTCAGGTGCTGACCACCTTCGTTGGATGCtt                                                <  1:2838288/36‑1 (MQ=255)

                                                         tAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAc                                 >  1:448236/1‑51 (MQ=255)

                                                          aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatcc                            >  1:1893065/1‑55 (MQ=255)

                                                          aaCGTTCAGGGGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTa                                   <  1:3781202/48‑1 (MQ=255)

                                                           aCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt           <  1:3703610/71‑1 (MQ=255)

                                                           aCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt           >  1:1799271/1‑71 (MQ=255)

                                                            cGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTg          <  1:3223304/71‑1 (MQ=255)

                                                            cGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTg          <  1:668188/71‑1 (MQ=255)

                                                            cGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTg          <  1:2821012/71‑1 (MQ=255)

                                                            cGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTg          <  1:2482210/71‑1 (MQ=255)

                                                             gTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt           >  1:3239501/1‑69 (MQ=255)

                                                              ttCAGGTGCTGACCACCTTCGATGGATGCTTCgtgg                                           <  1:2320480/36‑1 (MQ=255)

                                                              ttCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                          <  1:1175710/53‑1 (MQ=255)

                                                              ttCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                     <  1:3874000/58‑1 (MQ=255)

                                                                gaGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGttt      >  1:2508660/2‑69 (MQ=255)

                                                                cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTc       <  1:3159160/70‑1 (MQ=255)

                                                                cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTc       >  1:1792913/1‑70 (MQ=255)

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

                                                                cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGT‑GAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt      >  1:2402537/1‑70 (MQ=255)

                                                                 aGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGtttt     <  1:942347/71‑1 (MQ=255)

                                                                 aGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGtttt     <  1:2231230/71‑1 (MQ=255)

                                                                 aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa                                  >  1:2491314/1‑42 (MQ=255)

                                                                 aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt      <  1:2804426/70‑1 (MQ=255)

                                                                 aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt      <  1:2700611/70‑1 (MQ=255)

                                                                 aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt      <  1:2689224/70‑1 (MQ=255)

                                                                 aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt      <  1:2362760/70‑1 (MQ=255)

                                                                 aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt      <  1:180736/70‑1 (MQ=255)

                                                                 aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAACCTGGCTGGTCt      <  1:3119832/70‑1 (MQ=255)

                                                                  ggTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTa    <  1:3824925/71‑1 (MQ=255)

                                                                  ggTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTa    <  1:3409193/71‑1 (MQ=255)

                                                                  ggTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACAATCCATCAGACCAGCCAGGTTGGTCTTa    <  1:3880655/71‑1 (MQ=255)

                                                                   gtgtTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACccatcca                           <  1:2149349/47‑1 (MQ=255)

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

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

                                                                   gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc   >  1:923689/1‑71 (MQ=255)

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

GCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACG  >  minE/608072‑608210

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