#提取比对上的reads的ID,根据ID从原始测序数据中提取这些reads cut -f1 $line.TDNA.sam |sort|uniq > $line.TDNA.ID seqtk subseq $read1 $line.TDNA.ID > ${line}_1.TDNA.fq seqtk subseq $read2 $line.TDNA.ID > ${line}_2.TDNA.fq #spades组装到contig水平,contigs.fasta即为最终结果 spade...