viralintegration
viralintegration copied to clipboard
Add MetaPhlAn 3.0 module
Description of feature
Last summer, @parnaljoshi came up with the idea of using MetaPhlAn to look for viral reads in tumor data. The general approach was:
- align tumor fastq to human genome with STAR
- save unmapped reads
- run unmapped reads through MetaPhlAn:
metaphlan $out_prefix"Unmapped.out.mate1",$out_prefix"Unmapped.out.mate2" --bowtie2out metagenome.bowtie2.bz2 --add_viruses --ignore_eukaryotes --nproc 16 --input_type fastq --read_min_len 30 > "analysis/star_results/"$substring.txt
First, we need to simply add the metaphlan3 module.