blastread with local blast output v2.2.31+

5 visualizaciones (últimos 30 días)
Dwight Chambers
Dwight Chambers el 14 de En. de 2016
Comentada: Dwight Chambers el 15 de En. de 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 comentario
Dwight Chambers
Dwight Chambers el 15 de En. de 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?

Iniciar sesión para comentar.

Respuestas (0)

Categorías

Más información sobre Cell Arrays en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

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

Start Hunting!

Translated by