blastn with short sequences gives no results in the corresponding genome
1
2
Entering edit mode
11 weeks ago
Audrey ▴ 50

Hi all,

I've designed primers to identify a strain of interest.

Just to be sure that the primers matche only the genome of this strain and no other genomes, I am using blastn analysis.

My input file is :

>Seq_1 
ATGTGCTAGCTATGGCG
>Seq_2 
GTCGGTAATCGATTCGA

And my database is a multifasta file containing several genomes including the genome of my strain of interest.

My command line is the following :

> blastn -query Sequences.fasta -db genomes.fasta -out Results_blastn.txt -outfmt 6 -evalue 0.1

The results file is blank.

I do not understand how I've got no results even on the genome of my strain of interest ...

Do someone have an idea on how I should proceed to fix this issue?

Thank you so much in advance for your much appreciated help!

Audrey

blastn output primers • 162 views
ADD COMMENT
0
Entering edit mode

You might want to take a look specifically at PrimerBlast (https://www.ncbi.nlm.nih.gov/tools/primer-blast/) which is designed specifically for this.

ADD REPLY
4
Entering edit mode
11 weeks ago

Just add -task blastn-short to your command. By default, blastn uses megablast (optimized for highly similar long sequences).

ADD COMMENT
0
Entering edit mode

Indeed, it works better now ... Thank you so much for your help a.zielezinski!

ADD REPLY

Login before adding your answer.

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