pizzly icon indicating copy to clipboard operation
pizzly copied to clipboard

Error, could not find any transcript sequences

Open ndaniel opened this issue 6 years ago • 5 comments

Hello,

when using Pizzly 0.37.3 (SeqAn 2.2.0) and Kallisto 0.43.1 with Ensembl 81 and one gets this error message from Pizzly:

Error, could not find any transcript sequences check that the ids in the FASTA file and GTF file match

when running this:

wget ftp://ftp.ensembl.org/pub/release-81/fasta/homo_sapiens/cdna/Homo_sapiens.GRCh38.cdna.all.fa.gz

wget ftp://ftp.ensembl.org/pub/release-81/gtf/homo_sapiens/Homo_sapiens.GRCh38.81.gtf.gz

wget https://sourceforge.net/projects/fusioncatcher/files/test/reads_1.fq.gz

wget https://sourceforge.net/projects/fusioncatcher/files/test/reads_2.fq.gz

kallisto index -k 31 -i transcriptome.idx  Homo_sapiens.GRCh38.cdna.all.fa.gz

kallisto quant -i transcriptome.idx  -o output_kallisto --fusion  reads_1.fq.gz reads_2.fq.gz

pizzly -k 31 --gtf Homo_sapiens.GRCh38.81.gtf.gz --cache cache.txt  --align-score 2 --insert-size 400 --fasta Homo_sapiens.GRCh38.cdna.all.fa.gz --output output_pizzly output_kallisto/fusion.txt

which worked just fine with: pizzly version: 0.37.1, SeqAn version: 2.2.0, kallisto 0.43.1 (as shown here: https://github.com/pmelsted/pizzly/issues/7 )

ndaniel avatar Aug 14 '17 06:08 ndaniel