Skip to content

07_ANNOTATION

Jump to bottom
eolesin edited this page Jun 1, 2021 · 15 revisions

Gene annotation

cd 03_INDIV_ASSEMBLY/
for SET in `cat /export/dahlefs/work/Metagenomes_chimneys_2020_workfolder/AMOR_2020_Good`; \
    do anvi-get-sequences-for-gene-calls -c $SET/$SET.prefixed.contigs.db \
    --get-aa-sequences \ 
    -o $SET/$SET.prefixed.aminoacidseqs.fa; \
done

cd ..

# run kegg - kofams using the dev version of Anvi'o from Achim's anvio-dev conda env.

for SET in `cat /export/dahlefs/work/Metagenomes_chimneys_2020_workfolder/AMOR_2020_Good`; 
    do anvi-run-kegg-kofams -c 03_INDIV_ASSEMBLY/$SET/$SET.prefixed.contigs.db \
    --kegg-data-dir /export/dahlefs/work/Databases/2021-04-22_KOFAM-DB -T 40 \
    >> /export/dahlefs/work/Metagenomes_chimneys_2020_workfolder/06_ANNOTATION/2021-05-26_anvi-run-kegg-kofams.log;
done

# Estimate metabolism

Clone this wiki locally