u014210048的专栏

编译技巧

samtools flags 1294$samtools flags 12940x50e 1294 PROPER_PAIR,UNMAP,MUNMAP,SECONDARY,DUP$samtools flags 2560x100 256 SECONDARY$samtools flags MUNMAP0x8 8 MUNMAP$samtools flags UNMAP0x4 4 UNMAP$samtools fl

msisensor-pro: error while loading shared libraries: libhts.so.3: cannot open shared object file: No such file or directory在测试msisensor-pro时候，发现libhts.so.3找不到。判断libhts.so.3应该属于htslib库。于是安装一个。https://leoatchina.github.io/2017/08/26/install_samtool

$conda info --envs# conda environments:#base * /home/wubin/miniconda3DECoN /home/wubin/miniconda3/envs/DECoNSnakeMakeEnv /home/wubin/miniconda3/envs/SnakeMakeEnvcnvkit /home/wubin/mi.

今天被自己坑了一把。samtools index a.bam b.bam上面这样写，b.bam将是a.bam的bai文件。如果写作samtools index a.bam则会生成一个a.bam.bai文件，这个a.bam.bai文件的md5值会与b.bam一模一样

建立一个a.pl#!/usr/bin/perl -wuse strict;foreach my $i(0..10){ print "$i\n"; sleep(1);}建立一个b.pl#!/usr/bin/perl -wuse strict;foreach my $i(0..10){ print "$i==b\n"; sleep(1);}建立一个a.shperl a.pl >a.o 2>.