BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,193	A→G	11.6%	A687A (GCT→GCC) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,193	0	A	G	11.6%	252.6 / 20.7	95	A687A (GCT→GCC) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (38/46); new base G (3/8); total (41/54)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 106 total aligned reads displayed.

TCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGG  >  minE/608123‑608260

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

            tCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTc                                                         >  1:1994647/1‑71 (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:2231230/71‑1 (MQ=255)

              aGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGtttt                                                       <  1:942347/71‑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:2343644/70‑1 (MQ=255)

              aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt                                                        <  1:2804426/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)

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

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

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

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

                      accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAAc                                                 <  1:3631156/69‑1 (MQ=255)

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

                      accaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                                                >  1:1538295/1‑70 (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:924198/71‑1 (MQ=255)

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

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

                              cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAGGAACTTcgtcgt                                         <  1:2293151/69‑1 (MQ=255)

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

                              cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                                         <  1:474139/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)

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

                              cGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt                                       >  1:1495472/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)

                                       cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                                      <  1:334105/47‑1 (MQ=255)

                                       cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                                         <  1:463558/60‑1 (MQ=255)

                                       cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGt                              <  1:2497312/71‑1 (MQ=255)

                                       cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGt                              <  1:3467922/71‑1 (MQ=255)

                                       cTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACTGGATCTTCTTTGCCCAgc                              <  1:504485/71‑2 (MQ=255)

                                         tCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGTTTTACGTACtgg                                              >  1:3358305/1‑51 (MQ=255)

                                          cGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGaa                                                  >  1:2059806/1‑48 (MQ=255)

                                           gtggtgGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGt                          >  1:1672782/1‑71 (MQ=255)

                                            ttgtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACt                                                <  1:824844/46‑1 (MQ=255)

                                              gtggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTc                        <  1:546230/70‑1 (MQ=255)

                                               tggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                      <  1:2630236/71‑1 (MQ=255)

                                                ggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCtt                                      <  1:2843820/54‑1 (MQ=255)

                                                ggAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCgg                      >  1:677342/1‑70 (MQ=255)

                                                 gAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAccc                                 >  1:637500/1‑58 (MQ=255)

                                                  aaGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAAc                   >  1:571959/1‑71 (MQ=255)

                                                   aGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAAc                   >  1:2924643/1‑70 (MQ=255)

                                                     tAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                              >  1:3751140/1‑41 (MQ=255)

                                                     tAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTtc                                              >  1:3342509/1‑41 (MQ=255)

                                                     tAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTcgtcgt                                         >  1:622982/1‑46 (MQ=255)

                                                     tAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTAc                                   >  1:1036054/1‑52 (MQ=255)

                                                      aaCCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCgtcgtc                                        >  1:856909/1‑46 (MQ=255)

                                                      aaCCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCt                                       >  1:2703262/1‑47 (MQ=255)

                                                        ccatccatCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCg                           <  1:2832041/57‑1 (MQ=255)

                                                        ccatccatCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCg                           <  1:1975745/57‑1 (MQ=255)

                                                        ccatccatCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATa               >  1:1611263/1‑69 (MQ=255)

                                                        ccatccatCAGACCACCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCg                           <  1:1263972/57‑1 (MQ=255)

                                                         catccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag            <  1:1141675/71‑1 (MQ=255)

                                                           tccatcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATagag            >  1:203586/1‑69 (MQ=255)

                                                             catcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATa               >  1:3905764/1‑64 (MQ=255)

                                                              atcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTg                             <  1:126335/49‑1 (MQ=255)

                                                               tcAGACCAGCCAGGTTGGTCTTACGAACTtcgtcg                                          <  1:2829526/35‑1 (MQ=255)

                                                               tcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAGCGATAGAGAAGgtg      >  1:1607805/1‑71 (MQ=255)

                                                               tcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATa               <  1:1007353/62‑1 (MQ=255)

                                                               tcAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg      <  1:1191281/71‑1 (MQ=255)

                                                                cAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTc                        >  1:2140048/1‑52 (MQ=255)

                                                                cAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg      >  1:3031071/1‑70 (MQ=255)

                                                                 aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtgt     <  1:3074812/70‑1 (MQ=255)

                                                                 aGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGgtg      >  1:151327/1‑69 (MQ=255)

                                                                   aCCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAgg  <  1:2461678/71‑1 (MQ=255)

                                                                    ccagccagGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGtt                         >  1:1256042/1‑47 (MQ=255)

TCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGG  >  minE/608123‑608260

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