Question: Run CIRI program directly from BWA-mem alignment
0
gravatar for onjinakashima
3 days ago by
onjinakashima0 wrote:

Hello, I'm trying to find circular ARN using the program CIRI2. CIRI2 takes as input SAM file from BWA-mem, but I would like to go straight to CIRI2 output without saving SAM file. Is there a way to achieve this?

rna-seq alignment • 97 views
ADD COMMENTlink written 3 days ago by onjinakashima0
1

You could try process substitution for the input sam perl CIRI.pl -I <(put your bwa mem call (or even bwa mem | samtools pipeline) here) -O outfile -F ref.fa (-R ref_dir/) -P. It is not guaranteed to work, though, especially if CIRI makes multiple passes over the sam file.

ADD REPLYlink written 3 days ago by cschu181670
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1266 users visited in the last hour