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.
1 answer
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!
Log in to answer this question.
-
feature_ref files in upstream analyses related to Perturb-seq
written by 2443422660 •
-
Identifying primer sequences from raw FASTQ files
written by Chandini •
-
Percentage of samples covered per region in IGV
written by monica •
-
what is the output of 'samtools coverage'
written by monica •
-
Diversity analysis
written by anjalivasudevan98 •
-
fastqc html file
written by Cihan •
-
how to generate .bed files for methylation (RRBS) analysis
written by arshad1292 •
-
Problems with fastq-dump generating a single fastq file
written by pomodoro_sinensis •
-
Cellranger cnv analysis bin size
written by enho •
-
Do downloaded fastq files from NCBI (SRA) needs preprocessing before used in bowtie2 and featureCo…
written by RS •
