Hello, How do I remove a factor level for the sequence names of a GRanges object?
Sometimes I work with data sets that have not been aligned to all the same chromosomes. To make them comparable I have to remove or ignore reads for the odd chromosome. So if I have a GRanges object for which I want to remove all reads matching a given chromosome, e.g.: gr <- gr[seqnames(gr) != "chrXHet"] levels(seqnames(gr)) will return all original seqnames, but I would like to exclude "chrXHet" since I've removed all the reads matching that chromosome. How do I do this? -Chris _______________________________________________ Bioc-sig-sequencing mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
