Hi
Newbie here so please forgive me if this is a silly question.
I'm trying to use the Viral Fusion Seq (VFS) program to determine if there are any fusion products between the human genome and a herpes virus.
I got the VFS installed and running correctly according to their manual. But they don't give a description of what to do with your own data sets.
I have fastq files from human tissue samples (3) and I have the herpes virus genome.
Here are what I believe are the relevant settings in my vfs.conf file:
I "prepped" my human fastq files by running them through the prep.fastq.2.Illumina.1.8.gz.pl script. This produced a file called hs.vfs_in.fq.gz.
I tried to run VFS as follows:
Unfortunately I don't get any fusion products predicted. I'm not expecting 1000's of fusion product sequences. But the fact that I got nothing whatsoever makes me think I messed up some how.
Can anyone point me in the right direction? Are there other options in the vfs.conf file that I need to set/adjust?
Any and all help is welcome.
Thanks
Newbie here so please forgive me if this is a silly question.
I'm trying to use the Viral Fusion Seq (VFS) program to determine if there are any fusion products between the human genome and a herpes virus.
I got the VFS installed and running correctly according to their manual. But they don't give a description of what to do with your own data sets.
I have fastq files from human tissue samples (3) and I have the herpes virus genome.
Here are what I believe are the relevant settings in my vfs.conf file:
Code:
viralFa = /data/rna/vzv/vzv_genome.fa humanFa = /usr/local/vfs/references/human/hg38.fa humanDecoy = /usr/local/vfs/references/human/hs37d5.fa ntDB = /usr/local/vfs/annotation/nt/nt
I tried to run VFS as follows:
Code:
viral.fusion.pl hs_vzv /abs/path/to/hs.vfs_in.fq.gz
Can anyone point me in the right direction? Are there other options in the vfs.conf file that I need to set/adjust?
Any and all help is welcome.
Thanks