GeneFuse icon indicating copy to clipboard operation
GeneFuse copied to clipboard

No fusion generated using demo fastq files

Open foxchase opened this issue 7 years ago • 4 comments

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

foxchase avatar Mar 05 '18 15:03 foxchase

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

sfchen avatar Mar 06 '18 02:03 sfchen

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

foxchase avatar Mar 06 '18 21:03 foxchase

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.

sfchen avatar Mar 07 '18 02:03 sfchen

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

foxchase avatar Mar 08 '18 19:03 foxchase