blastread with local blast output v2.2.31+

5 vues (au cours des 30 derniers jours)
Dwight Chambers
Dwight Chambers le 14 Jan 2016
Commenté : Dwight Chambers le 15 Jan 2016
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

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!

Translated by