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

Two Bigwig files per sample are generated representing each strand #1275

Open
drimran87 opened this issue Mar 26, 2024 · 3 comments
Open

Two Bigwig files per sample are generated representing each strand #1275

drimran87 opened this issue Mar 26, 2024 · 3 comments

Comments

@drimran87
Copy link

Description of feature

The current rna-seq pipeline, generates two bigwig files one each for each strand. it would be great if like chipseq pipeline, it generate only one bigwig files, would be helpful when comparing chipseq and rnaseq data on IGV.

@MatthiasZepper
Copy link
Member

Conversely, the strand-specific coverage information is what most people need, so I think there are good reasons to keep it like it is.

Consider using the IGV functionality Overlay Data Tracks to combine your tracks when viewing, using WiggleTools to combine them upfront or creating coverage tracks afterwards from your BAM files. You can also directly view your BAM files in IGV.

@drpatelh
Copy link
Member

Yep, I agree. Most users will want to know what the coverage looks like based on the strandedness of the gene itself. We could have an implementation in the pipeline that generates all 3 bigwig files e.g. 1 for each strand and 1 combined. Be useful to know what command we would run with WiggleTools to do this with the current files created by the pipeline and we could make it an opt-in parameter to run this step.

@drpatelh drpatelh added this to the 3.15.0 milestone May 13, 2024
@drpatelh drpatelh modified the milestones: 3.15.0, 3.16.0 May 29, 2024
@tdsone
Copy link

tdsone commented Jun 18, 2024

@drimran87 in case you don't need the processing steps downstream of the read mapping you can set --outWigStrand Unstranded to collapse the wig files into one.

Reference: https://github.com/alexdobin/STAR/blob/master/doc/STARmanual.pdf

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

No branches or pull requests

4 participants