Skip to content
Snippets Groups Projects
Commit 545e8186 authored by Shaman Narayanasamy's avatar Shaman Narayanasamy
Browse files
parents 9d32ae2e 609c21ae
No related branches found
No related tags found
No related merge requests found
......@@ -137,19 +137,19 @@ rule ANALYSIS_MG_VARIANT_CALLING:
echo "Indexing bam: {input[1]}"
samtools index {input[1]}
fi
if [[ ! -f {input[0]}.fai ]]
then
echo "Fasta index doesn't exist, Creating one..."
echo "Indexing fasta: {input[0]}"
samtools faidx {input[0]}
fi
#temporary directory and files
VCF_MPU=$(mktemp --tmpdir={TMPDIR} -t "XXXXXX.mpu.vcf")
VCF_FRB=$(mktemp --tmpdir={TMPDIR} -t "XXXXXX.frb.vcf")
VCF_PLT=$(mktemp --tmpdir={TMPDIR} -t "XXXXXX.plt.vcf")
### run_mpileup {input[0]} {input[1]} {output[1]}
echo "Running samtools mpileup"
samtools mpileup -uf {input[0]} {input[1]} |\
......@@ -208,19 +208,19 @@ rule ANALYSIS_MT_VARIANT_CALLING:
echo "Indexing bam: {input[1]}"
samtools index {input[1]}
fi
if [[ ! -f {input[0]}.fai ]]
then
echo "Fasta index doesn't exist, Creating one..."
echo "Indexing fasta: {input[0]}"
samtools faidx {input[0]}
fi
#temporary directory and files
VCF_MPU=$(mktemp --tmpdir={TMPDIR} -t "XXXXXX.mpu.vcf")
VCF_FRB=$(mktemp --tmpdir={TMPDIR} -t "XXXXXX.frb.vcf")
VCF_PLT=$(mktemp --tmpdir={TMPDIR} -t "XXXXXX.plt.vcf")
### run_mpileup {input[0]} {input[1]} {output[1]}
echo "Running samtools mpileup"
samtools mpileup -uf {input[0]} {input[1]} |\
......
......@@ -220,6 +220,18 @@ rule CHECK_TOOL_VERSION:
# python {SRCDIR}/make.pwy.hierarchy.kegg.py -o {output[1]}
# """
rule _DOWNLOAD_TRIMMOMATIC_ADAPTERS:
output:
"{DBPATH}/adapters/adapters.done"
shell:
"""
wget http://www.usadellab.org/cms/uploads/supplementary/Trimmomatic/Trimmomatic-Src-0.32.zip
unzip Trimmomatic-Src-0.32.zip
cp -r trimmomatic-0.32/adapters {DBPATH}
rm Trimmomatic-Src-0.32.zip && rm -rf trimmomatic-0.32
touch {output}
"""
rule _DOWNLOAD_KEGG_INFORMATION:
log:
U_LOG
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment