Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Gonzalo S Nido
wgbs-pipeline
Commits
f91f0e2d
Commit
f91f0e2d
authored
Aug 25, 2020
by
Gonzalo S Nido
Browse files
Added pigz instead of gzip
parent
5d46a9f4
Changes
2
Hide whitespace changes
Inline
Side-by-side
prepare_vm.sh
View file @
f91f0e2d
...
...
@@ -4,7 +4,7 @@ echo "#### Updating Ubuntu dist..."
sudo
apt update
&&
sudo
apt
-y
upgrade
echo
"#### Installing some Ubuntu packages..."
sudo
apt
-y
install
unzip build-essential libncurses5-dev zlib1g-dev libbz2-dev liblzma-dev p7zip-full
sudo
apt
-y
install
unzip build-essential libncurses5-dev zlib1g-dev libbz2-dev liblzma-dev p7zip-full
pigz
echo
"#### Installing samtools 1.10..."
cd
~
...
...
wgbs-pipeline.sh
View file @
f91f0e2d
...
...
@@ -307,21 +307,21 @@ fi
echo
-ne
" 1.5 Extracting MT and Lambda methylation... "
CMD
=
"cat
${
output_dir
}
/CpG_context_b38_
${
sample_id
}
.dedup.txt | python3
${
ext_chr_srt
}
${
MT_chr_name
}
| gz
ip
-nc >
\
CMD
=
"cat
${
output_dir
}
/CpG_context_b38_
${
sample_id
}
.dedup.txt | python3
${
ext_chr_srt
}
${
MT_chr_name
}
|
pi
gz -nc >
\
${
output_dir
}
/CpG_context_chrM_
${
sample_id
}
.dedup.txt.gz &&
\
cat
${
output_dir
}
/CpG_context_b38_
${
sample_id
}
.dedup.txt | python3
${
ext_chr_srt
}
${
LAMBDA_chr_name
}
| gz
ip
-nc >
\
cat
${
output_dir
}
/CpG_context_b38_
${
sample_id
}
.dedup.txt | python3
${
ext_chr_srt
}
${
LAMBDA_chr_name
}
|
pi
gz -nc >
\
${
output_dir
}
/CpG_context_LAMBDA_
${
sample_id
}
.dedup.txt.gz &&
\
gz
ip
${
output_dir
}
/CpG_context_b38_
${
sample_id
}
.dedup.txt
\n\
cat
${
output_dir
}
/CHH_context_b38_
${
sample_id
}
.dedup.txt | python3
${
ext_chr_srt
}
${
MT_chr_name
}
| gz
ip
-nc >
\
pi
gz
${
output_dir
}
/CpG_context_b38_
${
sample_id
}
.dedup.txt
\n\
cat
${
output_dir
}
/CHH_context_b38_
${
sample_id
}
.dedup.txt | python3
${
ext_chr_srt
}
${
MT_chr_name
}
|
pi
gz -nc >
\
${
output_dir
}
/CHH_context_chrM_
${
sample_id
}
.dedup.txt.gz &&
\
cat
${
output_dir
}
/CHH_context_b38_
${
sample_id
}
.dedup.txt | python3
${
ext_chr_srt
}
${
LAMBDA_chr_name
}
| gz
ip
-nc >
\
cat
${
output_dir
}
/CHH_context_b38_
${
sample_id
}
.dedup.txt | python3
${
ext_chr_srt
}
${
LAMBDA_chr_name
}
|
pi
gz -nc >
\
${
output_dir
}
/CHH_context_LAMBDA_
${
sample_id
}
.dedup.txt.gz &&
\
gz
ip
${
output_dir
}
/CHH_context_b38_
${
sample_id
}
.dedup.txt
\n\
cat
${
output_dir
}
/CHG_context_b38_
${
sample_id
}
.dedup.txt | python3
${
ext_chr_srt
}
${
MT_chr_name
}
| gz
ip
-nc >
\
pi
gz
${
output_dir
}
/CHH_context_b38_
${
sample_id
}
.dedup.txt
\n\
cat
${
output_dir
}
/CHG_context_b38_
${
sample_id
}
.dedup.txt | python3
${
ext_chr_srt
}
${
MT_chr_name
}
|
pi
gz -nc >
\
${
output_dir
}
/CHG_context_chrM_
${
sample_id
}
.dedup.txt.gz &&
\
cat
${
output_dir
}
/CHG_context_b38_
${
sample_id
}
.dedup.txt | python3
${
ext_chr_srt
}
${
LAMBDA_chr_name
}
| gz
ip
-nc >
\
cat
${
output_dir
}
/CHG_context_b38_
${
sample_id
}
.dedup.txt | python3
${
ext_chr_srt
}
${
LAMBDA_chr_name
}
|
pi
gz -nc >
\
${
output_dir
}
/CHG_context_LAMBDA_
${
sample_id
}
.dedup.txt.gz &&
\
gz
ip
${
output_dir
}
/CHG_context_b38_
${
sample_id
}
.dedup.txt
\n\
pi
gz
${
output_dir
}
/CHG_context_b38_
${
sample_id
}
.dedup.txt
\n\
zcat
${
output_dir
}
/C??_context_LAMBDA_
${
sample_id
}
.dedup.txt.gz | python3
${
meth_prop
}
>
\
${
output_dir
}
/Lambda_meth_prop_
${
sample_id
}
.txt"
echo
-e
${
CMD
}
>
${
log_dir
}
/1_5_extract_MT.cmd
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment