#!/usr/bin/bash
find /Data/.../bam_leafcutter/ -name "*bam" > bam_files.txt
process_bam() {
bamfile=$1
echo Converting $bamfile to $bamfile.junc
samtools index $bamfile
regtools junctions extract -a 8 -m 50 -M 500000 -s XS $bamfile -o $bamfile.junc
echo $bamfile.junc >> test_juncfiles.txt
}
export -f process_bam
parallel --j 20 process_bam {} :::: bam_files.txt
parallel并行运算
于 2024-09-23 10:47:14 首次发布