Blastn is not working for short sequences?
1
0
Entering edit mode
5.0 years ago
Kumar ▴ 120

I am trying to do blastn of two fasta files,which includes short sequence containing query fasta file (prom.fasta) and a bacterial draft genome (bac.fasta) as subject as follows,

prom.fasta
>og1_1
TTCGGTTCTAACGGCTGGGATTCGG

and

bac.fasta
>test
ATCAGCGTGGATTTGCCGGCGCCGGAATGGCCGATGATGCCGAATACTTCGCCGGGGCCG
ATCCGCAGATCGAGCGGATGCAGCGCCACGATCTGGCGACCGTCGACGGAATAGGACTTG
TGCAGGCGCTGAAACTGGATCACCGCGCGTGGCCGATTGGGGGGTCGTGAAGCCTACCAG
TCAATGTTTGGGGCCGTCATTCCATTCGGTTCTAACGGCTGGGATTCGGGAGTCGGGATT
GGGGATTCGCAAGAGCGGACGCCCTGGCCTTTGCGAATCCCGAGCCTCCAATCCCTAACG

I have used the following command for blastn analysis

blastn -query prom.fasta -subject bac.fasta -outfmt 6 -out result.xls

blastn gives no output though, the bac.fasta file having the exact similar sequence of prom.fasta. But, I do not know, why it is not giving output. could anyone help me to fix it.

software-error fasta blast sequence • 1.3k views
ADD COMMENT
2
Entering edit mode
5.0 years ago
GenoMax 147k

You will have to use an additional parameter -task blastn-short for short sequences.

blastn-short    BLASTN program optimized for sequences shorter than 50 bases
ADD COMMENT

Login before adding your answer.

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