在docker中用cleandata与宿主序列比对得到sam文件这一步出现错误,


输入的命令是for i in `tail -n+2 $datafile/metadata.txt|cut -f1`;do bowtie2 --threads 6 \ -x $host_index/Canis_lupus/wolf \ -1 $workdir/1.host_del_and_qc/cleandata/${i}_1.fastq \ -2 $workdir/1.host_del_and_qc/cleandata/${i}_2.fastq \ -S $workdir/1.host_del_and_qc/cleandata/${i}.bowtie.sam \ --very-sensitive-local \ 2>$workdir/1.host_del_and_qc/cleandata/${i}.bowtie.map.log  # Filter and output as SAM(or BAM),-f 4 to Remove reads from the alignment,  # if you want to keep the -F on the alignment (12 for both ends), add -b before the output to output in BAM format samtools view --threads 6 \ -f 12 -h \ $workdir/1.host_del_and_qc/cleandata/${i}.bowtie.sam \ > $workdir/1.host_del_and_qc/cleandata/${i}.fin.sam done

得到的结果是无法打开attachments-2025-03-U1HS1aia67d3c5a02b0d3.png

请先 登录 后评论

1 个回答

xun - 电路元件工程师

你的上一步,bowtie应该报错了,你看看这个文件,
$workdir/1.host_del_and_qc/cleandata/${i}.bowtie.map.log
里面咋说的,是不是你的宿主文件写的不对?


请先 登录 后评论
  • 1 关注
  • 0 收藏,40 浏览
  • ZXX111 提出于 7小时前

相似问题