Question: Run CIRI program directly from BWA-mem alignment
0
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
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 8 months ago by cschu1811.4k
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: 1327 users visited in the last hour