项目作者: tomwhite

项目描述 :
Produce a report on a BAM file listing useful information
高级语言: Java
项目地址: git://github.com/tomwhite/bamreport.git
创建时间: 2018-04-05T11:14:55Z
项目社区:https://github.com/tomwhite/bamreport

开源协议:

下载


  1. BAM=/Users/tom/workspace/gatk/src/test/resources/large/CEUTrio.HiSeq.WGS.b37.NA12878.20.21.bam
  2. ~/sw/bedtools2/bin/bedtools genomecov -ibam $BAM -bga > bam.bedgraph

From http://seqanswers.com/forums/showpost.php?p=116995&postcount=11

  1. DICT=/Users/tom/workspace/gatk/src/test/resources/large/human_g1k_v37.20.21.dict
  2. grep '@SQ' $DICT | tr ':' '\t' | awk -v OFS='\t' '{print $3, $5}' > $DICT.genome
  3. ~/sw/bedtools2/bin/bedtools makewindows -g $DICT.genome -w 1000000 \
  4. | ~/sw/bedtools2/bin/bedtools coverage -a - -b $BAM -counts > $BAM.counts.txt
  5. `

bedtools coverage -hist -abam $BAM


  1. mvn compile
  2. mvn exec:java -Dexec.mainClass=com.cloudera.bamreport.BamReport -Dexec.args="/Users/tom/workspace/gatk/src/test/resources/large/CEUTrio.HiSeq.WGS.b37.NA12878.20.21.bam"
  3. # unmapped BAM
  4. mvn exec:java -Dexec.mainClass=com.cloudera.bamreport.BamReport -Dexec.args="/Users/tom/workspace/gatk/src/test/resources/org/broadinstitute/hellbender/metrics/analysis/QualityScoreDistribution/unmapped.bam"

Sparklines

  1. mvn compile
  2. mvn -q exec:java -Dexec.mainClass=com.cloudera.bamreport.BamCoverageSparkline -Dexec.args="/Users/tom/workspace/gatk/src/test/resources/large/CEUTrio.HiSeq.WGS.b37.NA12878.20.21.bam" > sparkline.html
  3. open sparkline.html