Blastn for gRNA characterization in ANopheles
1
0
Entering edit mode
3 months ago
Sbrillo • 0

Hi!

I'm trying to blast a dataset of gRNA against all the sequences present on the NCBI coming from Anopheles gambiae species. This is the script that I used:

blastn -db nt -query candidategRNA.fasta -task blastn-short -remote -evalue 10 -entrez_query 'Anopheles gambiae [organism]' -outfmt 6 -out G3-ill-G3-Nano-results.table 

Blast is creating the output file but it's empty... I already tried to blast my gRNAs manually and I always found matches with the web-browser Blastn version....

There is a way to use blast against all the sequences belonging to a specific organism and not just against the reference genome assembly (in my case Agamp3) ?

Or do you know how I can create my local database importing all the sequences belonging to Anopheles in the NCBI archive?

Any suggestions?

Thanks in advance!

Anopheles NCBI gRNA blastn • 275 views
ADD COMMENT
2
Entering edit mode
3 months ago
GenoMax 108k

Instead of entrez_query can you try limiting search using -taxids 7165 (which is the taxID of Anopheles gambiae). If that does not work you can try -taxids 7164 (which is Anopheles genus).

ADD COMMENT
0
Entering edit mode

Thanks a lot!

I have the output file with all the info that I needed, I will try to expand my blast result to Anopheles genus as you suggested :)

I have another question...

My input file contains more than 200 possible target sites and the output is quite big.

Do you have any suggestions on how to create a nice Plot of the results using R?

Thanks in advance

ADD REPLY
0
Entering edit mode

Don't what kind of plot you are looking for since some of it may depend on output format of your blast results. Came across this via a search see if it is useful: https://environmentalmicrobiome.biomedcentral.com/articles/10.1186/s40793-020-00361-y

ADD REPLY

Login before adding your answer.

Traffic: 1826 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