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?
blastread with local blast output v2.2.31+
5 vues (au cours des 30 derniers jours)
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
Réponses (0)
Voir également
Catégories
En savoir plus sur Cell Arrays dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!