How to generate input wig files for HMMcopy R package?
1
2
Entering edit mode
12 months ago

Hello everyone, I am trying to use the HMMcopy R package for low pass copy number analysis. HMMcopy needs .wig files as input for the analysis. Unfortunately, there is no documentation available about how to obtain these files. Can someone please provide any info about how to generate those wig files for HMMcopy. I have fastq files and i have already generated bam and sam files for those fastq files. Thanks in advance.

R HMMcopy Copy number analysis genome Tool • 486 views
ADD COMMENT
0
Entering edit mode
11 months ago
Alex • 0

Not sure if you figured it out, but I have been looking around on how to do this as well and found this:

http://bioinformatics-ca.github.io/bioinformatics_for_cancer_genomics_mod5_lab_data_prep_2015/

I haven't had the chance to try it yet, but while looking for how to do it I stumbled across your question. I think what we are interested in is towards the bottom of the page for the GC and mapability files. Towards the top is the readCounter for the third wig file needed.

Hope this helps!

ADD COMMENT

Login before adding your answer.

Traffic: 938 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6