python - 遺伝子座タグ ID に基づいてシーケンスを取得する - SnakeMake

Question

この種のプログラミングに関する質問が StackOverflow ガイドライン内で許可されることを願っています。そうでない場合は、間違っている場合は修正してください。

私は、バイオインフォマティクスの目的で、snakemake のワークフローに取り組んできました。私が今持っているファイルは、さまざまな遺伝子座タグを含むテキスト ファイルです。これらの遺伝子座タグからシーケンスを取得する必要があります。私はsnakemakeのドキュメントを研究してきましたが、遺伝子座タグと生物に基づいてシーケンスを取得するためのlibaryの正しい機能をまだ見つけることができません...私が今持っているコード:

rule get_id:
    input:
        a="RNA-Seq-counts.txt"
    output:
        b="output.txt"
    shell:
        "python RetrieveFirstColumn.py {input} {output}"

output.txt の出力

lp_0001
lp_0002
lp_0004
lp_0005
lp_0009
lp_0010
lp_0011
lp_0012
lp_0013
lp_0014
lp_0015
lp_0016
lp_0017