blastread with local blast output v2.2.31+
Afficher commentaires plus anciens
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
le 15 Jan 2016
Réponses (0)
Catégories
En savoir plus sur Cell Arrays dans Centre d'aide et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!