Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Fragment length parameter not so necessary for bigWig generation #161

Open
cnluzon opened this issue Jan 18, 2023 · 0 comments
Open

Fragment length parameter not so necessary for bigWig generation #161

cnluzon opened this issue Jan 18, 2023 · 0 comments

Comments

@cnluzon
Copy link
Collaborator

cnluzon commented Jan 18, 2023

Since the pipeline runs on paired-end data and according to deepTools documentation on --extendReads parameter:

This parameter allows the extension of reads to fragment size. If set, each read is extended, without exception. NOTE: This feature is generally NOT recommended for spliced-read data, such as RNA-seq, as it would extend reads over skipped regions. Single-end: Requires a user specified value for the final fragment length. Reads that already exceed this fragment length will not be extended. Paired-end: Reads with mates are always extended to match the fragment size defined by the two read mates. Unmated reads, mate reads that map too far apart (>4x fragment length) or even map to different chromosomes are treated like single-end reads. The input of a fragment length value is optional. If no value is specified, it is estimated from the data (mean of the fragment size of all mate reads).

This is already estimated from the BAM alignment as a mean, should be accurate enough and it only affects a small subset of the reads, we could consider removing the parameter from the bigwig rules. However we would still need it for the scaling factors, so we cannot remove the parameter from minute.yaml unless #10 is solved first.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant