BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,139,485	A→G	K253K (AAA→AAG)	yghU →	putative S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,139,485	0	A	G	100.0%	125.1 / NA	34	K253K (AAA→AAG)	yghU	putative S‑transferase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (18/16); total (18/16)

GGATATGGCGATTTGGCCGTGGTTTGGCAACGTGGTGTTAGGTGGTGTGTATGATGCCGCTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATA  >  NC_000913/3139350‑3139632

ggATATGGCGATTTGGCCGTGGTTTGGCAACGTGGTGTTAGGTGGTGTGTATGATGCCGCTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGtt                                                                                                                                      <  1:409484/151‑1 (MQ=255)

            ttGGCCGTGGTTTGGCAACGTGGTGTTAGGTGGTGTGTATGATGCCGCTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAAc                                                                                                                          >  1:462631/1‑151 (MQ=255)

              ggCCGTGGTTTGGCAACGTGGTGTTAGGTGGTGTGTATGATGCCGCTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCTCCAACg                                                                                                                         <  2:363515/150‑1 (MQ=255)

                    ggTTTGGCAACGTGGTGTTAGGTGGTGTGTATGATGCCGCTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGc                                                                                                                   <  2:231316/150‑1 (MQ=255)

                    ggTTTGGCAACGTGGTGTTAGGTGGTGTGTATGATGCCGCTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCt                                                                                                                  >  1:652652/1‑151 (MQ=255)

                             aCGTGGTGTTAGGTGGTGTGTATGATGCCGCTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCa                                                                                                         <  1:349630/151‑1 (MQ=255)

                               gtggtgTTAGGTGGTGTGGATGATGCCGCTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGAGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCa                                                                                                         >  2:568721/1‑149 (MQ=255)

                                       aGGTGGTGTGTATGATGCCGCTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGa                                                                                                <  2:44105/150‑1 (MQ=255)

                                           ggtgTGTATGATGCCGCTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGcc                                                                                           >  2:184407/1‑151 (MQ=255)

                                             tgtgtATGATGCCGCTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCAt                                                                                         >  2:269988/1‑151 (MQ=255)

                                                     atgCCGCTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGcc                                                                                   >  2:129432/1‑149 (MQ=255)

                                                           cTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGAttt                                                                           >  1:7876/1‑151 (MQ=255)

                                                              aGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCg                                                                         >  2:605725/1‑150 (MQ=255)

                                                              aGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCg                                                                         >  1:433383/1‑150 (MQ=255)

                                                               gTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCg                                                                         >  1:671483/1‑149 (MQ=255)

                                                                     tGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAAt                                                                 >  1:293947/1‑151 (MQ=255)

                                                                      gATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATa                                                                <  1:253818/151‑1 (MQ=255)

                                                                       aTGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATAc                                                               >  2:334051/1‑151 (MQ=255)

                                                                         gCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACgg                                                             <  1:539947/151‑1 (MQ=255)

                                                                         gCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACg                                                              <  2:590719/150‑1 (MQ=255)

                                                                                   ataAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCg                                                   <  1:208767/151‑1 (MQ=255)

                                                                                     aaGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTc                                                 >  1:436261/1‑151 (MQ=255)

                                                                                             aCAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGAGGCCAGTGATTTCTAGACGAATACGGAAGATAAGCGTCAGGGGt                                           <  2:122147/149‑1 (MQ=255)

                                                                                             aCAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTa                                          >  1:146515/1‑150 (MQ=255)

                                                                                                  gCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGtt                                    <  2:420822/151‑1 (MQ=255)

                                                                                                     gggCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTg                                   <  2:24336/149‑1 (MQ=255)

                                                                                                         gAAAGAAGTAGGCGAACGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTtc                             >  2:343879/1‑151 (MQ=255)

                                                                                                                        aCGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCAt              <  1:664805/151‑1 (MQ=255)

                                                                                                                        aCGTCCGGCGGTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCAt              >  1:139140/1‑151 (MQ=255)

                                                                                                                               gcggTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGcc       >  1:29322/1‑151 (MQ=255)

                                                                                                                               gcggTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGcc       <  2:525338/151‑1 (MQ=255)

                                                                                                                               gcggTGAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGcc       >  1:524978/1‑151 (MQ=255)

                                                                                                                                    gAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGAt   <  1:592856/150‑1 (MQ=255)

                                                                                                                                    gAAGCGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATa  <  2:175309/151‑1 (MQ=255)

GGATATGGCGATTTGGCCGTGGTTTGGCAACGTGGTGTTAGGTGGTGTGTATGATGCCGCTGAGTTTCTTGATGCGGGCAGTTATAAGCATGTACAACGCTGGGCGAAAGAAGTAGGCGAACGTCCGGCGGTGAAACGTGGGCGTATTGTTAACCGCACCAACGGACCGCTGAATGAGCAGTTGCATGAGCGCCATGACGCCAGTGATTTCGAGACGAATACGGAAGATAAGCGTCAGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATA  >  NC_000913/3139350‑3139632

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