Skip to content

07_ANNOTATION

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

Run HMMs on for single-copy core genes for each sample contig set.

for i in `cat AMOR_2020_Good`; do anvi-run-hmms -c 03_INDIV_ASSEMBLY/${i}/${i}.prefixed.contigs.db -T 50; done

Gene annotation (on kjempefuru)

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