kma
kma copied to clipboard
Error 11 in running generate_introns.py
Hi, I encountered an error when running generate_introns.py:
(base) [xfu1@sjgenappprdn01 hg38]$ python2 /export/home/xfu1/Profiler/Profiler/kma/inst/pre-process/generate_introns.py --genome /export/home/xfu1/hg38/Homo_sapiens.GRCh38.dna.alt.fa --gtf /export/home/xfu1/hg38/Homo_sapiens.GRCh38.108.gtf --extend 20 --out /export/home/xfu1/Profiler/Profiler/kma/out_dir
INFO: Reading in GTF: /export/home/xfu1/hg38/Homo_sapiens.GRCh38.108.gtf
INFO: Grouping transcripts by gene
INFO: Writing intron BED file: /export/home/xfu1/Profiler/Profiler/kma/out_dir/introns.bed
INFO: Computing intron-to-transcript compatability
INFO: Opening FASTA: /export/home/xfu1/hg38/Homo_sapiens.GRCh38.dna.alt.fa
INFO: Note: will take a while the first time it is opened.
Traceback (most recent call last):
File "/export/home/xfu1/Profiler/Profiler/kma/inst/pre-process/generate_introns.py", line 154, in
I used the human reference genome downloaded from ensembl. Can someone help me to fix it? Thanks in advance! -Eddie