viralintegration icon indicating copy to clipboard operation
viralintegration copied to clipboard

Add MetaPhlAn 3.0 module

Open allaway opened this issue 2 years ago • 1 comments

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.

allaway avatar Apr 22 '22 18:04 allaway