STAR RNA - Generate Wiggle Plots from already aligned BAM files
2
1
Entering edit mode
4.4 years ago
sarajeanne08 ▴ 20

Hi All,

Is it possible to re-run STAR RNA aligner to obtain wiggle plots on already aligned and sorted BAM files? Or do I need to re-align again to obtain this type of output? I really do not want to have to re-align - so I am open to other options.

Thank you!

RNA-Seq • 3.6k views
ADD COMMENT
0
Entering edit mode

Maybe this post can help.

ADD REPLY
2
Entering edit mode
4.4 years ago

bamCoverage from deepTools would be the quickest method to do this. You can also have it output bigWig rather than wiggle files, which is generally preferable (you can directly load those in IGV, for example).

ADD COMMENT
0
Entering edit mode

Do you know if Deeptools can handle paired end reads and strandedness? as I know Bedtools can't handle these properly.

Thank you!

ADD REPLY
0
Entering edit mode

It can handle that without a problem (full disclosure, I'm one of the developers).

ADD REPLY
0
Entering edit mode

Hi Devon,

Thank you so much for the Deeptools recommendation! It works so well and is easy.

Thanks again!

ADD REPLY
2
Entering edit mode
4.4 years ago
lcscs12345 ▴ 30

See https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4631051/

STAR --runMode inputAlignmentsFromBAM \
--inputBAMfile Aligned.sortedByCoord.out.bam \
--outWigType bedGraph \
--outWigStrand Stranded
ADD COMMENT
0
Entering edit mode

Hi Icscs,

Thank you for your help. I tried running those commands but ended up with a bug:

$ STAR --runMode inputAlignmentsFromBAM --inputBAMfile Aligned.sorted.out.trimmed.bam.dupesmarked.bam --outWigType bedGraph --outWigStrand Stranded
Feb 07 11:15:48 ..... Reading from BAM, output wiggle
BUG: alignment extends past chromosome in signalFromBAM.cpp
ADD REPLY

Login before adding your answer.

Traffic: 2478 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6