metagenlab / zAMP

zAMP is a bioinformatic pipeline designed for convenient, reproducible and scalable amplicon-based metagenomics
https://zamp.readthedocs.io/en/latest/
MIT License
8 stars 5 forks source link

Updated cutadapt parameters #56

Closed vicasze closed 4 weeks ago

vicasze commented 1 month ago

Summary

I added a new parameter to be able to provide additional cutadapt arguments for the forward and/or reverse primer. This is needed for the ITS database preparation, because the presence of the primer needs to be set as "optional" according to their presence/absence on the DB. On the other hand, I remove this optional primer trimming in zamp run and zamp insilico, because the read data should always contain primers and be trimmed.

In zamp run, cutadapt is called directly on the rules file, and not as an script.

I added usage example to prepare a fungal ITS database.

Details

Additions

Changes

Fixes

farchaab commented 1 month ago

Can you add commit hashes and their description in the Details section ?

vicasze commented 4 weeks ago

@farchaab as suggested, cutadapt is called directly on the snakemake rules files, and not as an outside script. Tested on data/sra-samples.tsv