My collection of code snippets.
Bash date format to name output files.
date +%Y%m%d%-H%M%S
Output stdout and stderr to separate files from a bash script.
exec 1>>${LOG_FILE}
exec 2>>${LOG_FILE}
Other options:
exec 3>&1 4>&2
trap 'exec 2>&4 1>&3' 0 1 2 3
exec 1>log.out 2>&1
Checking if file and directory exists
check_file_exists() {
if [ ! -f "$1" ]; then
echo "[trim.sh] error: file '$1' does not exist." >&2
exit
fi
}
check_dir_exists() {
if [ ! -d "$1" ]; then
echo "[trim.sh] error: directory '$1' does not exist, creating it." >&2
mkdir -p $1
fi
}
Subset a fasta file:
samtools faidx genome.fa chr1 chr2:1:2000 > genome_subset.fa