Entering edit mode
6.0 years ago
novicebioinforesearcher
▴
70
we have RNA seq bam aligned using Hisat2, since hisat provides hisat2_extract_splice_sites.py
script which gives co ordinates for splice junctions was wondering if some one could point us to how we can obtain reads mapping to these junctions something similar to Htseq where it provides genes mapping to genes.
Thanks.
Past useful threads:
How To Extract Spliced Rnaseq Reads
Filtering out the spliced reads from bam file. (@Pierre likely has a version that will let you capture those reads).
So is there any way to do this outside of using those splicing tools? we have curated list of splice junctions that have been seen to change in certain conditions, (around 200 junctions), we did match those 200 with hisat2's script output and now just need reads mapping to it will bedtools perform such operation? sorry i really do not understand @Pierre's solutions although i am sure it works primarily due to lack of understanding towards java programming
duplicate of Extracting Intron-Exon Reads from bam files ?
Out put from your program here can be captured to get the splices reads: A: Filtering out the spliced reads from bam file. ?
will this work just did a bit or more research