This is a test version of Biostars. For the public version, visit https://www.biostars.org.
BLAST masks repetitive sequences - but where: database or query sequence or both?

I understand that BLAST+ masks repetitive sequences. However, I cannot find out if the repetitive sequences are masked in the database or the query sequence or both? [reference would be most appreciated].

Secondly, would BLAST mask a sequence like the following:

QFSAGKRQFSAGKRQFSAGKRQFSAGKRQFSAGKRQWLGGEEEYDPEENLNMETRQFSAGKRQFSAGKRQFSAGKRDWEEELTPEELMDMFQAPETRQFSAGKRQFSAGKRQFSAGKRQWVGGEEEYDPEEMLNMATRQFSAGKRQFSAGKRQFSAGKRQWVGGEEAFLPEMDTRQFSAGKRQFSAGKRQFSAGKRQFSAGKR

where QFSAGKR is repeated 18 times...

BTW the above sequence is thyrotropin releasing hormone where the functional bioactive neuropeptide can be as small as only three amino acids long which is repeated multiple times.

Also, would be nice, if someone could clarify what masks means - i.e. does this mean BLAST totally ignore this sequence (i.e. replace with N/X) or simply give such sequences a lower e-value...

blast
• 5,660 views • link
updated by Ram Ram • written by IsmailM IsmailM

In fact you should be able to find this information on the standalone blast manual page. Masking could be soft/lowercase so preserving sequence or hard - replacing nucleotides with Ns or aminoacids with Xs. Algorithm is presented friendly in "BLAST" by Korf, Yandell and Bedell, however it is quite old book.

• 0 views • link
written by boczniak767 boczniak767

1 answer

Blast may mask both query and database:

-dust for masking query

and -db_soft_mask for masking the database

It is very easy for you to check if blast would mask thyrotropin, just go to Blast page and play with it. If you go to the bottom of the page, clicking on "+Algorithm parameters" will expand several option, including masking ones.

• 0 views • link
written by h.mon h.mon

Log in to answer this question.

More posts like this