Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Can not get the fusions with the testdata #31

Open
kingyang728 opened this issue Oct 20, 2020 · 1 comment
Open

Can not get the fusions with the testdata #31

kingyang728 opened this issue Oct 20, 2020 · 1 comment

Comments

@kingyang728
Copy link

kingyang728 commented Oct 20, 2020

Can not detect the fusions with the test data,is anyone can help me? all the testdata is from the git and
http://opengene.org/dataset.html

Here is my shell command:
./genefuse -r Homo_sapiens_assembly19.fasta -f genes/cancer.hg19.csv -1 ./testdata/genefuse.R1.fq.gz -2 ./testdata/genefuse.R2.fq.gz -h report.html > result
12:17:45 start with 4 threads
12:18:1 mapper indexing done
12:18:6 sequence number before filtering: 0
12:18:6 removeByComplexity: 0
12:18:6 removeByDistance: 0
12:18:6 removeIndels: 0
12:18:53 matcher indexing done
12:18:53 removeAlignables: 0
12:18:53 found 0 fusions
12:18:53 done

@kingyang728 kingyang728 changed the title Can not get Can not get the fusions with the testdata Oct 20, 2020
@granek
Copy link

granek commented Jan 6, 2021

I had the same problem! The fusion files in the repo are not compatible with the NCBI version of hg19 that is linked in the README.
The fusion files have chromosome names that start with "chr", but the genome from NCBI just uses the chromosome number (e.g. "chr15" vs "15").

The easiest fix is to download hg19 from UCSC: http://hgdownload.cse.ucsc.edu/goldenPath/hg19/bigZips/hg19.fa.gz

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants