Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Makeblast db with more sequences in fasta format
1
gravatar for marghe08
4.2 years ago by
marghe08 • 10
marghe08 • 10 wrote:

Hello! I am a beginner and I need help. I downloaded from NCBI two distinct fasta files, which contain both a nucleotide sequence, and are saved in a folder on my pc; these are the names of the files: NG_006672.1.fasta NG_006682.3.fasta using locally blast would like to create a database, with makeblastdb command, which contains my 2 fasta sequences. How can I do?
blast blastdb • 1.8k views
ADD COMMENT • link •
modified 15 months ago by david.prihoda • 10 • written 4.2 years ago by marghe08 • 10
4
gravatar for Ram
4.2 years ago by
Ram ♦ 32k
Baylor College of Medicine, Houston, TX
Ram ♦ 32k wrote:

Just cat the fasta file and use the concatenated FASTA file as the input to the -in argument. Here's NCBI manual on how to use makeblastdb https://www.ncbi.nlm.nih.gov/books/NBK279688/

Either cat them as a preliminary step or use the syntax -in <(cat file1.fa file2.fa) to use a subshell in the makeblastdb
ADD COMMENT • link modified 15 months ago • written 4.2 years ago by Ram ♦ 32k
1

-in <(cat file1.fa file2.fa) seems not to work with Blast 2.7.1, it shows this error: BLAST options error: File /dev/fd/11 does not exist . I found out you can actually use a regular pipe like so: cat file1.fa file2.fa | makeblastdb -title myDBtitle -out myDBpath -parse_seqids

ADD REPLY • link written 15 months ago by david.prihoda • 10

That could be a shell restriction. If your shell allows process substitution, <() should work fine.

Also, this is a 3 year old post, so there is quite the chance the blast executable could be a different version.

ADD REPLY • link written 15 months ago by Ram ♦ 32k
1
gravatar for Sej Modha
4.2 years ago by
Sej Modha ♦ 4.8k
Glasgow, UK
Sej Modha ♦ 4.8k wrote:

https://www.ncbi.nlm.nih.gov/books/NBK279688/
ADD COMMENT • link written 4.2 years ago by Sej Modha ♦ 4.8k
0
gravatar for david.prihoda
15 months ago by
david.prihoda • 10
david.prihoda • 10 wrote:

I found out you can just omit the -in parameter, add a -title parameter and pipe the sequences using: cat file1.fa file2.fa | makeblastdb -title myDBtitle -out myDBpath -parse_seqids
ADD COMMENT • link written 15 months ago by david.prihoda • 10
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1001 users visited in the last hour
_