pizzly
pizzly copied to clipboard
Error, could not find any transcript sequences
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 )