Skip to content
GitLab
Explore
Sign in
Primary navigation
Search or go to…
Project
B
Blockkurs MicrobialBioinfo
Manage
Activity
Members
Labels
Plan
Issues
Issue boards
Milestones
Wiki
Code
Merge requests
Repository
Branches
Commits
Tags
Repository graph
Compare revisions
Snippets
Build
Pipelines
Jobs
Pipeline schedules
Artifacts
Deploy
Releases
Package Registry
Container Registry
Model registry
Operate
Environments
Terraform modules
Monitor
Incidents
Service Desk
Analyze
Value stream analytics
Contributor analytics
CI/CD analytics
Repository analytics
Issue analytics
Model experiments
Help
Help
Support
GitLab documentation
Compare GitLab plans
Community forum
Contribute to GitLab
Provide feedback
Keyboard shortcuts
?
Snippets
Groups
Projects
Show more breadcrumbs
AMR
Blockkurs MicrobialBioinfo
Merge requests
!53
added script comments
Code
Review changes
Check out branch
Download
Patches
Plain diff
Merged
added script comments
master-patch-944e
into
master
Overview
0
Commits
1
Pipelines
0
Changes
1
Merged
Vanni Benvenga
requested to merge
master-patch-944e
into
master
2 years ago
Overview
0
Commits
1
Pipelines
0
Changes
1
Expand
0
0
Merge request reports
Compare
master
master (base)
and
latest version
latest version
bdbdd9c0
1 commit,
2 years ago
1 file
+
9
−
1
Inline
Compare changes
Side-by-side
Inline
Show whitespace changes
Show one file at a time
exercises/Pangenome_analysis/panaroo.sh
+
9
−
1
Options
@@ -6,24 +6,32 @@
#SBATCH --output=panaroo_%j.out
#SBATCH --error=panaroo_%j.err
# Load singularity module
module load singularityce
# Declare variables
# input data folder
datadir
=
# results folder
wd
=
# Path to singularity images folder
singpath
=
/shares/amr.imm.uzh/bioinfo/singularity/
# Use the input data folder to generate an array containing the sample names
array
=(
$(
ls
$datadir
|
cut
-f
1
-d
'_'
|
uniq
)
)
# Move to the results folder, create input and output folders
cd
$wd
mkdir
-p
panaroo/annotation_files
mkdir
-p
panaroo/panaroo_output
# g
ff file soft link
s
# g
enerate soft links to the .gff files while looping through the sample name
s
for
sample_id
in
"
${
array
[@]
}
"
do
ln
-s
$wd
/
$sample_id
/prokka_annotation/
$sample_id
.gff
$wd
/panaroo/annotation_files/
done
# Runs singularity on the annotations and send the output to the output folder
echo
Starting panaroo
$(
date
+
"%T"
)
cd
panaroo
$singpath
/panaroo_1.3.0--pyhdfd78af_0.sif panaroo
-i
$wd
/panaroo/annotation_files/
*
.gff
-o
./panaroo_output/
--clean-mode
strict
--remove-invalid-genes
Loading