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%	74.0 / NA	21	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 (11/10); total (11/10)

GGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACAC  >  NC_000913/2200048‑2200192

ggTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGctc                                                                  <  1:399755/81‑1 (MQ=255)

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

    gCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCtt                                                                >  2:399793/1‑79 (MQ=255)

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

    gCTGGTCGTTAATCATGGAAATCGATGAAGGATTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTaa                                                              <  1:399757/81‑1 (MQ=255)

       ggTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTaagc                                                            >  2:399773/1‑80 (MQ=255)

           gTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCt                                                       >  1:399795/1‑81 (MQ=255)

                      aaaTCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTCGAGCAGCTCTTAAGCAAGCTGTGGCTATgg                                             >  2:399782/1‑80 (MQ=255)

                       aaTCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATggg                                            >  2:399806/1‑80 (MQ=255)

                          cGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAgg                                        >  2:399792/1‑81 (MQ=255)

                             tGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCaa                                     <  1:399759/81‑1 (MQ=255)

                               aaGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAg                                   <  1:399760/81‑1 (MQ=255)

                               aaGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAg                                   >  2:399798/1‑81 (MQ=255)

                                 ggCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAg                                   >  2:399789/1‑79 (MQ=255)

                                 ggCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCa                                 <  1:399761/81‑1 (MQ=255)

                                           tAGGCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATg                       <  1:399762/81‑1 (MQ=255)

                                              gCATGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCt                    <  1:399763/81‑1 (MQ=255)

                                                aTGTCGCCAGCCGAACTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTgg                  >  2:399802/1‑81 (MQ=255)

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

                                                               cTCAGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAaca   >  1:399803/1‑81 (MQ=255)

                                                                  aGCGTTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAacac  <  1:399764/79‑1 (MQ=255)

GGTCGCTGGTCGTTAATCATGGAAATCGATGAAGGCTTTGCGGTAGGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACAC  >  NC_000913/2200048‑2200192

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