#samtools建索引,提取比对上的信息,保存sam格式 samtools index -@ 12 $line.sorted.bam samtools view $line.sorted.bam TDNA > $line.TDNA.sam #提取比对上的reads的ID,根据ID从原始测序数据中提取这些reads cut -f1 $line.TDNA.sam |sort|uniq > $line.TDNA.ID seqtk subseq $read1 $line.TDNA.ID...