rseg segmentation fault core dumped

68 views
Skip to first unread message

cc1...@essex.ac.uk

unread,
Jul 15, 2017, 6:38:42 AM7/15/17
to RSEG Users
I am interested in using your package- 'rseg' for differential peak calling between treated and control read files that I have.
I seem to be having some trouble with it:
Although I'm told that everything installed fine it doesn't seem to work:
I have two sam files (which I converted to bed files) 
and I also have a file that specifies the size of the chromosomes as follows:

chrom   size
chr1    197195432
chr2    181748087
chr3    159599783
chr4    155630120
chr5    152537259
chr6    149517037
chr7    152524553
chr8    131738871
chr9    124076172
chr10   129993255
chr11   121843856
chr12   121257530
chr13   120284312
chr14   125194864
chr15   103494974
chr16   98319150
chr17   95272651
chr18   90772031
chr19   61342430
chrX    166650296
chrY    15902555
chrM    16299
chr13_random    400311
chr16_random    3994
chr17_random    628739
chr1_random     1231697
chr3_random     41899
chr4_random     160594
chr5_random     357350
chr7_random     362490
 
chr8_random     849593

NOTE: I should mention that in the documentation, it is not at all clear how the actual input files are specified i.e. the treated and control BED files. Should they be precluded by '-input' and -control'?

However while everything seems to be in order I get an error:

rseg-diff -c ~/mm9.genome -o out.bed -i 20 -v -mode 2  test.bed control.bed  #-c the last two arguments are the input treated and control files

Error:
[LOADING_DATA] chromosomes
[LOADING_DATA] reads
Segmentation fault (core dumped)
Reply all
Reply to author
Forward
0 new messages