Quantcast
Channel: Post Feed
Viewing all articles
Browse latest Browse all 41826

Problem using a custom blast database in tblastx

$
0
0
Hello. I tried to create a custom DB for blast. I downloaded fasta files from NCBI from an SRA project and then run this command  makeblastdb -dbtype nucl -in sra_data-DB.fasta -input_type fasta -out plant_Transcriptome_DB -title plant_Transcriptome_DB and all seemed that went okey. Building a new DB, current time: 06/18/2014 17:22:49 New DB name:   plant_Transcriptome_DB New DB title:  plant_transcriptome_DB Sequence type: Nucleotide Keep Linkouts: T Keep MBits: T Maximum file size: 1000000000B Adding sequences from FASTA; added 411965 sequences in 55.2244 seconds. Three files created. (nhr,nin,nsq) After this, i run a test with blastdbcmd to see if everything went right. blastdbcmd -db plant_Transcriptome_DB.nsq -info  and returned me this: Database: plant_transcriptome_DB 411,965 sequences; 107,369,081 total bases Date: Jun 18, 2014  5:22 PM    Longest sequence: 999 bases Volumes:   /*/*/*/*/*/*/*/plant_Transcriptome_DB   Ok! lets run tblastx and see what is happening: tblastx -query contig2-query.fasta -db plant_Transcriptome_DB -out output -html  Return : Error: (1431.1) FASTA-Reader: Warning: FASTA-Reader: Ignoring invalid residues at position(s): On line 2: 6, 8, 13, 33, 43-44, 47-48, 50, 55-56, 58-59, 68-70, 72, 81 Error: (1431.1) FASTA-Reader: Warning: FASTA-Reader: Ignoring invalid residues at position(s): On line 3: 7, 16, 18, 20-22, 26-28, 30, 35, 39, 41-43, 45, 47, 49, 54, 56-57, 59-60, 63-64, 66-67, 69, 71, ...

Viewing all articles
Browse latest Browse all 41826

Trending Articles