Questions tagged [strandedness]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
9 votes
0 answers

Run cuffcompare in strand-agnostic mode

Is there a way to run Cufflinks' cuffcompare in a strand-agnostic mode? I would like to do this because I have some RNA-seq datasets derived from an unstranded run, that should be compared to a ...
aechchiki's user avatar
  • 2,676
4 votes
0 answers

R package equivalent to RSeQC infer_experiment to get strandedness of RNA-Seq

I am currently writing an R package that includes a module to run featureCounts (gene quantification tool) from Rsubread. I wanted to be able to specify the correct strandedness option to ...
hmgeiger's user avatar
3 votes
1 answer

Discrepancy with featurecounts analysis using a forward stranded and reverse stranded protocol

My RNAseq analysis pipeline is as follows: fastqc (read quality is good, some overrepresentation of adaptor sequence) → trimmomatic (trimmed adaptor sequence, qc report after trimming suggests the ...
xtian's user avatar
  • 31
2 votes
1 answer

Why is there antisense sequence in RNAseq data

I'm looking at RNAseq data from CCLE. The data is paired-end. Take the cell line Hs578T and the gene HRAS as an example. The cell line carries a G12D mutation (c.35G>A), so the change in cds is: <...
geom_na's user avatar
  • 237
1 vote
4 answers

At what processing step should library strandedness type be taken into account?

Suppose I have single-end RNA-seq data for which the reads in the fastq file are reversed with respect to the original extracted RNAs. Suppose I have the following workflow: Map the reads on the ...
bli's user avatar
  • 3,130
1 vote
1 answer

How to detect bedtools version used by pybedtools (in order to correctly preserve strand information when merging gtf records)

I have issues working with pybedtools due to strand information being lost after extracting and merging transcript coordinates from a gtf file. It seems that the solution to preserve strand ...
bli's user avatar
  • 3,130
0 votes
1 answer

HISAT2: RNA strandedness

My library is unstranded and the code I'm trying to use is this: ...
Neophytos Kouphou's user avatar
0 votes
2 answers

Does rRNA depletion protocol give higher number of mapped reads in Intronic regions?

Recently, I have downloaded a publicly available dataset, which are 350 tumor samples. I see the following information from the published paper. They used Ribo Zero Gold and rRNA was depleted. Strand ...
maven's user avatar
  • 1