blastread with local blast output v2.2.31+

Hello, I'm trying to read in the results of an alignment that I performed using a local blast program (v2.2.31+) and against the supplied human genome from NCBI. I used the a discontinuous megablast ('dc-megablast') algorithm and obtained an output text file. I didn't modify the output at all but when I try to read it into matlab using the command:
data = blastread('SRF_mouse_human.txt')
I get the error:
Attempted to access stats_start(0); index must be a positive integer or logical.
Error in blastread (line 125) statstext = blasttext(stats_start(size(stats_start,1)):size(blasttext,2));
Is the problem that I'm using dc-megablast as opposed to blastn? Am I not asking my local instance of blast for the right output format? Wrong version of BLAST? Any thoughts very much appreciated.
Dwight

1 commentaire

Dwight Chambers
Dwight Chambers le 15 Jan 2016
As a followup, I'm able to successfully read in single query entries into MATLAB either from a locally sourced version of BLAST or NCBI's website. Is it that I cannot read multiple BLAST results simultaneously?

Connectez-vous pour commenter.

Réponses (0)

Catégories

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by