Skip to content

Instantly share code, notes, and snippets.

@MADscientist314
Created April 28, 2021 20:02
Show Gist options
  • Save MADscientist314/a236fccca84ade7c1632c3107f43f1c0 to your computer and use it in GitHub Desktop.
Save MADscientist314/a236fccca84ade7c1632c3107f43f1c0 to your computer and use it in GitHub Desktop.
fastqscreener wrapper
threads=48
for n in $(cat manifest_02);
do
nohup fastq_screen
time
fastq_screen \
--nohits \
--threads $threads \
--subset 0 \
--outdir screened/ \
--conf fastq_screen.tmp.conf \
trimmed/$n"_1_sequence_kneaddata.trimmed.1.fastq.gz
echo "$n complete";
rm -v screened/$n"_1_sequence_kneaddata.trimmed.1.tagged.fastq.gz" > $n.fastqscreen.log;
done
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment