Hi, I'm trying to convert some wig tracks to bigWig format. Actually I have some issues with wigToBigWig. First of all I'm not able to use it piping the output of a previous command. Since I have to remove "track" lines, it would be great to do something like this:
grep -v track MyWig.wig | wigToBigWig - chromsizesfile output.bw I guess piping would be nice to merge several chromosome-specific wig to a single stream for bigWig conversion. Also, some tools that generate wig files apparently add some points beyond the end of chromosomes. Tipically this happens when a stepsize is specified. let's say chrA is 103 bp long, the wig gives point with a span of 5 bp, if the last point of this chromosome is, say: chrA 100 8 wigToBigWig apparently extends the data to position 105 (which is beyond 103...). Wouldn't be better to cut data (with a warning message) instead of giving up with an error? Thanks d /* Davide Cittaro Cogentech - Consortium for Genomic Technologies via adamello, 16 20139 Milano Italy tel.: +39(02)574303007 e-mail: [email protected] */ _______________________________________________ Genome maillist - [email protected] https://lists.soe.ucsc.edu/mailman/listinfo/genome
