赞
踩
SRR10428360_1_val_1.fq.gz
SRR10428360_2_val_2.fq.gz
index=/home/data/vip30/database/genome/Ensembl/Homo_sapiens/Homo_sapiens.GRCh38_release101.genomeinputdir=/home/data/vip30/data/rnaseq/GSE140202/cleandata/trim_galoreoutdir=/home/data/vip30/data/rnaseq/GSE140202/mapping/Hisacat /home/data/vip30/data/rnaseq/GSE140202/rawdata/SraAccList.txt | while read iddo echo "hisat2 -p 16 -x ${index} -1 ${inputdir}/${id}_1_val_1.fq.gz -2 ${inputdir}/${id}_2_val_2.fq.gz 2>${id}.log | samtools sort -@ 1 -o ${outdir}/${id}.Hisat_aln.sorted.bam - && samtools index ${outdir}/${id}.Hisat_aln.sorted.bam ${outdir}/${id}.Hisat_aln.sorted.bam.bai"done >Hisat.shnohup sh Hisat.sh >Hisat.log &
index:索引
in
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。