Skip to content

Instantly share code, notes, and snippets.

@nuada
Last active August 29, 2015 14:03
Show Gist options
  • Select an option

  • Save nuada/e0c813ca5f353656ab48 to your computer and use it in GitHub Desktop.

Select an option

Save nuada/e0c813ca5f353656ab48 to your computer and use it in GitHub Desktop.
Quick and dirty way to assess number of PhiX reads in samples
#!/bin/bash
DATA='/data/project'
for sample in $(cd ${DATA}; ls -1 *.r1.fastq.gz | cut -d . -f 1); do
bowtie2 -x /resources/phix/bowtie2/phix -1 ${DATA}/${sample}.r1.fastq.gz -2 ${DATA}/${sample}.r2.fastq.gz -S /dev/null --al-conc-gz phix-${sample}.r%.fastq.gz -p 4 2>&1 | tee -a phix.log
done
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment