Finding the right database for blast CLI
0
0
Entering edit mode
10 months ago
cf556 ▴ 10

We are looking into the possibility of off-target mapping for some sequences that were added to our alignment reference. We are only concerned about off-target mapping to genes (including introns). We are also only concerned about off-target mapping in mouse.

To look at this I want to do a blast search on ~1000 sequences that mapped to our added sequence, only to the mouse genome and only to genes (including introns). I have a fasta file for mm10 and a gtf file describing the genes, but I am unsure how to put this into blast + (blast CLI) so that I know that I am searching genes.

Thanks!

blast • 197 views
ADD COMMENT
0
Entering edit mode

You can make your blast database using makeblastdb and align your sequences against that database. Use -outfmt 6 for the output. Convert tabular blast+ output to BED format (Convert blast tabular output to bed format ). Then you can intersect BED output with a BED file for gene models.

A command line manual for blast+ is available here.

ADD REPLY

Login before adding your answer.

Traffic: 2047 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6