User:Karmella Haynes/Notebook/PcTF Genomics/2012/05/16

From OpenWetWare
Jump to: navigation, search
Owwnotebook icon.pngPc-TF Genomics Report.pngMain project page
Resultset previous.pngPrevious entry      Next entryResultset next.png


  • Peak-finding: BAM > BED using MACS

Peak-finding: BAM > BED using MACS

  • Done by Mark Teng
  • MACS: Model-based Analysis of ChIP-Seq (Zhang et al. 2008) - Given a sonication size (bandwidth) and a high-confidence fold-enrichment (mfold), MACS slides 2 bandwidth windows across the genome to find regions with tags more than mfold enriched relative to a random tag genome distribution.
  • Arguments and parameters in Excel files from Mark
    • About this run: Attempted to analyze IP of hPc-TF (1_aln.bam) and H3K27me3 (6_aln.bam)
    • 04-03/14 Note: no ChIP input sequencing sample was used; I did the input HiSeq experiment later in 2013
Parameter Sample 1 Sample 6
name NA NA
format AUTO AUTO
ChIP-seq file /scratch/mwteng/haynes/aligned/bam/1_aln.bam /scratch/mwteng/haynes/aligned/bam/6_aln.bam
Control file None None
Effective genome size 2.7e+09 2.7e+09
Band width 300 300
Model fold 10,30 10,30
P value cut-off 1.00e-05 1.00e-05
Range for calc. regional lambda 10,000 bps 10,000 bps
Tag size determined as 100 bps 100 bps
Total tags 8,207,576 106,569,524
Tags after filtering 1,076,608 4,947,906
Maximum duplicate tags 1 1
Redundant rate 0.87 0.95
d 104 233