16S-pipeline/0.setup at master · BeckResearchLab/16S-pipeline · GitHub

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
#!/bin/bash

source globals.sh

clean=1
if [ "$#" -ge 1 ]
then
	if [ "$1" == "clean" ]
	then
		echo "forcing regeneration of all files"
		clean=1
	else
		echo "unknown argument: $1"
		exit 1
	fi
fi

mkdir -p data
cd data

echo downloading gold.fa for chimera detection
if [ ! -e gold.fa ]
then
	wget http://drive5.com/uchime/gold.fa
fi

echo creating symlinks to data
ln -sf ../data_raw/fasta-qual-mapping-files/*-full.* .
ln -sf ../data_raw/fasta-qual-mapping-files/*-mapping2.txt* .

echo making fastq
for run in $RUNS
do
	if [ ! -e $run-full.fastq -o "$clean" -eq 1 ]
	then
		echo making fastq for $run
		echo python $FQ2FASTQ $run-full.fasta $run-full.qual
		python $FQ2FASTQ $run-full.fasta $run-full.qual > $run-full.fastq
	fi
done

FQ=""
echo strip barcode and primer and relabel
for run in $RUNS
do
	if [ ! -e $run-clean.fastq -o "$clean" -eq 1 ]
	then
		awk '{ if (substr($0, 1, 1) != "#") { id = $1; printf(">%s\n%s\n", id, $2); } }' ../sample_info.xls > $run-barcodes.fasta
		echo processing $run
		echo python $FQRELABEL $run-full.fastq $PRIMER $run-barcodes.fasta ${run}_
		python $FQRELABEL $run-full.fastq $PRIMER $run-barcodes.fasta ${run}_ > $run-clean.fastq
	fi
	FQ="$FQ $run-clean.fastq"
done

echo combining $FQ into combined.fastq
if [ ! -e combined.fastq -o "$clean" -eq 1 ]
then
	cat $FQ > combined.fastq
fi