@Trivas a trick is to simply paste the link for biostar thread you want to refer to. Biostar expands the title automatically. e.g. How To Extract Spliced Rnaseq Reads
Find exon-exon spanning reads
1 answer
There are a ton of threads with the same/similar question. Here's a good place to start but with a bit of googling you can find several more threads.
I want to generate counts for exon-exon spanning reads only and want to use HTSeq count for this if possible. I'm on the Galaxy website and can provide a BAM file of my sample, but as far as I can see, none of the GENCODE GTF files contains information on whether certain reads are exon-exon spanning reads. So how should I do this?
Log in to answer this question.
More posts like this
-
RNA-seq: identifying pericentromeric region
written by Polly •
-
Which zUMIs output to use for differential gene expression?
written by firestar
-
Determining transcripts for repetitive elements in RNA-seq data
written by rodd
-
Using dropEst output for Seurat
written by annadv •
-
Transcript Specif qPCR primer
written by ravihansa82
-
Identify reads that span intron-exon junctions in RNA-Seq
written by oriolebaltimore
-
Several exonerate options; which to choose?
written by mforthman •
-
Filtering a Tophat BAM file
written by Mehmet
-
Count splice junction reads from SAM/BAM files
written by juncheng
-
How To Understand Some Annotated Terms Made By Annovar
written by zengtony743
Are you worried about DNA contamination?
yes exactly
one hacky way to do it could be to filter only reads that have a skip e.g. a N in the cigar string...grep on the CIGAR column for N
Yes I think this is the solution because spliced reads are the same as exon-exon spanning reads! I hadn't thought of that, thanks!