SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
How to fix fastq files from a bacteria sample? ymc Bioinformatics 2 05-30-2015 08:47 AM
How to fix Illumina FASTQ files with read length varies errors antgomo Illumina/Solexa 5 08-19-2014 11:28 PM
cummeRbund: combine data sets / multiple experiments on one graph? orthodoc RNA Sequencing 1 02-14-2014 05:10 PM
Making a Blast Script for FASTA files utagenomics Bioinformatics 8 04-26-2013 11:29 AM
Making bam files in Pysam Martino Bioinformatics 1 06-28-2012 12:09 AM

Reply
 
Thread Tools
Old 10-15-2015, 11:05 AM   #1
rbruenn
Junior Member
 
Location: Bay Area, California

Join Date: Oct 2015
Posts: 5
Default makeblastdb is making multiple sets of files, how can I fix it?

Hello!

I just realized I should have posted this to this forum, so here it is again, sorry!

I am a graduate student at UC Berkeley currently working with raw reads of several transcriptomes in an attempt to find and assemble reads that match a couple of genes I'm studying. This site has already been very useful to me (Thank you!!!!), but I haven't found any answers pertaining to my current problem.

I was hoping I could get some help with a BLAST problem I'm having. I am working with standalone blast, and am building blast databases from fasta files of transcriptome raw reads. I have successfully used the command:

makeblastdb -in COST_1_final.fasta -input_type fasta -dbtype nucl -out COST_1_final for several fasta files of raw reads (with varying names of course) but a few of the files result in multiple sets of database files, marked for example <filename>.00.nhr and <filename>.01.nhr with what I believe is an alias file <filename>.nal

The command as it runs give the usual message:
Building a new DB, current time: 10/15/2015 10:49:55
New DB name: COST_u_final
New DB title: COST_u_final.fasta
Sequence type: Nucleotide
Keep Linkouts: T
Keep MBits: T
Maximum file size: 1000000000B
Adding sequences from FASTA; added 16313457 sequences in 937.429 seconds.

but then results in multiple sets of files.

Any ideas about how I can make just one set of database files, as has successfully happened with the rest of my fasta files? I would really appreciate any help!!
rbruenn is offline   Reply With Quote
Old 10-15-2015, 11:14 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,022
Default

Solved here: http://seqanswers.com/forums/showthread.php?t=63477

Last edited by GenoMax; 10-15-2015 at 11:35 AM.
GenoMax is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 11:45 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO