Question: Run CIRI program directly from BWA-mem alignment
gravatar for onjinakashima
8 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 • 344 views
ADD COMMENTlink written 8 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 8 months ago by cschu1811.4k
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: 1327 users visited in the last hour