File: quast_all.sh

package info (click to toggle)
spades 3.13.1+dfsg-2
  • links: PTS, VCS
  • area: main
  • in suites: bullseye, sid
  • size: 22,172 kB
  • sloc: cpp: 136,213; ansic: 48,218; python: 16,809; perl: 4,252; sh: 2,115; java: 890; makefile: 507; pascal: 348; xml: 303
file content (21 lines) | stat: -rwxr-xr-x 1,018 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
#!/bin/bash

############################################################################
# Copyright (c) 2015 Saint Petersburg State University
# Copyright (c) 2011-2014 Saint Petersburg Academic University
# All Rights Reserved
# See file LICENSE for details.
############################################################################

# input_folder references_folder output_folder threads

mkdir $3
mkdir $3/tmp
mkdir $3/all_reports
ls $1 | xargs -P $4 -I {} -t quast --min-contig 1000 --contig-thresholds 5000,8000,12000 -e -R $2/{}.fasta $1/{}/scaffolds.fasta -o $3/tmp/{}
ls $1 | xargs -P $4 -I {} -t cp $3/tmp/{}/report.tsv $3/all_reports/{}.tsv
python3 quast_all.py $3/all_reports $3
grep "( inversion )" $3/tmp/*/contigs_reports/contigs_report_*.stdout > $3/inversions.txt
grep "( translocation )" $3/tmp/*/contigs_reports/contigs_report_*.stdout > $3/translocations.txt
grep "( relocation," $3/tmp/*/contigs_reports/contigs_report_*.stdout > $3/relocations.txt
wc -l $3/*.txt > $3/mis_classification.info