BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,200,118	C→T	intergenic (+2820/‑161)	yehF → / → yehI	putative protein YehF/DUF4132 domain‑containing protein YehI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,200,118	0	C	T	100.0%	89.8 / NA	26	intergenic (+2820/‑161)	yehF/yehI	putative protein YehF/DUF4132 domain‑containing protein YehI
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (10/16); total (10/16)

GGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGG  >  NC_000913/2200047‑2200195

gggTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGct                                                                      >  1:373841/1‑81 (MQ=255)

  gTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGctct                                                                    <  2:373802/81‑1 (MQ=255)

     gCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTa                                                                  <  2:373804/80‑1 (MQ=255)

     gCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTa                                                                  <  1:373803/80‑1 (MQ=255)

                aTCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCtgtg                                                       <  2:373805/80‑1 (MQ=255)

                 tCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGc                                                     <  1:373806/81‑1 (MQ=255)

                 tCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGc                                                     <  2:373807/81‑1 (MQ=255)

                   aTGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCt                                                    <  1:373808/80‑1 (MQ=255)

                   aTGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCt                                                    <  2:373809/80‑1 (MQ=255)

                     ggAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATg                                                 >  1:373838/1‑81 (MQ=255)

                      gAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATg                                                 <  1:373810/80‑1 (MQ=255)

                            gATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAgg                                           >  1:373862/1‑80 (MQ=255)

                             aTGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGc                                          >  2:373842/1‑80 (MQ=255)

                                         cGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGat                             <  1:373811/81‑1 (MQ=255)

                                          ggTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGata                            >  1:373850/1‑81 (MQ=255)

                                                 aTGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTgg                     <  2:373812/81‑1 (MQ=255)

                                                    tCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTggggg                  <  2:373813/81‑1 (MQ=255)

                                                    tCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTgggg                   <  2:373814/80‑1 (MQ=255)

                                                     cGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTgggg                   >  1:373844/1‑79 (MQ=255)

                                                      gCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGa                 <  1:373815/80‑1 (MQ=255)

                                                         aGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTaa             <  2:373816/81‑1 (MQ=255)

                                                          gCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAAt            >  1:373849/1‑81 (MQ=255)

                                                          gCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAAt            >  2:373847/1‑81 (MQ=255)

                                                             gAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTc          >  2:373851/1‑80 (MQ=255)

                                                               aCTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAac       <  1:373817/81‑1 (MQ=255)

                                                                   aGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAg   >  2:373857/1‑81 (MQ=255)

                                                                    gcgtTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAgg  <  1:373818/81‑1 (MQ=255)

GGGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGG  >  NC_000913/2200047‑2200195

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