...
Expand |
---|
|
Warning |
---|
title | Submit to the TACC queue |
---|
| Move into your scratch directory and create a new directory: Code Block |
---|
cds
mkdir day2
cd day2
|
Then create a commands file with the commands below and use launcher_creator.py followed by qsub to run them. |
- Mapping
- For BWA, the commands are:
Code Block |
---|
title | ONE LINE command for ALL bwa operations |
---|
|
bwa aln $BI/ngs_course/human_variation/ref/hs37d5.fa $BI/ngs_course/human_variation/allseqs_R1.fastq > r1.sai && bwa aln $BI/ngs_course/human_variation/ref/hs37d5.fa $BI/ngs_course/human_variation/allseqs_R2.fastq > r2.sai && bwa sampe $BI/ngs_course/human_variation/ref/hs37d5.fa r1.sai r2.sai $BI/ngs_course/human_variation/allseqs_R1.fastq $BI/ngs_course/human_variation/allseqs_R1.fastq > test.sam
|
- For bowtie2:
Code Block |
---|
title | One-line command for bowtie |
---|
|
bowtie2 -t -x $BI/ngs_course/human_variation/ref/hs37d5_bowtie2 -1 $BI/ngs_course/human_variation/allseqs_R1.fastq -2 $BI/ngs_course/human_variation/allseqs_R1.fastq -S test.sam
{code:title= |
- Convert, sort, and index output: }# Convert, sort, and index output:
samtools view -S -b
view \-S \-b test.sam > test.bam && samtools sort test.bam test.sorted && samtools index test.sorted.bam | code
|
For a bit of extra challenge, you could make two separate commands
files, one for each mapper, qsub them. Watch out - you'll need to make sure they both write output to different files, or run in different directories!