1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
|
#!/bin/bash
usage(){
echo "
Written by Brian Bushnell
Last modified September 12, 2018
Description: Sorts reads by name or other keys such as length,
quality, mapping position, flowcell coordinates, or taxonomy.
Intended to merge temp files produced by SortByName if the program
ran out of time during merging.
Usage: mergesorted.sh sort_temp* out=<file>
Input may be fasta, fastq, or sam, compressed or uncompressed.
Parameters:
in=<file,file,...> Input files. Files may be specified without in=.
out=<file> Output file.
delete=t Delete input files after merging.
name=t Sort reads by name.
length=f Sort reads by length.
quality=f Sort reads by quality.
position=f Sort reads by position (for mapped reads).
taxa=f Sort reads by taxonomy (for NCBI naming convention).
sequence=f Sort reads by sequence, alphabetically.
flowcell=f Sort reads by flowcell coordinates.
shuffle=f Shuffle reads randomly (untested).
list=<file> Sort reads according to this list of names.
ascending=t Sort ascending.
memmult=.35 Write a temp file when used memory drops below this
fraction of total memory.
Taxonomy-sorting parameters:
tree= Specify a taxtree file. On Genepool, use 'auto'.
gi= Specify a gitable file. On Genepool, use 'auto'.
accession= Specify one or more comma-delimited NCBI accession to
taxid files. On Genepool, use 'auto'.
Java Parameters:
-Xmx This will set Java's memory usage, overriding autodetection.
-Xmx20g will specify 20 gigs of RAM, and -Xmx200m will specify 200 megs.
The max is typically 85% of physical memory.
-eoom This flag will cause the process to exit if an
out-of-memory exception occurs. Requires Java 8u92+.
-da Disable assertions.
Please contact Brian Bushnell at bbushnell@lbl.gov if you encounter any problems.
"
}
#This block allows symlinked shellscripts to correctly set classpath.
pushd . > /dev/null
DIR="${BASH_SOURCE[0]}"
while [ -h "$DIR" ]; do
cd "$(dirname "$DIR")"
DIR="$(readlink "$(basename "$DIR")")"
done
cd "$(dirname "$DIR")"
DIR="$(pwd)/"
popd > /dev/null
#DIR="$( cd "$( dirname "${BASH_SOURCE[0]}" )" && pwd )/"
CP="$DIR""current/"
z="-Xmx2g"
z2="-Xms2g"
set=0
if [ -z "$1" ] || [[ $1 == -h ]] || [[ $1 == --help ]]; then
usage
exit
fi
calcXmx () {
source "$DIR""/calcmem.sh"
setEnvironment
parseXmx "$@"
if [[ $set == 1 ]]; then
return
fi
freeRam 2000m 84
z="-Xmx${RAM}m"
z2="-Xms${RAM}m"
}
calcXmx "$@"
mergesorted() {
local CMD="java $EA $EOOM $z $z2 -cp $CP sort.MergeSorted $@"
echo $CMD >&2
eval $CMD
}
mergesorted "$@"
|
