tourmaline
tourmaline copied to clipboard
Amplicon sequence processing workflow using QIIME 2 and Snakemake
Development branch is here: https://github.com/aomlomics/tourmaline/tree/develop Planning doc for Tourmaline 2.0 is here: https://docs.google.com/document/d/1A7rriQ-hwhYD5ZguarRtQAKJDl_MvyfUCa0e00mSxSs/edit @ksilnoaa
Update [fastqc_per_base_sequence_quality_dropoff.py](https://github.com/aomlomics/tourmaline/blob/develop/scripts/fastqc_per_base_sequence_quality_dropoff.py) to work with fastq_summary.qzv output
- Trimming stats - \# of reads per sample - Sequence quality dropoff
Have a script that will take non-demultiplexed fastq files and demultiplex them by sample barcodes and/or marker (e.g., from a big NextSeq run).
DADA2 in R provides this and it's very convenient to see what each step is doing
Instead of cloning the github directory every time to start a brand new Tourmaline run, have a Bash script that will generate all of the required input config files.
If user is running Rosetta (to enable QIIME 2 to run on Apple Silicon Macs), the instructions to install multiqc conda environment will not work. User first needs to create...
- Option to remove samples with a file and/or below a certain read depth - Option to filter ASVs by length, minimum frequency - Is there a way to cache...