Skip to content
GitLab
Explore
Sign in
Primary navigation
Search or go to…
Project
A
atacseq_analysis
Manage
Activity
Members
Labels
Plan
Issues
7
Issue boards
Milestones
Iterations
Wiki
Requirements
Code
Merge requests
0
Repository
Branches
Commits
Tags
Repository graph
Compare revisions
Snippets
Locked files
Build
Pipelines
Jobs
Pipeline schedules
Test cases
Artifacts
Deploy
Releases
Container Registry
Operate
Environments
Monitor
Incidents
Service Desk
Analyze
Value stream analytics
Contributor analytics
CI/CD analytics
Repository analytics
Code review analytics
Issue analytics
Insights
Model experiments
Help
Help
Support
GitLab documentation
Compare GitLab plans
Community forum
Contribute to GitLab
Provide feedback
Keyboard shortcuts
?
Snippets
Groups
Projects
BICF
Astrocyte
atacseq_analysis
Merge requests
!12
Resolve "Fix Pool and Pseudoreps"
Code
Review changes
Check out branch
Download
Patches
Plain diff
Merged
Resolve "Fix Pool and Pseudoreps"
20-FixPoolandPseudo
into
master
Overview
0
Commits
15
Pipelines
1
Changes
11
Merged
Holly Ruess
requested to merge
20-FixPoolandPseudo
into
master
5 years ago
Overview
0
Commits
15
Pipelines
1
Changes
11
Expand
Closes
#20 (closed)
Closes
#19 (closed)
Edited
5 years ago
by
Venkat Malladi
0
0
Merge request reports
Compare
master
master (base)
and
latest version
latest version
586ddd0c
15 commits,
5 years ago
11 files
+
241
−
209
Expand all files
Preferences
File browser
List view
Tree view
Compare changes
Inline
Side-by-side
Show whitespace changes
Show one file at a time
Search (e.g. *.vue) (Ctrl+P)
docs/xcor_header.txt
0 → 100644
+
17
−
0
Options
See https://github.com/crazyhottommy/phantompeakqualtools for more details
COL1: Filename: tagAlign/BAM filename
COL2: numReads: effective sequencing depth i.e. total number of mapped reads in input file
COL3: estFragLen: comma separated strand cross-correlation peak(s) in decreasing order of correlation.
The top 3 local maxima locations that are within 90% of the maximum cross-correlation value are output.
In almost all cases, the top (first) value in the list represents the predominant fragment length.
If you want to keep only the top value simply run
sed -r 's/,[^\t]+//g' <outFile> > <newOutFile>
COL4: corr_estFragLen: comma separated strand cross-correlation value(s) in decreasing order (col2 follows the same order)
COL5: phantomPeak: Read length/phantom peak strand shift
COL6: corr_phantomPeak: Correlation value at phantom peak
COL7: argmin_corr: strand shift at which cross-correlation is lowest
COL8: min_corr: minimum value of cross-correlation
COL9: Normalized strand cross-correlation coefficient (NSC) = COL4 / COL8
COL10: Relative strand cross-correlation coefficient (RSC) = (COL4 - COL8) / (COL6 - COL8)
COL11: QualityTag: Quality tag based on thresholded RSC (codes: -2:veryLow,-1:Low,0:Medium,1:High,2:veryHigh)