BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,175	G→A	13.9%	Y693Y (TAC→TAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,175	0	G	A	13.9%	299.2 / 31.9	115	Y693Y (TAC→TAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (41/58); new base A (3/13); total (44/71)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

Only 100 of 116 total aligned reads displayed.

CAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAAC  >  minE/608106‑608243

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

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

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

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

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

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

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

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

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

                  caCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAg                                                                      >  1:3902437/1‑52 (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:462972/1‑70 (MQ=255)

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

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

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

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

                     gTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagcca                                                <  1:10364/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)

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

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

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

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

                         aCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGtt                                            >  1:1799271/1‑71 (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)

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

                            ttCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                                                      <  1:3874000/58‑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:942347/71‑1 (MQ=255)

                               aGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGTTGGtttt                                      <  1:2231230/71‑1 (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:2700611/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:2149349/47‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                        ccaccTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc                                    >  1:2649502/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)

                                               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)

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

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

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

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

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

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

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

                                                        cTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa                                     <  1:151775/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)

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

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

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

CAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAAC  >  minE/608106‑608243

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