我如何运行MetaVelvet? “入门”说明无法按预期工作

时间:2015-05-06 17:37:50

标签: bioinformatics

我正在尝试通过在输出上运行MetaVelvet来改进宏基因组程序集。 'Getting Started' instructions for the MetaVelvet extension不生成Graph2文件,这是MetaVelvet扩展所必需的。

如何生成Graph2文件以便运行MetaVelvet?

1 个答案:

答案 0 :(得分:2)

'入门'指令在velvetg命令中缺少一个标志。如果需要生成Graph2文件,则需要该标志。

以下是这些说明的更正版本:

  

〜$ velvetg out-dir -read_trkg yes -exp_cov auto -ins_length 260

提示: 如果你使用多个k-mers运行velveth,因此有多个输出文件夹,你可以使用find命令组合每一个而不必每次都运行命令。

想象一下,多个k-mer velveth输出文件夹以' pooled_reads'开始:

  找到。 -name' pooled_reads _ *' -exec velvetg {} -read_trkg yes -exp_cov auto -ins_length 200 \;

您也可以使用xargs或while循环来执行此操作,但我喜欢find命令。