Changes

From Genome Analysis Wiki
Jump to navigationJump to search
no edit summary
Line 30: Line 30:  
  ${BIN}/bwa aln -q 15 ${REF}/human_g1k_v37_chr20.fa ${IN}/NA12878.exon.sample.unpaired.fastq.gz ${OUT}/NA12878.exon.sample.unpaired.fastq.gz.sai
 
  ${BIN}/bwa aln -q 15 ${REF}/human_g1k_v37_chr20.fa ${IN}/NA12878.exon.sample.unpaired.fastq.gz ${OUT}/NA12878.exon.sample.unpaired.fastq.gz.sai
   −
  ${BIN}/bwa samse ${REF}/human_g1k_v37_chr20.fa \
+
  ${BIN}/bwa samse ${REF}/human_g1k_v37_chr20.fa ${OUT}/NA12878.exon.sample.unpaired.fastq.gz.sai ${IN}/NA12878.exon.sample.unpaired.fastq.gz | ${BIN}/samtools-hybrid view -uhS - | ${BIN}/samtools-hybrid sort -m 10000000 - ${OUT}/NA12878.exon.sample.unpaired.bwa.sorted
  ${OUT}/NA12878.exon.sample.unpaired.fastq.gz.sai \
  −
  ${IN}/NA12878.exon.sample.unpaired.fastq.gz \
  −
  | ${BIN}/samtools-hybrid view -uhS - \
  −
  | ${BIN}/samtools-hybrid sort -m 10000000 \
  −
  - ${OUT}/NA12878.exon.sample.unpaired.bwa.sorted
   
   
 
   
  ${BIN}/bwa sampe ${REF}/human_g1k_v37_chr20.fa \
+
  ${BIN}/bwa sampe ${REF}/human_g1k_v37_chr20.fa ${OUT}//NA12878.exon.sample.read1.fastq.gz.sai ${OUT}/NA12878.exon.sample.read2.fastq.gz.sai ${IN}/NA12878.exon.sample.read1.fastq.gz ${IN}/NA12878.exon.sample.read2.fastq.gz ${BIN}/samtools-hybrid view -uhS - ${BIN}/samtools-hybrid sort -m 10000000 - ${OUT}/NA12878.exon.sample.paired.bwa.sorted
  ${OUT}//NA12878.exon.sample.read1.fastq.gz.sai \
  −
  ${OUT}/NA12878.exon.sample.read2.fastq.gz.sai \
  −
  ${IN}/NA12878.exon.sample.read1.fastq.gz \
  −
  ${IN}/NA12878.exon.sample.read2.fastq.gz \
  −
  | ${BIN}/samtools-hybrid view -uhS - \
  −
  | ${BIN}/samtools-hybrid sort -m 10000000 \
  −
  - ${OUT}/NA12878.exon.sample.paired.bwa.sorted
      
3. Merge multiple BAMs into one
 
3. Merge multiple BAMs into one

Navigation menu