Skip to content

Commit

Permalink
Update README.md
Browse files Browse the repository at this point in the history
  • Loading branch information
terrimporter authored Sep 2, 2021
1 parent 2da838e commit 9232489
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion README.md
Original file line number Diff line number Diff line change
Expand Up @@ -502,7 +502,7 @@ ESVtable <- reshape2::dcast(df, SampleName ~ GlobalESV, value.var = "ESVsize", f
If you are targeting a broad group, such as Metazoa using COI primers, you can still filter out pseudogenes using removal method 1 that uses ORFfinder. This can be done in two steps, for example by first processing invertebrate phyla, then processing phylum chordata that includes the vertebrata clade (see NCBI taxonomy). Note that pseudogene removal method 2 that uses HMMer is currently only available for COI arthropoda at this time.

1. Edit the config_ESV.yaml file as follows:
Line 156, set the taxonomy filter to '-e Annelida -e Arthropoda -e Cnidaria -e Echinodermata -e Mollusca -e Nematoda -e Platyhelminthes -e Porifera' to target invertebrate animal phyla
Line 156, set the taxonomy filter to '-e Annelida -e Arthropoda -e Bryozoa -e Cnidaria -e Echinodermata -e Mollusca -e Nematoda -e Nematomorpha -e Nemertea -e Platyhelminthes -e Porifera' to target invertebrate animal phyla
Line 163, set to removal method 1 (uses ORFfinder)
Line 177, set to '5' to use the invertebrate mitochondrial genetic code for translation
Run snakemake. Move invertebrate outfiles into their own directory so they do not get over-written: taxon.zotus, chimera.denoised.nonchimeras.taxon, orf.fasta*, rdp.csv.tmp, results.csv .
Expand Down

0 comments on commit 9232489

Please sign in to comment.