BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	609,063	A→G	14.1%	S397S (TCT→TCC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	609,063	0	A	G	14.1%	203.8 / 22.2	85	S397S (TCT→TCC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (33/40); new base G (6/6); total (39/46)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 111 total aligned reads displayed.

CGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGA  >  minE/608993‑609132

cGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                                       <  1:2738835/71‑1 (MQ=255)

  cagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAgaga                                                                     <  1:2118396/71‑1 (MQ=255)

   agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGATCTCATACTGCAGagaa                                                                    <  1:1570117/71‑1 (MQ=255)

    gcaggcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAggaag                                                                   <  1:2366258/71‑1 (MQ=255)

    gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                                                   <  1:2007039/71‑1 (MQ=255)

     caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                                       <  1:154232/66‑1 (MQ=255)

     caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                                       <  1:312338/66‑1 (MQ=255)

      aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                                                    >  1:515667/1‑68 (MQ=255)

      aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                                                    >  1:1216229/1‑68 (MQ=255)

      aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                                                    >  1:3710429/1‑68 (MQ=255)

      aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                                                 <  1:316184/71‑1 (MQ=255)

      aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                                                 >  1:3078264/1‑71 (MQ=255)

       agcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                                                 <  1:3330683/70‑1 (MQ=255)

        gcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgt                                                               <  1:4835971/71‑1 (MQ=255)

        gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                                               <  1:1461097/71‑1 (MQ=255)

        gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                                               <  1:419717/71‑1 (MQ=255)

        gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                                               >  1:4337394/1‑71 (MQ=255)

          aaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgtga                                                             <  1:4618897/71‑2 (MQ=255)

           aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                            <  1:1082422/71‑1 (MQ=255)

           aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                            >  1:2302400/1‑71 (MQ=255)

           aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                            <  1:5588921/71‑1 (MQ=255)

           aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                            <  1:1551045/71‑1 (MQ=255)

            tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                                       <  1:1820379/59‑1 (MQ=255)

            tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                                       <  1:1803017/59‑1 (MQ=255)

            tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCg                                                           <  1:5383566/71‑1 (MQ=255)

              gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGAt                                                         <  1:2751295/71‑1 (MQ=255)

              gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGa                                                          <  1:940828/70‑1 (MQ=255)

              gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGAt                                                         >  1:2137280/1‑71 (MQ=255)

                    tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGacac                                                   <  1:1115463/68‑1 (MQ=255)

                    tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                                      >  1:2152665/4‑68 (MQ=255)

                    tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                                      >  1:1970741/4‑68 (MQ=255)

                    tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                                      >  1:3663707/4‑68 (MQ=255)

                    tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                                                      <  1:5667610/52‑1 (MQ=255)

                      aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                 >  1:705335/1‑71 (MQ=255)

                      aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                 >  1:3179307/1‑71 (MQ=255)

                      aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                 >  1:3379225/1‑71 (MQ=255)

                       tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                                                   <  1:2456467/52‑1 (MQ=255)

                       tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                                                   <  1:1334793/52‑1 (MQ=255)

                       tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                                                   <  1:5640415/52‑1 (MQ=255)

                       tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                 >  1:3110437/1‑70 (MQ=255)

                       tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                 >  1:3116851/1‑70 (MQ=255)

                        cGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                                                   <  1:2280545/51‑1 (MQ=255)

                        cGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                 >  1:2781174/1‑69 (MQ=255)

                         gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa                                               <  1:845802/70‑1 (MQ=255)

                         gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAg                                              >  1:5548385/1‑71 (MQ=255)

                           tgttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                                            >  1:456444/1‑71 (MQ=255)

                             ttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaaga                                                                  <  1:5177378/47‑1 (MQ=255)

                              tgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                                                      <  1:4411910/42‑1 (MQ=255)

                              tgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                                                      <  1:2901027/42‑1 (MQ=255)

                               gAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAgg                                                                >  1:2996215/1‑47 (MQ=255)

                                aaGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                                                    >  1:2942945/1‑42 (MQ=255)

                                    ccGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAAt                                      <  1:1150458/68‑1 (MQ=255)

                                     cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                                   <  1:5337889/70‑1 (MQ=255)

                                     cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                                   <  1:2198931/70‑1 (MQ=255)

                                       gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                <  1:3077275/71‑1 (MQ=255)

                                       gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                <  1:5231204/71‑1 (MQ=255)

                                        aCGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa                               <  1:3909367/71‑1 (MQ=255)

                                         cGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                >  1:1036772/1‑69 (MQ=255)

                                         cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGacac                                                   >  1:1559307/1‑50 (MQ=12)

                                           cATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                                            >  1:1500114/1‑55 (MQ=255)

                                           cATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTGAAgg                            <  1:1224327/71‑2 (MQ=255)

                                            aTCAGATCGTCGTTTTCATACTGCAAAGAGGATGTGACGATAGACACCTGCGCGGCATATTTTTTGaa                              >  1:3717520/1‑68 (MQ=255)

                                             tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                            <  1:16782/37‑1 (MQ=255)

                                             tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                            <  1:611434/37‑1 (MQ=255)

                                             tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                          <  1:2073074/71‑1 (MQ=255)

                                             tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                          <  1:4503201/71‑1 (MQ=255)

                                             tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                          <  1:4737692/71‑1 (MQ=255)

                                              cAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGtt                           >  1:4818782/1‑69 (MQ=255)

                                               aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                            <  1:3710638/35‑1 (MQ=255)

                                               aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGa                                                          >  1:5544640/1‑37 (MQ=255)

                                               aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                              >  1:1385118/1‑65 (MQ=255)

                                               aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                              >  1:2111126/1‑65 (MQ=255)

                                               aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa                               >  1:2873580/1‑64 (MQ=255)

                                                  tCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg                     <  1:4179864/71‑1 (MQ=255)

                                                   cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                 >  1:3452629/1‑42 (MQ=255)

                                                   cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg                     >  1:5694392/1‑70 (MQ=255)

                                                    aTCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGc                    <  1:5530304/70‑1 (MQ=255)

                                                     tCGTTCTCATTCTGCAGAGAAGAGGTGTCGATGGACACTTTAGGGGCGAATttgttg                                <  1:3695301/57‑1 (MQ=255)

                                                     tCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAg                  <  1:1416829/71‑1 (MQ=255)

                                                     tCGTTCTCATACTGCAGAGAAGAGCTGTCGATGGa                                                      <  1:4657500/35‑1 (MQ=255)

                                                       gTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGGGGCAGtt                >  1:145938/1‑71 (MQ=255)

                                                       gTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtt                >  1:3218500/1‑71 (MQ=255)

                                                        ttCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttt               <  1:701365/71‑1 (MQ=255)

                                                        ttCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttt               <  1:5655573/71‑1 (MQ=255)

                                                          ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                                       >  1:2721395/1‑45 (MQ=255)

                                                          ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                                       >  1:3134255/1‑45 (MQ=255)

                                                          ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt              >  1:1640732/1‑70 (MQ=255)

                                                          ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt              <  1:5499503/70‑1 (MQ=255)

                                                           tcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                <  1:2005458/51‑1 (MQ=255)

                                                           tcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                          <  1:2449680/57‑1 (MQ=255)

                                                            cATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttttt             >  1:3973005/1‑69 (MQ=255)

                                                            cATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAAATTCTTGAAGTTCAGCGGCAGttttt             >  1:3973945/1‑69 (MQ=255)

                                                             aTACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttttt             <  1:2829398/68‑1 (MQ=255)

                                                             aTACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttttt             <  1:4208587/68‑1 (MQ=255)

                                                                   cAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGt                            >  1:3659735/1‑47 (MQ=255)

                                                                   cAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACca      >  1:4997603/1‑69 (MQ=255)

                                                                    agagAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAg   <  1:2205120/71‑1 (MQ=255)

                                                                     gagaAGAGGTGTCGATGGGCACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGa  <  1:2248125/71‑1 (MQ=255)

                                                                     gagaAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCa                   >  1:439162/1‑54 (MQ=255)

                                                                      agaagaGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAg   <  1:1388042/69‑1 (MQ=255)

CGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGA  >  minE/608993‑609132

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