wigToBigWig: items longer than chrom.size,how to solve it? please help me !
1
0
Entering edit mode
2.6 years ago

1.get GSM2258884: from https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE85158

2.Obtained after decompression: 13.TNBC-Basal.MB468.H3K27ac.rep1.wig

3.Remove track:grep -v track 13.TNBC-Basal.MB468.H3K27ac.rep1.wig >13.TNBC-Basal.MB468.H3K27ac.rep1-out.wig

4.wigToBigWig 13.TNBC-Basal.MB468.H3K27ac.rep1-out.wig hg19.chrom.sizes 13.TNBC-Basal.MB468.H3K27ac.rep1-out.bw

showing: line 279589177 of 13.TNBC-Basal.MB468.H3K27ac.rep1-out.wig: chromosome chrM has 16571 bases, but item ends at 16580

UCSC • 1.1k views
ADD COMMENT
0
Entering edit mode

Perhaps the data uses hg38 not hg19?

ADD REPLY
0
Entering edit mode

The hg19 is correct. In addition, I have used hg38, but showing the same problem.

ADD REPLY
1
Entering edit mode

It it only a chrM problem? Quick-n-dirty, modify the chromsizes file and set chrM to 16580, see if that works. Maybe some oddity with the chrM in that particular reference fasta they used? See if things make sense on genome browser of hg19.

ADD REPLY
0
Entering edit mode
2.6 years ago

Use the bedClip UCSC utility to cut the features down to chrom sizes. But double check before that that you're really on the correct assembly and that you have the correct file.

For hg19, for example. the UCSC browser now has both the original chrM sequence (which they picked way back before NCBI had picked a different reference chrM sequence) and the official mitochondria reference called "chrMT" that was selected by the GRCh. So maybe replacing "chrM" with "chrMT" will fix your problem here.

It's possible that this particular file was made with a genome assembly downloaded from e.g. NCBI or Ensembl and someone added the "chr" to the file, which could lead to this particular problem, a length difference on the mitochondrial sequence.

For many applications, you can also remove all the chrM annotations. For a lot of assays, e.g. RNA-seq and most Chip-seq experiments, the biologists don't care about the mitochondrial genome results so they can be removed.

ADD COMMENT

Login before adding your answer.

Traffic: 1837 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