rnaseq icon indicating copy to clipboard operation
rnaseq copied to clipboard

Two Bigwig files per sample are generated representing each strand

Open drimran87 opened this issue 10 months ago • 3 comments

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.

drimran87 avatar Mar 26 '24 22:03 drimran87

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.

MatthiasZepper avatar Apr 09 '24 15:04 MatthiasZepper

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 avatar May 13 '24 11:05 drpatelh

@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

tdsone avatar Jun 18 '24 12:06 tdsone