Merge pull request #46 from stijndcl/fix/lz4-decompression-leftover

Fix bugs related to LZ4 decompression in non-database steps

scripts/build_database.sh

@@ -732,15 +732,15 @@ create_kmer_index() {
 	log "Started the construction of the $KMER_LENGTH-mer index."
 	for PREFIX in A C D E F G H I K L M N P Q R S T V W Y; do
 		pv -N $PREFIX "$OUTPUT_DIR/uniprot_entries.tsv.lz4" \
-			| gunzip \
+			| $CMD_LZ4CAT \
 			| cut -f4,7 \
 			| grep "^[0-9]*	[ACDEFGHIKLMNPQRSTVWY]*$" \
 			| umgap splitkmers -k"$KMER_LENGTH" \
 			| sed -n "s/^$PREFIX//p" \
 			| LC_ALL=C $CMD_SORT \
 			| sed "s/^/$PREFIX/"
 	done \
-			| umgap joinkmers "$(guz "$OUTPUT_DIR/taxons.tsv.lz4")" \
+			| umgap joinkmers "$(luz "$OUTPUT_DIR/taxons.tsv.lz4")" \
 			| cut -d'	' -f1,2 \
 			| umgap buildindex \
 			> "$OUTPUT_DIR/$KMER_LENGTH-mer.index"
@@ -755,7 +755,7 @@ create_tryptic_index() {
 	have "$TABDIR/sequences.tsv.lz4" || return
 	log "Started the construction of the tryptic index."
 	pv "$TABDIR/sequences.tsv.lz4" \
-		| gunzip \
+		| $CMD_LZ4CAT \
 		| cut -f2,3 \
 		| grep -v "\\N" \
 		| umgap buildindex \

scripts/helper_scripts/.gitignore

@@ -2,6 +2,7 @@
 dat-parser
 functional-analysis
 lcas
+substitute-aas
 taxa-by-chunk
 taxons-lineages
 taxons-uniprots-tables