你的上一步,bowtie应该报错了,你看看这个文件,
$workdir/1.host_del_and_qc/cleandata/${i}.bowtie.map.log
里面咋说的,是不是你的宿主文件写的不对?
输入的命令是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
得到的结果是无法打开