From Genome Analysis Wiki
Jump to navigationJump to search
1 byte added
, 14:47, 9 March 2011
Line 13: |
Line 13: |
| setenv OUT ~/seq/wednesday/output | | setenv OUT ~/seq/wednesday/output |
| mkdir --p ${OUT} | | mkdir --p ${OUT} |
− |
| |
− | or run
| |
− |
| |
− | csh /home/hyun/wed/bin/sh/step0.sh
| |
| | | |
| 1. Align using BWA | | 1. Align using BWA |
Line 29: |
Line 25: |
| | | |
| ${BIN}/bwa sampe ${REF}/human_g1k_v37_chr20.fa ${OUT}//NA12878.exon-targetted.ILLUMINA.chr20.19986kb-20281kb.read1.fastq.gz.sai ${OUT}/NA12878.exon-targetted.ILLUMINA.chr20.19986kb-20281kb.read2.fastq.gz.sai ${IN}/NA12878.exon-targetted.ILLUMINA.chr20.19986kb-20281kb.read1.fastq.gz ${IN}/NA12878.exon-targetted.ILLUMINA.chr20.19986kb-20281kb.read2.fastq.gz | ${BIN}/samtools-hybrid view -uhS - | ${BIN}/samtools-hybrid sort -m 10000000 - ${OUT}/NA12878.exon-targetted.ILLUMINA.chr20.19986kb-20281kb.paired.bwa.sorted | | ${BIN}/bwa sampe ${REF}/human_g1k_v37_chr20.fa ${OUT}//NA12878.exon-targetted.ILLUMINA.chr20.19986kb-20281kb.read1.fastq.gz.sai ${OUT}/NA12878.exon-targetted.ILLUMINA.chr20.19986kb-20281kb.read2.fastq.gz.sai ${IN}/NA12878.exon-targetted.ILLUMINA.chr20.19986kb-20281kb.read1.fastq.gz ${IN}/NA12878.exon-targetted.ILLUMINA.chr20.19986kb-20281kb.read2.fastq.gz | ${BIN}/samtools-hybrid view -uhS - | ${BIN}/samtools-hybrid sort -m 10000000 - ${OUT}/NA12878.exon-targetted.ILLUMINA.chr20.19986kb-20281kb.paired.bwa.sorted |
| + | |
| + | or run |
| + | |
| + | csh /home/hyun/wed/bin/sh/step1.sh |
| + | |
| | | |
| 2. MERGE ALIGNED BAMS INTO A SINGLE BAM | | 2. MERGE ALIGNED BAMS INTO A SINGLE BAM |