This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Removal of gene sequences from genome.

I designed a bait set to caputure genes but i would like to know whether my baits bind on other places other than the designed gene sequences. My idea is to blast it back to the reference genome and identify the "background noise" of the non gene sequences. An easy solution i think is to remove all gene sequences from my reference genome. Leaving only the intergenic and repetitive regions.

Does someone have a easy solution for this? Or other suggestion i could try?

blast genome gff
• 1,695 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by Candah Candah

1 answer

get a BED for your genes and then use bedtools mask fasta http://bedtools.readthedocs.io/en/latest/content/tools/maskfasta.html to replace the bases with N.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

I like it, thanks for your suggestion

• 0 views • link
written by Candah Candah

Log in to answer this question.

More posts like this