Consider the following (in branch blast_to_gtf of pygtftk). This branch contains the latest ligtftk release (develop branch).
from pygtftk.gtf_interface import GTF
def native_str(x):
return bytes(x.encode())
blast = GTF._dll.load_blast(native_str('blastn.blast'))
blast.size
[1] 22657 segmentation fault ipython
NB: in this branch, blastn.blast can be found in pygtftk/data/mini_real_10M/ (as a gzip file). This file was obtained by blasting the computed sequence of transcripts (mini_real_10M_tx_seq.fa) against the genome region (chr1_hg38_10M.fa).
Cheers
Consider the following (in branch blast_to_gtf of pygtftk). This branch contains the latest ligtftk release (develop branch).
NB: in this branch, blastn.blast can be found in pygtftk/data/mini_real_10M/ (as a gzip file). This file was obtained by blasting the computed sequence of transcripts (mini_real_10M_tx_seq.fa) against the genome region (chr1_hg38_10M.fa).
Cheers