BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,922	C→G	26.4%	L444L (CTG→CTC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,922	0	C	G	26.4%	210.4 / 79.1	121	L444L (CTG→CTC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (46/43); new base G (18/14); total (64/57)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 127 total aligned reads displayed.

GGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTA  >  Exported/608860‑608992

 gTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtt                                                               <  1:65963/71‑1 (MQ=255)

  ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                                              <  1:5431602/71‑1 (MQ=255)

  ttCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGttt                                                              <  1:3358694/71‑1 (MQ=255)

    cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                                             <  1:2553028/70‑1 (MQ=255)

    cAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTc                                                            >  1:4130640/1‑71 (MQ=255)

      gACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGc                                                          >  1:1343391/1‑71 (MQ=255)

       aCTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                                                         <  1:1821597/71‑1 (MQ=255)

         ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCa                                                                    <  1:3575453/58‑1 (MQ=255)

         ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCa                                                                    <  1:2527513/58‑1 (MQ=255)

         ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCa                                                                    <  1:1190732/58‑1 (MQ=255)

         ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCa                                                                    <  1:4165177/58‑1 (MQ=255)

         ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                                                  <  1:1342268/60‑1 (MQ=255)

         ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                                                  <  1:4931906/60‑1 (MQ=255)

           cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                                                  >  1:4370299/1‑58 (MQ=255)

           cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGc                                                          >  1:2451882/1‑66 (MQ=255)

           cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                                                      <  1:3010404/70‑1 (MQ=255)

           cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAgg                                                      <  1:1832260/70‑1 (MQ=255)

           cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAg                                                       <  1:808652/69‑1 (MQ=255)

           cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                                     >  1:3961435/1‑71 (MQ=255)

           cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGt                                                     >  1:4920491/1‑71 (MQ=255)

                cgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcg                                                >  1:5607572/3‑70 (MQ=255)

                cAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                                >  1:2086234/1‑71 (MQ=255)

                 gaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgc                                               >  1:4024147/2‑69 (MQ=255)

                  aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGcc                                                         >  1:2503781/1‑60 (MQ=255)

                  aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                              <  1:2942887/71‑1 (MQ=255)

                  aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                              <  1:4610772/71‑1 (MQ=255)

                  aCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                                              <  1:62444/71‑1 (MQ=255)

                    cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                                            >  1:2418760/1‑71 (MQ=255)

                    cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                                            <  1:3899536/71‑1 (MQ=255)

                    cTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcgtg                                            <  1:604767/71‑3 (MQ=255)

                       cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAc                                         >  1:664257/1‑71 (MQ=255)

                         ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAcc                                        >  1:3300909/1‑70 (MQ=255)

                           ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAg                                                                      <  1:3972192/38‑1 (MQ=255)

                            tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                       >  1:4919861/1‑68 (MQ=255)

                             cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcgtg                                            >  1:5664726/1‑60 (MQ=255)

                             cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                    <  1:5057062/70‑1 (MQ=255)

                             cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                    <  1:3515728/70‑1 (MQ=255)

                             cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                    <  1:5081114/70‑1 (MQ=255)

                             cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                    <  1:3217052/70‑1 (MQ=255)

                              aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                  >  1:366141/1‑71 (MQ=255)

                              aGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                  <  1:2485728/71‑1 (MQ=255)

                                tttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACt                                <  1:3820837/71‑1 (MQ=255)

                                 ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtttt                                                             >  1:1781154/1‑41 (MQ=255)

                                 ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                               >  1:4708645/1‑71 (MQ=255)

                                 ttctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTg                               >  1:5586913/4‑71 (MQ=255)

                                  tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAg                                                       >  1:2557281/1‑46 (MQ=255)

                                  tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAg                                                       >  1:42513/1‑46 (MQ=255)

                                  tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGtt                                                    >  1:3815029/1‑49 (MQ=255)

                                  tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                               >  1:2813492/1‑70 (MQ=255)

                                  tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                               >  1:3031028/1‑70 (MQ=255)

                                  tttCGTCAACGCCACCGTTGATTGCGTACAGCATGGTTTTCGCCAGGtt                                                    >  1:1673985/1‑49 (MQ=255)

                                   ctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                              <  1:5301396/69‑1 (MQ=255)

                                    tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                            >  1:2057723/1‑71 (MQ=255)

                                    tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                            >  1:1163239/1‑71 (MQ=255)

                                    tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                  >  1:4648764/1‑65 (MQ=255)

                                    tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                  >  1:4202879/1‑65 (MQ=255)

                                     cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                            <  1:3273245/70‑1 (MQ=255)

                                       tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTg                         <  1:5268998/71‑1 (MQ=255)

                                       tCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTg                         <  1:1694310/71‑1 (MQ=255)

                                         aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                       >  1:2669789/1‑71 (MQ=255)

                                         aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                       >  1:2402722/1‑71 (MQ=255)

                                         aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                       <  1:898099/71‑1 (MQ=255)

                                          aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACTCGCACCGAAGAACTGCATTTGtt                       <  1:2118390/70‑1 (MQ=255)

                                          aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                       >  1:3459613/1‑70 (MQ=255)

                                           ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                             <  1:1155144/61‑1 (MQ=255)

                                           ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                             <  1:2393552/61‑1 (MQ=255)

                                            cccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAc                    >  1:1557742/2‑71 (MQ=255)

                                            cccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAc                    >  1:3756992/2‑71 (MQ=255)

                                             ccgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTAc                    >  1:5557805/1‑70 (MQ=255)

                                              cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                    <  1:4665827/53‑1 (MQ=255)

                                              cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTgct                       >  1:1081659/1‑64 (MQ=255)

                                              cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTgct                       >  1:3426410/1‑64 (MQ=255)

                                              cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                   >  1:4221099/1‑70 (MQ=255)

                                              cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                   >  1:5186119/1‑70 (MQ=255)

                                                ccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                     >  1:5516307/1‑66 (MQ=255)

                                                ccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaat                <  1:2375314/71‑2 (MQ=255)

                                                 cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg               <  1:775996/71‑1 (MQ=255)

                                                 cgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACg               <  1:2771617/71‑1 (MQ=255)

                                                   ttGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGat             >  1:4718237/1‑71 (MQ=255)

                                                    taaTTGCTTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                             <  1:438023/52‑1 (MQ=255)

                                                    tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaa                                     >  1:3720903/1‑46 (MQ=255)

                                                     gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc            <  1:4320014/70‑1 (MQ=255)

                                                     gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc            <  1:2818379/70‑1 (MQ=255)

                                                     gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc            <  1:4810051/70‑1 (MQ=255)

                                                      aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca           <  1:3114830/70‑1 (MQ=255)

                                                      aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca           <  1:261160/70‑1 (MQ=255)

                                                        tGCGTACAGCATGGTTTTCCCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTa                     >  1:3366110/1‑58 (MQ=255)

                                                         gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                                       >  1:4149446/1‑38 (MQ=9)

                                                         gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                                       >  1:1466690/1‑38 (MQ=9)

                                                         gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaagaa                                  >  1:4326718/1‑44 (MQ=12)

                                                         gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg       >  1:5340523/1‑71 (MQ=255)

                                                         gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg       >  1:3865652/1‑71 (MQ=255)

                                                          catACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                  >  1:2913451/3‑43 (MQ=255)

                                                          cGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                  >  1:2067710/1‑43 (MQ=255)

                                                            tACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt    <  1:4018357/71‑1 (MQ=255)

                                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                 >  1:3498844/1‑40 (MQ=255)

                                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                 >  1:136434/1‑40 (MQ=255)

                                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa  >  1:5663115/1‑71 (MQ=255)

                                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa  >  1:2240311/1‑71 (MQ=255)

                                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa  >  1:5572364/1‑71 (MQ=255)

GGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTA  >  Exported/608860‑608992

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