售后

编辑本页后请点击“保存”。请参阅 syntax 了解维基语法。只有在您能改进该页面的前提下才编辑它。如果您想尝试一些东西，请先到 playground 热身。

====bam文件按照比对到参考基因组正负链拆分====

针对双端reads
  /TJPROJ6/RNA_SH/personal_dir/zhangxin/jiaoben/get_forward_reverse_from_bam_for_paired.sh <in.bam> <prefix>
  
  ####代码
  #!/usr/bin/bash
  
  if [ ! $# -eq 2 ];then
          echo
          echo "Split paired-end reads according to the positive and negative strands based on their alignment to the reference genome"
          echo
          echo "Please note the input format of the file: for example ...."
          echo
          echo "get_forward_reverse_from_bam_for_paired.sh <in.bam> <prefix>"
          echo
          exit 1
  fi
  
  #对双端reads 按正负链拆分
  ###### forward strand
  samtools view -b -f 128 -F 16 $1 > $1.fwd1.bam
  samtools view -b -f 80 $1 > $1.fwd2.bam
  samtools merge -f $2_forward.bam $1.fwd1.bam $1.fwd2.bam
  samtools index $2_forward.bam
  
  ###### reverse strand
  samtools view -b -f 144 $1 > $1.rev1.bam
  samtools view -b -f 64 -F 16 $1 > $1.rev2.bam
  samtools merge -f $2_reverse.bam $1.rev1.bam $1.rev2.bam
  samtools index $2_reverse.bam
  
针对单端数据
  /TJPROJ6/RNA_SH/personal_dir/zhangxin/jiaoben/get_forward_reverse_from_bam_for_single.sh <in.bam> <prefix>
  
  ####代码
  #!/usr/bin/bash
  
  if [ ! $# -eq 2 ];then
          echo
          echo "Split single reads according to the positive and negative strands based on their alignment to the reference genome"
          echo
          echo "Please note the input format of the file: for example ...."
          echo
          echo "get_forward_reverse_from_bam_for_single.sh <in.bam> <prefix>"
          echo
          exit 1
  fi
  ###### forward strand
  samtools view -bF 20 $1 > $2_forward.bam
  samtools index $2_forward.bam
  ###### reverse strand
  samtools view -bf 16 $1 > $2_reverse.bam
  samtools index $2_reverse.bam

编辑摘要

当您选择开始编辑本页，即寓示你同意将你贡献的内容按下列许可协议发布： CC Attribution-Share Alike 4.0 International

售后

用户工具

站点工具

页面工具