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

No fusion generated using demo fastq files #9

Open
foxchase opened this issue Mar 5, 2018 · 4 comments
Open

No fusion generated using demo fastq files #9

foxchase opened this issue Mar 5, 2018 · 4 comments

Comments

@foxchase
Copy link

foxchase commented Mar 5, 2018

No fusion detected with the demo data
genefuse -r hg19.fasta -f genes/druggable.hg19.csv -1 R1.fq -2 R2.fq -h report.html > result

@sfchen
Copy link
Member

sfchen commented Mar 6, 2018

I have updated the demo dataset, which contains a valid EML4-ALK fusion with ultra-high frequency.

Please download it again from: http://opengene.org/dataset.html

@foxchase
Copy link
Author

foxchase commented Mar 6, 2018

Hi Shifu,
I downloaded the new demo data and tried again with no fusion generated, the same issue. Additionally, I tried the gen_fusion_file.jl file and got a lot errors, please see attachments. Would you please check the version that you use is the same as that on the github. Many Thanks.
error_Julia.txt
TRFPgenes.txt

@sfchen
Copy link
Member

sfchen commented Mar 7, 2018

Could you let me know how many lines do R1.fq and R2.fq have? They should be 15000000 lines for each file.

For the second issue. Which genes do you want to include? I can generate the CSV file for you. Julia language is evolving so fast, so the package may encounter compatible problems.

@foxchase
Copy link
Author

foxchase commented Mar 8, 2018

Hi Sifu,
Yes, R1 and R2 each has 15,000,000 lines. I ran fusion analysis for your R1 and R2, and there are ALK-EML4 fusion and others in the R1/R2 fastq files; that is no problem for the R1 and R2 fastq file. Is the version of GeneFuse that you use different from the one you posted?
I work on illumina RNA Fusion Panel, attached is the gene list of the Panel, I need hg19 for the CSV file. It is greatly appreciated for doing that. Additionally, please let me know if you need test (positive) fastq files of that panel.
TRFPgenes.txt

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