Hello! I am using ChIPseeker to analyze the BED data from SEACR (a CUT & RUN analysis tool). When I used readPeakFile() function in ChIPseeker to read the BED file, it gave out such a list:
GRanges object with 18965 range and 3 metadata columns:
seqnames ranges strand | V4 V5 V6
<rle> <iranges> <rle> | <numeric> <numeric> <factor>
[1] chr10 42527196-42538077 * | 3500650 3840.71 chr10:42529404-42529405
... According to the description from SEACR (https://github.com/FredHutch/SEACR), V4 is "Total signal contained within denoted coordinates", V5 is "Maximum bedgraph signal attained at any base pair within denoted coordinates". I wonder if V4 and V5 satisfy the requirements of Grange objects for ChIPseeker. And which column was taken as the real score for ChIPseeker, V4 or V5? I am puzzled since the example data from ChIPseeker gave out a list with different metadata columns.
GRanges object with 1331 ranges and 2 metadata columns:
seqnames ranges strand | V4 V5
<rle> <iranges> <rle> | <character> <numeric>
[1] chr1 815093-817883 * | MACSpeak1 295.76
So my last question is, if my V4 and V5 are not suitable, how can I convert it into right order?
Thank you very much!