|Thread||Thread Starter||Forum||Replies||Last Post|
|processing pindel output files||odoyle81||Bioinformatics||9||01-17-2017 10:10 PM|
|coverageBed -g option error||cmccabe||Bioinformatics||0||09-29-2015 01:46 PM|
|BEDTools coverageBed error||JMorlan||Bioinformatics||0||04-04-2014 08:29 AM|
|coverageBed error: terminate called after throwing an instance of 'std::bad_alloc'||hari_iyer16||Bioinformatics||4||04-11-2013 06:12 PM|
|Yesterday, 12:52 AM||#1|
Join Date: Apr 2017
SplicingCompass: error in processing CoverageBed files
I want to do differential splicing analysis using SplicingCompass (http://www.leibniz-hki.de/files/cont...ssTutorial.pdf).
- bam files (ideally from TopHat mapping)
- counted read numbers using CoverageBed, aligned to all exons in union transcript (script is provided by the authors, annotation file downloaded from UCSC database; hg38, annotation file: genes and gene predictions, NCBI_Refseq, GTF; id to symbol mapping file: genes and gene predictions, NCBI_Refseq, selected fields from primary and related tables --> name and name2)
- junction files from mapping.
I used mapping with STAR, so I converted the SJ_out files to junction.bed files using this script on github:
to convert them into the right format for SplicingCompass.
I had no problems until line 38, where I got the following error message:
Number of input files ok?
Processing CoverageBed output files:
|=============== | 12%
Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec, :
scan() expected 'a real', got 'D00614:183:CABNAANXX:3:1107:7896:66038'
It is unable to process the gff files, I couldn't identify where I made a mistake. Could you help me with it, please?
See the examples attached.
Thank you for the help!
|coveragebed, differential splicing, splicing, splicingcompass|