BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,913	G→A	21.5%	I447I (ATC→ATT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,913	0	G	A	21.5%	240.6 / 56.1	121	I447I (ATC→ATT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (50/45); new base A (14/12); total (64/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

GGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATC  >  minE/608843‑608982

ggACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTg                                                                       <  1:3891587/71‑1 (MQ=255)

  aCATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAt                                                                     >  1:1881430/1‑71 (MQ=255)

    aTCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGa                                                                      <  1:2356034/68‑1 (MQ=255)

     tCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGc                                                                  <  1:2574888/71‑1 (MQ=255)

     tCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGc                                                                  <  1:4092283/71‑1 (MQ=255)

      cGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                                 >  1:3047628/1‑71 (MQ=255)

      cGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCg                                                                 >  1:2707379/1‑71 (MQ=255)

        ccTTTGATCGGTTCAGACTTCTGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                                               <  1:2881538/71‑1 (MQ=255)

        ccTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGAt                                                                     >  1:4518497/1‑65 (MQ=255)

        ccTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                                               <  1:2868979/71‑1 (MQ=255)

        ccTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                                               <  1:4510550/71‑1 (MQ=255)

        ccTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                                               <  1:4843710/71‑1 (MQ=255)

          tttGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                                               >  1:1517222/1‑69 (MQ=255)

             gATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGa                                                                      >  1:1843127/1‑59 (MQ=255)

             gATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGa                                                                      >  1:2231155/1‑59 (MQ=255)

                 ggTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAAAGCCGCCGTTGATTGCGTACAGCATGGt                                                      <  1:3284944/71‑1 (MQ=255)

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

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

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

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

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

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

                          ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGa                                                                      <  1:4029474/46‑1 (MQ=255)

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

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

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

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

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

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

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

                          ttCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGtt                                                     >  1:2264402/1‑63 (MQ=255)

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

                            cGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATg                                                        >  1:3448639/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)

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

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

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

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

                                       gCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTa                                                               >  1:2643916/1‑40 (MQ=255)

                                       gCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                                >  1:2833316/1‑71 (MQ=255)

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

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

                                            ttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                             >  1:55847/1‑69 (MQ=255)

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

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

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

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

                                              cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                          <  1:5057062/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:899933/1‑41 (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:2433046/1‑46 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                      cGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                 >  1:462022/1‑71 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                 ccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtt                                          <  1:3790739/35‑1 (MQ=16)

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

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

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

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

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

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

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

GGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATC  >  minE/608843‑608982

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