Mercurial > repos > artbio > repenrich2

view test-data/tool_wrapper.sh @ 9:2b61c6407efb draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/main/tools/repenrich2 commit 961a14cd73dd8e40b2f30d5e4df136b98cd8f07a
author artbio
date Thu, 25 Apr 2024 16:22:34 +0000
parents 4905a332a094
children
line wrap: on
line source

input_base=Samp
baseReference=chrM

bowtie-build ${baseReference}.fa ${baseReference}

python RepEnrich_setup.py ${baseReference}_repeatmasker.txt ${baseReference}.fa \
       setup_folder_${baseReference}

bowtie $baseReference -p 16 -t -m 1 -S --max ${input_base}_multimap.fastq \
       ${input_base}.fastq ${input_base}_unique.sam

samtools view -bS ${input_base}_unique.sam > ${input_base}_unique.bam
samtools sort ${input_base}_unique.bam ${input_base}_unique_sorted
mv ${input_base}_unique_sorted.bam ${input_base}_unique.bam
samtools index ${input_base}_unique.bam
rm ${input_base}_unique.sam

python RepEnrich.py ${baseReference}_repeatmasker.txt ${input_base} ${input_base} \
        setup_folder_${baseReference} ${input_base}_multimap.fastq ${input_base}_unique.bam --cpus 16