PureCLIP
PureCLIP copied to clipboard
pureclip error in plant's CLIP-Seq
Hi: I have tried to use pureclip to identify cross linked sites from CLIP-Seq. But when I used bam and fa from Arabidopsis, it always show errors (I have used different chromosome notations: chr1, Chr1, and 1). I have use the same command line in example data successfully. Thanks! Linhua
pureclip -i chr.bam -bai chr.bam.bai -g chrNucTAIR10.fa -o test.bed -nt 10 -iv 'chr1'
Protein-RNA crosslink site detection
===============
Created look-up table for values from -2000 to 0 with step size 0.00333333 (size: 600000).
Loading reference ...
[1] 61202 abort pureclip -i chr.bam -bai chr.bam.bai -g chrNucTAIR10.fa -o test.bed -nt 10 -i