Biostar Beta. Not for public use.
How to create bigwig using mm10.fasta and mm10.fasta bowtie index files
0
Entering edit mode
13 months ago
United States

Hi,

I am trying to call copy number from mouse data using HMMcopy.

Created the wig files for my BAM samples. Also created mm10_gc.wig using gcCounter (HMMcopy). I am facing issues while creating the big wig file.

generateMap.pl -w 40 -i mm10 mm10.fa > mm10.bigwig

Getting the following error

HMMcopy_0.1.1/util/mappability/internal/fastaToRead missing, try recompiling...

ADD COMMENTlink
1
Entering edit mode

It looks like your HMMcopy install is incomplete. Did you compile it yourself? If the said executable (fastaToRead) is there then try adding that directory you $PATH.

ADD REPLYlink
0
Entering edit mode

Hi Genomax, I reinstalled the tool it worked.

ADD REPLYlink
1
Entering edit mode
12 weeks ago
ATpoint 17k
Germany

Get kentUtils, and use wigTobigWig. It requires a file with the chromosome sizes, which you can get here. Save the chrSizes as a text file. Then use:

./wigToBigWig in.wig chromSizes.txt out.bigwig
ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1