find ~/input/Yunnan15nation/ -print | grep 'bam$' > ~/input/Yunnan_dir.txt #把目标文件名写进dir.txt
cat ~/input/Yunnan_dir.txt |while read line
do
dir_prefix=$(dirname $line) #将函数结果返回给变量的时候用$() 包住函数
midname=$(basename $line .bam)
samtools mpileup -q 1 -d 8000 -f ~/input/GRCh37.75_reference/human_g1k_v37_decoy.fa $line | varscan mpileup2snp --min-coverage 3 --min-reads2 1 --min-avg-qual 20 --p-value 0.99 --min-var-freq 0.01 --output-vcf 1>${dir_prefix}${midname}.varscan.snp.vcf
echo "hava done ${dir_prefix}${midname}" >>"~/bash/Yunnan15.log"
done
Linux对批量文件执行某操作
最后编辑于 :
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 一、待解决问题 A文件如下所示。tip:ID列指的是图中带红框的那列,A大小73.1MB、B大小约1GBimage...
- 一、在任意文件夹下面创建形如 1/2/3/4/5/6/7/8/9 格式的文件夹系列。 二、在创建好的文件夹下面,里...
- 参考基因组下载有三大全文网站提供参考基因组下载,它们分别是:1.NCBI (https://www.ncbi.nl...