BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,159	A→C	16.3%	S699A (TCC→GCC) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,159	0	A	C	16.3%	225.8 / 48.0	116	S699A (TCC→GCC) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (42/55); new base C (2/17); total (44/72)
Fisher's exact test for biased strand distribution p-value = 8.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

TTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCC  >  minE/608089‑608229

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)

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

       tCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg                                                                 >  1:136118/1‑71 (MQ=255)

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

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

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

              cagcagcATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGtggtgg                                                           <  1:3562883/69‑1 (MQ=255)

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

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

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

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

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

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

                     aTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc                                                   <  1:3187811/71‑1 (MQ=255)

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

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

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

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

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

                           cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                                             <  1:1082883/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)

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

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

                                  tcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAcca                                      <  1:3408674/71‑1 (MQ=255)

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

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

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

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

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

                                     cGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGaa                                                         >  1:2506238/1‑49 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                           cGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTg                             <  1:668188/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)

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

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

                                               cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTc                          <  1:3159160/70‑1 (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:2343644/70‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                        accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                 >  1:1538295/1‑70 (MQ=255)

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

                                                         ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                      >  1:2649502/1‑64 (MQ=255)

                                                         ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                      >  1:1479166/1‑64 (MQ=255)

                                                         ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc               >  1:2907297/1‑71 (MQ=255)

                                                          caccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcg              >  1:2952013/1‑71 (MQ=255)

                                                            ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcgtc            <  1:2284422/71‑1 (MQ=255)

                                                            ccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtcgtc            <  1:924198/71‑1 (MQ=255)

                                                               tCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc         >  1:3731260/1‑71 (MQ=255)

                                                                cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt          <  1:474139/69‑1 (MQ=255)

                                                                cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt          <  1:2823118/69‑1 (MQ=255)

                                                                cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc         <  1:3724881/70‑1 (MQ=255)

                                                                cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc         <  1:1797071/70‑1 (MQ=255)

                                                                cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt        >  1:3396910/1‑71 (MQ=255)

                                                                 gatggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt       >  1:3933899/1‑71 (MQ=255)

                                                                   tggatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                         <  1:2676588/51‑1 (MQ=255)

                                                                     gatgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                      >  1:2492549/1‑52 (MQ=255)

                                                                      atgCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc  <  1:1021903/71‑1 (MQ=255)

TTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCC  >  minE/608089‑608229

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