handle many SAM files （批量处理SAM文件）

发表于2019年11月6日由daizao 
#sam2bam,bam2sort,bam2index
screen -S yourname
for i in `ls *.sam`;
do
j=${i/\.sam/};
echo "samtools view -b -S $j.sam > $j.bam
samtools sort $j.bam sorted_$j
samtools index sorted_$j.bam" >> $j.sh
done

for i in `ls *.sh`;do echo "bash $i" >> total.sh;done

ParaFly -c total.sh -CPU 9
此条目发表在Linux分类目录。将固定链接加入收藏夹。

发表评论

邮箱地址不会被公开。 必填项已用*标注

To create code blocks or other preformatted text, indent by four spaces:

    This will be displayed in a monospaced font. The first four 
    spaces will be stripped off, but all other whitespace
    will be preserved.
    
    Markdown is turned off in code blocks:
     [This is not a link](http://example.com)

To create not a block, but an inline code span, use backticks: 

Here is some inline `code`.

For more help see http://daringfireball.net/projects/markdown/syntax

Protected with IP Blacklist CloudIP Blacklist Cloud