...
| Code Block |
|---|
| title | solution |
|---|
| collapse | true |
|---|
|
#In the 2bRAD library prep the adapters can be ligated to the fragment in either
#fragment in either orientation (note the mirrored sticky ends left by bcgl).
#So the genomic fragment may have been read from either direction
#check for reverse complement of cut site
grep "GCA......TCG" A_cat.fastq | wc -l
grep "GCA......TCG" B_cat.fastq | wc -l
#201574
#202908
#(216159 + 201574)/426295 = 0.98
#(220686 + 202908)/444053 = 0.95 |
Demultiplexing
Separate the pooled fastq files based on the ligation barcodes found at the end of the reads.
| Code Block |
|---|
| title | solution |
|---|
| collapse | true |
|---|
|
#The command to run trim2bRAD_2barcodes_dedup.pl is complicated, so
#another script -- 2bRAD_trim_launch_dedup.pl -- is used to generate the command for us:
2bRAD_trim_launch_dedup.pl cat.fastq > dedupCommands
#look at the commands file
cat dedupCommands
#returns our next commands to execute:
trim2bRAD_2barcodes_dedup.pl input=A_CTGCAG_R1.fastq site=".{12}CGA.{6}TGC.{12}|.{12}GCA.{6}TCG.{12}" adaptor="AGATC" sampleID=100
trim2bRAD_2barcodes_dedup.pl input=B_GAAGTT_R1.fastq site=".{12}CGA.{6}TGC.{12}|.{12}GCA.{6}TCG.{12}" adaptor="AGATC" sampleID=100 |
...