Skip to content

Snakemake pipeline to produce DRAM product refined based on temporal paper metabolic definitions

Notifications You must be signed in to change notification settings

emerge-bii/DRAM_product_refined_pipeline

 
 

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

42 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

DRAM product refined pipeline

Description

  • Input: DRAM annotations file.
  • Output: distillate (gene presence/absence), product (reaction % presence) and product refined (pathway calls).

Each pathway (e.g. Calvin-Benson carbon fixation) is defined as a combination of reactions (e.g. RubisCO). For a genome to have a positive pathway call, it must have >60% of the total pathway (accounting for alternate reaction paths), it must have some presence in >70% of the reactions (at least one gene) and, it must have the signature genes, where applicable (used to distinguish pathways with high overlap, e.g. the Wood-Ljungdahl methanogenic pathway and hydrogenotrophic methanogenesis).

How to Run

  • Create conda environment: conda env create --name DRAM_product_refined -f environment.yaml.
  • Activate environment: conda activate DRAM_product_refined.
  • Add path to DRAM annotations to the config in config/config.yaml.
  • Run pipeline: snakemake -j 30.

About

Snakemake pipeline to produce DRAM product refined based on temporal paper metabolic definitions

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages

  • Python 96.9%
  • HTML 3.1%