Hello, such an amazing GRO-seq mediated TU identification tool. When I used groHMM, I found that you bam file is different from mine so I received an error log:
Error in Fp[[i]] + 1 : non-numeric argument to binary operator
- my bam file format after import into R:
GRanges object with 6 ranges and 0 metadata columns:
seqnames ranges strand
<Rle> <IRanges> <Rle>
[1] chr1 13934-13984 -
[2] chr1 13935-13985 -
[3] chr1 13939-13989 -
[4] chr1 13948-13986 -
[5] chr1 13949-13999 -
[6] chr1 13951-13992 -
-------
seqinfo: 195 sequences from an unspecified genome
GRanges object with 210615 ranges and 0 metadata columns:
seqnames ranges strand
<Rle> <IRanges> <Rle>
[1] chr7 5936085 -
[2] chr7 23221745 +
[3] chr7 150801773 +
[4] chr7 101909724 +
[5] chr7 8429817 +
... ... ... ...
[210611] chr7 157011485 +
[210612] chr7 91477417 -
[210613] chr7 30639575 +
[210614] chr7 42202754 -
[210615] chr7 138754360 -
-------
seqinfo: 1 sequence from an unspecified genome
Have noticed the difference in "ranges" column, so do you mean that the ranges should be points rather than range?
Hello, such an amazing GRO-seq mediated TU identification tool. When I used groHMM, I found that you bam file is different from mine so I received an error log:
Have noticed the difference in "ranges" column, so do you mean that the ranges should be points rather than range?