Changes

From Genome Analysis Wiki
Jump to navigationJump to search
36 bytes removed ,  15:16, 14 February 2012
Line 103: Line 103:  
   foreach file (`ls fastq/SAMPLE*.fastq | cut -f 2 -d '/' | cut -f 1 -d '.'`)
 
   foreach file (`ls fastq/SAMPLE*.fastq | cut -f 2 -d '/' | cut -f 1 -d '.'`)
 
   echo $file
 
   echo $file
   bin/bwa aln -q 15 ref/human_g1k_v37_chr20.fa fastq/$file.fastq > bwa.sai/$file.sai
+
   bwa aln -q 15 ref/human_g1k_v37_chr20.fa fastq/$file.fastq > bwa.sai/$file.sai
   bin/bwa samse -r "@RG\tID:ILLUMINA\tSM:$file" ref/human_g1k_v37_chr20.fa bwa.sai/$file.sai fastq/$file.fastq | \
+
   bwa samse -r "@RG\tID:ILLUMINA\tSM:$file" ref/human_g1k_v37_chr20.fa bwa.sai/$file.sai fastq/$file.fastq | \
     bin/samtools view -uhS - | bin/samtools sort -m 2000000000 - bams/$file
+
     samtools view -uhS - | bin/samtools sort -m 2000000000 - bams/$file
   bin/samtools index bams/$file.bam
+
   samtools index bams/$file.bam
 
   end
 
   end
   Line 113: Line 113:  
  for file in `ls fastq/SAMPLE*.fastq | cut -f 2 -d '/' | cut -f 1 -d '.'`;
 
  for file in `ls fastq/SAMPLE*.fastq | cut -f 2 -d '/' | cut -f 1 -d '.'`;
 
   do echo $file;
 
   do echo $file;
   bin/bwa aln -q 15 ref/human_g1k_v37_chr20.fa fastq/$file.fastq > bwa.sai/$file.sai;
+
   bwa aln -q 15 ref/human_g1k_v37_chr20.fa fastq/$file.fastq > bwa.sai/$file.sai;
   bin/bwa samse -r "@RG\tID:ILLUMINA\tSM:$file" ref/human_g1k_v37_chr20.fa bwa.sai/$file.sai fastq/$file.fastq | \
+
   bwa samse -r "@RG\tID:ILLUMINA\tSM:$file" ref/human_g1k_v37_chr20.fa bwa.sai/$file.sai fastq/$file.fastq | \
     bin/samtools view -uhS - | bin/samtools sort -m 2000000000 - bams/$file;
+
     samtools view -uhS - | samtools sort -m 2000000000 - bams/$file;
   bin/samtools index bams/$file.bam;
+
   samtools index bams/$file.bam;
 
  done
 
  done
  
533

edits

Navigation menu