文件过滤_hisat2比对小不同:fastp过滤后和trim_galore过滤后文件

trim_galore过滤后的文件名:

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

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值