BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,117,047	G→T	100%	S59Y (TCC→TAC)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,047	0	G	T	100.0%	139.1 / NA	39	S59Y (TCC→TAC)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base T (20/19); total (20/19)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGC  >  NC_000913/4116905‑4117185

gCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCttt                                                                                                                                      >  2:187050/1‑149 (MQ=255)

 cACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTg                                                                                                                                     >  2:285118/1‑149 (MQ=255)

        aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGcc                                                                                                                              >  1:193966/1‑149 (MQ=255)

        aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGc                                                                                                                               <  2:255244/148‑1 (MQ=255)

        aGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGc                                                                                                                               >  1:255244/1‑148 (MQ=255)

            aaTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAg                                                                                                                           <  4:56687/148‑1 (MQ=255)

            aaTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAg                                                                                                                           <  4:283039/148‑1 (MQ=255)

             aTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGa                                                                                                                          <  3:298804/148‑1 (MQ=255)

             aTGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGAt                                                                                                                         <  1:332312/149‑1 (MQ=255)

               ggCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGAttt                                                                                                                       <  1:285118/149‑1 (MQ=255)

                                       ttCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCaa                                                                                                >  3:85975/1‑148 (MQ=255)

                                          tttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTg                                                                                            >  1:163034/1‑149 (MQ=255)

                                           ttttCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTgg                                                                                           <  1:233005/149‑1 (MQ=255)

                                                 cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                                                       <  4:26990/147‑1 (MQ=255)

                                                 cAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGtttt                                                                                       >  3:26990/1‑147 (MQ=255)

                                                  aGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTgg                                                                                                                                    >  4:7626/1‑101 (MQ=255)

                                                  aGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTgg                                                                                                                                    <  3:7626/101‑1 (MQ=255)

                                                          gTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGAttt                                                                             <  4:207646/148‑1 (MQ=255)

                                                                tCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCaaa                                                                      <  1:6477/149‑1 (MQ=255)

                                                                   cGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAAtt                                                                    >  3:19383/1‑148 (MQ=255)

                                                                         ttCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGt                                                                                                              >  2:315022/1‑100 (MQ=255)

                                                                         ttCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGt                                                                                                              <  1:315022/100‑1 (MQ=255)

                                                                          tCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTGCCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTg                                                             >  3:285229/1‑148 (MQ=255)

                                                                             tAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGa                                                          <  4:159286/148‑1 (MQ=255)

                                                                               aTACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGaca                                                        >  3:322803/1‑148 (MQ=255)

                                                                                       aGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTGCCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATg                                                <  4:285229/148‑1 (MQ=255)

                                                                                           gCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAg                                                                                 <  3:188855/111‑1 (MQ=255)

                                                                                           gCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAg                                                                                 >  4:188855/1‑111 (MQ=255)

                                                                                                      aaaaaatatatCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGAt                                >  2:241657/8‑149 (MQ=255)

                                                                                                       gAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATc                               <  2:231071/149‑1 (MQ=255)

                                                                                                          cTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCa                             <  4:19383/148‑1 (MQ=255)

                                                                                                           tGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCAtt                           >  2:269355/1‑149 (MQ=255)

                                                                                                                tCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATcc                              <  1:338903/141‑1 (MQ=255)

                                                                                                                tCGGCTTTCGCCAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATcc                              >  2:338903/1‑141 (MQ=255)

                                                                                                                          ccAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg            >  2:381245/1‑149 (MQ=255)

                                                                                                                           cAGCACTTCTACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGc           >  1:143014/1‑149 (MQ=255)

                                                                                                                                    tACCAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGc  >  1:141876/1‑149 (MQ=255)

                                                                                                                                      ccAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAAt                                                 >  2:329439/1‑100 (MQ=255)

                                                                                                                                      ccAGCGTGTAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAAt                                                 <  1:329439/100‑1 (MQ=255)

GCACTGCCAGACAATGGCGTTATAGATAGGCTTGCCGGTTTCTTTTTCCCAGACAATAGTGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGC  >  NC_000913/4116905‑4117185

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`