Skip to content

GitLab

  • Projects
  • Groups
  • Snippets
  • Help
    • Loading...
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
  • Sign in
A
atacseq_analysis
  • Project overview
    • Project overview
    • Details
    • Activity
    • Releases
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
    • Locked Files
  • Issues 7
    • Issues 7
    • List
    • Boards
    • Labels
    • Service Desk
    • Milestones
    • Iterations
  • Merge Requests 0
    • Merge Requests 0
  • Requirements
    • Requirements
    • List
  • CI / CD
    • CI / CD
    • Pipelines
    • Jobs
    • Schedules
    • Test Cases
  • Security & Compliance
    • Security & Compliance
    • Dependency List
    • License Compliance
  • Operations
    • Operations
    • Incidents
    • Environments
  • Packages & Registries
    • Packages & Registries
    • Container Registry
  • Analytics
    • Analytics
    • CI / CD
    • Code Review
    • Insights
    • Issue
    • Repository
    • Value Stream
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Members
    • Members
  • Collapse sidebar
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
  • BICF
  • Astrocyte
  • atacseq_analysis
  • Issues
  • #29

Closed
Open
Opened Jan 09, 2020 by Holly Ruess@s185797Maintainer

Add FRiP score

total reads

total_reads=$(samtools view -c ${sample}.bam)

reads in peaks

reads_in_peaks=$(bedtools sort -i ${sample}_peaks.narrowPeak
| bedtools merge -i stdin | bedtools intersect -u -nonamecheck
-a ${sample}.bam -b stdin -ubam | samtools view -c)

FRiP score

FRiP=$(awk "BEGIN {print "${reads_in_peaks}"/"${total_reads}"}")

Do this step after calling peaks; add in multiQC report

Colored cutoff: The fraction of reads in called peak regions (FRiP score) should be >0.3, though values greater than 0.2 are acceptable. For EN-TEx tissues (ENCODE GTEx tissue sample), FRiP scores will not be enforced as QC metric. TSS enrichment remains in place as a key signal to noise measure.

Assignee
Assign to
Version 2.0
Milestone
Version 2.0
Assign milestone
Time tracking
None
Due date
None
Reference: BICF/Astrocyte/atacseq_analysis#29