Wig file obtained from MACS14 is longer than chromosome size?
The wig file:
variableStep chrom=chr15 span=1
...
12968274 1
But my chromosome size is chr15 12968165. Wig file is longer than chromosome?
The parameters:
./MACS14 -t ...bed -c ...bed -s 49 -g 1.1e+9 -p 1e-5 --nomodel -w -S --space 1 -m 10,30
Could someone give me a hand?
• 3,644 views
•
link
updated
by
Ram
•
written
by
secretloong •
•
written
by
secretloong •
2 answers
I've run into the same problem. I can't explain why MACS produces wig files like this, but I have a quick fix. I used bedClip to trim off the ends of the wig files.
There is a nice brief blog post written by Sukhdeep Singh on how to use bedClip here.
Yes, MACS extend the reads to 300bp, so it will exceed the limits of the chromosome. bedgraph files generated from it also has the same problem. If you upload to UCSC, it will complain errors. As suggested, use bedClip first
• 0 views
•
link
updated
by
Ram
•
written
by
Ming Tommy Tang
•
written
by
Ming Tommy Tang
Log in to answer this question.
More posts like this
-
How to filter bam file for only reads that span certain regions of the genome?
written by kaufman.eli.j •
-
Low heterozygosity rates in simulated reads
written by all.lirik •
-
ATAC-Seq Data looking like RNA-Seq data...?
written by katywithawhy11 •
-
DNAcopy plot error
written by sunilpachakar •
-
Segmentation fault from wigToBigWig when using the -clip option
written by bowensamigger •
-
gc content file for human genome large and difficult to parse into separate chromosome files
written by chrisclarkson100
-
Converting sra to bed files inefficiently
written by chrisclarkson100
-
ChIP-seq peak calling & visualized in IGV
written by izzy.yichao.cai
-
How To Extract Sequence From Overlapping Regions From Fasta Files With Perl
written by speedgummi •
-
How Do I Get The Gene Names From The Chromosome Locations In A Wig File
written by pixie@bioinfo
