BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,766	G→A	13.7%	N496N (AAC→AAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,766	0	G	A	13.7%	322.3 / 34.1	117	N496N (AAC→AAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (51/50); new base A (4/12); total (55/62)
Fisher's exact test for biased strand distribution p-value = 6.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 128 total aligned reads displayed.

AACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCA  >  minE/608697‑608835

aaCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                                                                      >  1:2160874/1‑71 (MQ=255)

aaCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                                                                      >  1:2155921/1‑71 (MQ=255)

  cGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                     >  1:413721/1‑70 (MQ=255)

  cGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTc                                                                    <  1:1147769/71‑1 (MQ=255)

  cGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTc                                                                    <  1:1321775/71‑1 (MQ=255)

    tCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                                                                  >  1:2991471/1‑71 (MQ=255)

     cACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                                                                  <  1:811859/70‑1 (MQ=255)

      aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                                <  1:169534/71‑1 (MQ=255)

      aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                                <  1:3381467/71‑1 (MQ=255)

      aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                                >  1:3393377/1‑71 (MQ=255)

      aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                                <  1:1382315/71‑1 (MQ=255)

      aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagt                                                                 >  1:421674/1‑70 (MQ=255)

      aCGGTAGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                                <  1:3826744/71‑1 (MQ=255)

      aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg                                                                >  1:1016592/1‑69 (MQ=255)

          tCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagt                                                            <  1:3079021/71‑1 (MQ=255)

            gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                                                           <  1:190988/70‑1 (MQ=255)

            gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                                                           <  1:50034/70‑1 (MQ=255)

            gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                          >  1:796502/1‑71 (MQ=255)

            gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                          >  1:1512976/1‑71 (MQ=255)

               cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg                                                                >  1:1562948/1‑60 (MQ=255)

               cAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                                                        <  1:218100/70‑1 (MQ=255)

                   gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGatgatg                                                                       >  1:1242730/1‑51 (MQ=255)

                   gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgca                                                              >  1:320981/1‑60 (MQ=255)

                   gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgca                                                              >  1:1868695/1‑60 (MQ=255)

                     cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt                                                 <  1:2051919/71‑1 (MQ=255)

                      aTCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt                                                 >  1:1522671/1‑70 (MQ=255)

                       tCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGt                                                       <  1:3051200/63‑1 (MQ=255)

                       tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                     <  1:4145358/49‑1 (MQ=255)

                       tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTa                                               <  1:2708923/71‑1 (MQ=255)

                        cAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTag                                              >  1:3120778/1‑71 (MQ=255)

                         agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                                                           <  1:3079110/57‑1 (MQ=255)

                         agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                                                           <  1:1099994/57‑1 (MQ=255)

                         agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGAt                                                         >  1:1827846/1‑59 (MQ=255)

                         agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                                                        <  1:1085544/60‑1 (MQ=255)

                         agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGttt                                                >  1:3886515/1‑68 (MQ=255)

                         agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                             >  1:451347/1‑71 (MQ=255)

                         agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                             >  1:346124/1‑71 (MQ=255)

                         agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                             >  1:935191/1‑71 (MQ=255)

                         agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                             >  1:1732663/1‑71 (MQ=255)

                         agagagGCTTCGTAGCTGTACTTGTCGTGAATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                             >  1:2619530/1‑71 (MQ=255)

                           agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                     <  1:2126694/45‑1 (MQ=255)

                           agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                     <  1:861542/45‑1 (MQ=255)

                           agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                     <  1:1384714/45‑1 (MQ=255)

                           agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTa                                                      >  1:2418589/1‑60 (MQ=255)

                            gagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcca                                           >  1:4136021/1‑70 (MQ=255)

                             agGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTca                                                                   >  1:4096061/1‑45 (MQ=255)

                             agGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                         >  1:2684317/1‑71 (MQ=255)

                              ggCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagcc                                        <  1:3326718/71‑1 (MQ=255)

                               gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCgccagcca                                       <  1:4113818/71‑1 (MQ=255)

                               gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCgccagcca                                       <  1:3823915/71‑1 (MQ=255)

                               gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCgccagcca                                       <  1:1034688/71‑1 (MQ=255)

                                cTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCCGTGATGTACTGTTTAGccagccag                                      >  1:2053782/1‑71 (MQ=255)

                                    gTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCa                                  >  1:2380215/1‑71 (MQ=255)

                                    gTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTc                                    <  1:1198291/69‑1 (MQ=255)

                                     tAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                     <  1:2990876/35‑1 (MQ=255)

                                     tAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCAt                                 <  1:2160025/71‑1 (MQ=255)

                                      aGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACt                                                    >  1:4136330/1‑51 (MQ=255)

                                      aGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCAt                                 >  1:108145/1‑70 (MQ=255)

                                       gatGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGa                               <  1:2609227/69‑1 (MQ=255)

                                       gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                                <  1:1282908/70‑1 (MQ=255)

                                       gCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATg                                <  1:3559966/70‑1 (MQ=255)

                                         tGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAg                             >  1:4140331/1‑71 (MQ=255)

                                          gTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGt                            >  1:1612845/1‑71 (MQ=255)

                                          gTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGt                            >  1:1471044/1‑71 (MQ=255)

                                           tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAg                             <  1:2537923/69‑1 (MQ=255)

                                           tACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTg                           <  1:1098142/71‑1 (MQ=255)

                                           tACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTg                           <  1:1339310/71‑1 (MQ=255)

                                            aCTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTg                           <  1:2965741/70‑1 (MQ=255)

                                             cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACt                                                    <  1:1753592/44‑1 (MQ=255)

                                             cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGa                          <  1:2422392/70‑1 (MQ=255)

                                             cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                         >  1:3235142/1‑71 (MQ=255)

                                             cTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTg                           >  1:134661/1‑69 (MQ=255)

                                             cTTGTCCTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                         >  1:5771/1‑71 (MQ=255)

                                                gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                        <  1:1579295/69‑1 (MQ=255)

                                                gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                        <  1:717764/69‑1 (MQ=255)

                                                gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCa                      >  1:79584/1‑71 (MQ=255)

                                                   gTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATg                    <  1:402324/70‑1 (MQ=255)

                                                    tGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATgcg                  >  1:3656334/1‑71 (MQ=255)

                                                    tGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATg                                <  1:659908/57‑1 (MQ=255)

                                                     gCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAgccagcca                                       >  1:528735/1‑49 (MQ=255)

                                                     gCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAgccagcca                                       >  1:239678/1‑49 (MQ=255)

                                                       aTGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc               <  1:1369906/71‑1 (MQ=255)

                                                          tAGTGGATGATGTTCAGTGCATTGATTTACTGTTTAGCCAGCCAGTCCATGAAg                             >  1:920581/1‑54 (MQ=38)

                                                          tAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATgc                   >  1:218226/1‑64 (MQ=255)

                                                           aGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                        >  1:26146/1‑58 (MQ=255)

                                                           aGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCAATGAAGTGATc                        >  1:3660021/1‑58 (MQ=255)

                                                            gTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                                <  1:561172/49‑1 (MQ=255)

                                                             tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATgcgc                 <  1:1005433/63‑1 (MQ=255)

                                                              ggATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGa                               <  1:2923111/48‑1 (MQ=255)

                                                               gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACtt       >  1:510941/1‑71 (MQ=255)

                                                               gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACtt       <  1:153771/71‑1 (MQ=255)

                                                               gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCAATCACtt       <  1:3595939/71‑1 (MQ=255)

                                                               aatgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtt       >  1:3186688/2‑71 (MQ=255)

                                                                atgatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACtt       <  1:2269558/70‑1 (MQ=255)

                                                                  gatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCAt                     <  1:1056947/54‑1 (MQ=255)

                                                                  gatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACt        <  1:1927204/67‑1 (MQ=255)

                                                                  gatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcat    >  1:1997222/1‑71 (MQ=255)

                                                                    tgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca  <  1:4129407/71‑1 (MQ=255)

                                                                    tgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca  <  1:3868023/71‑1 (MQ=255)

                                                                    tgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca  <  1:2765494/71‑1 (MQ=255)

AACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCA  >  minE/608697‑608835

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