time megahit -t 18 -1 tail -n+2 metadatatxtcut -f1sed s^kneaddata;s$_R1_kneaddata_paired_1fastqtr n sed s$ -2 tail -n+2 metadatatxtcut -f1sed s^kneaddata;s$_R1_kneaddata_paired_2fastqtr n sed s$ -o mu
这是一个使用time megahit命令进行多样本组装的命令。
命令解读如下:
time megahit -t 18
:运行time megahit命令,使用18个线程进行组装。-1
:指定输入的第一组PE(paired-end)fastq文件。tail -n+2 metadata.txt
:从metadata.txt文件中提取第2行及之后的内容。cut -f1
:根据字段提取第1列。sed 's/^/kneaddata\//;s/$/_R1_kneaddata_paired_1.fastq/'
:在每行开头加上"kneaddata/",并在末尾加上"_R1_kneaddata_paired_1.fastq"。tr '\n' ','
:将每行以逗号分隔。sed 's/,$//'
:删除最后一个逗号。
-2
:指定输入的第二组PE fastq文件,处理方法与上述相同。-o mulltimegahit
:指定输出目录为mulltimegahit。
综上所述,该命令的作用是使用time megahit命令对多个样本的PE fastq文件进行组装,并将结果输出到mulltimegahit目录中
原文地址: http://www.cveoy.top/t/topic/iuNk 著作权归作者所有。请勿转载和采集!