Hi, I am working with dUTP bases stranded RNA-Seq paired end data. After aligning with STAR , I would like to convert my bam files to bigwig. I am trying to use bamCoverage but cannot use the --filterRNAstrand in Galaxy. Please can anyone help . Thanks
•
written
by
a.bafna •
1 answer
Hi, thanks for reply. I indeed used "Only include reads originating from fragments from the forward or reverse strand" option and selected reverse . Now I assumed it will only count the plus-strand reads, which originally came from genes on the reverse strand using a dUTP-based method) On viewing the result in UCSC Browser I get the coverage as http://genome.ucsc.edu/cgi-bin/hgTracks?db=mm10&lastVirtModeType=default&lastVirtModeExtraState=&virtModeType=default&virtMode=0&nonVirtPosition=&position=chr1%3A91405146%2D91470160&hgsid=1355475517_CvVtrd8za2fLHooDlDFACr2ZAp72. Is that what I should select and expect. I was expecting the reads to be in downward direction as given after wigtobigwig conversion for reverse reads.
Log in to answer this question.
-
bamCoverage with DNAse-seq shift
written by JMR
-
Strange offsetting in deeptools?
written by Elijah •
-
Can you generate a mappability mask with STAR (or another splice-aware aligner)?
written by mercerf93 •
-
filtering issue when using bamCoverage
written by GA •
-
Bam output of deduplication using UMItools
written by Ati •
-
Normalizing ChIP-seq and ATAC-seq data?
written by star
-
How to run eXpress RNA-seq quantification?
written by k.kathirvel93
-
Apply Deseq2 size factor to aligned reads
written by Bastien Hervé
-
Tophat2 Accepted hits Visualization and HT-seq Count not agreeing
written by Krish Srinivasan •
-
What is the best practice for BAM sorting?
written by Kirill Tsyganov
why are you unable to use
--filterRNAstrand? Do you get an error when trying to use this parameter? If so, what is the error?Hi , In usegalaxy.eu instance there is no option for --filterRNAstrand. Hence I cannot filter them.
Sounds like
deeptoolson that galaxy instance is out of date.All
--filterRNAStrandreally does is selects reads that have mapped to one strand and returns a bigwig file with coverage of those reads. (–filterRNAstrand=forward keeps minus-strand reads, which originally came from genes on the forward strand using a dUTP-based method)You can do this yourself using
samtools viewto select the minus-strand reads from your BAM file and use that as input tobamCoverageto generate your plus-strand coverage bigwig file, then repeat with plus-strand reads.Its not that simple with paired-end reads (because usually, the mates are on different strands), see How To Separate Illumina Based Strand Specific Rna-Seq Alignments By Strand