Skip to content

A question about bam file #9

@NJU-Bio-Info

Description

@NJU-Bio-Info

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
  • your bam file:
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?

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions