Question: Run CIRI program directly from BWA-mem alignment
gravatar for onjinakashima
5 months 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 • 264 views
ADD COMMENTlink written 5 months ago by onjinakashima0

You could try process substitution for the input sam perl -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 5 months ago by cschu1811.2k
Please log in to add an answer.


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