Get introns in bam format from alignment reads result
Hey all,
I have alignt rna-seq data with Star against the reference genome hg38. Now I have the results in bam format.
Does anyone know if there is a program that can that detect the introns from the alignment result. Finally I want the only introns in a bed file that are mapped against the reference genome. Before I start programing I want ask if somebody know about nice tools ?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Which Reference Hg38 or T2T?? Both??
written by Jon •
-
extracting data from single end aligned bam file - chromosome, strand, start and stop from cigar
written by Jalil Sharif •
-
Dowload Structural variant VCF file from 1000genome website
written by Chvatil
-
To determine the number of correct splice junctions from the alignment output from HISAT2
written by Chyguc •
-
How to generate start & end coordinates for each chromosome of human ref genome in bed. format (0-b…
written by Alewa
-
detect large indels
written by nakanomasayuki265 •
-
Trimmomatic -> STAR
written by xd_d
-
choosing a suitable max intron size for STAR (plant alignment)
written by Biogeek
-
Annotation file (gff3) for hg38 to be used with MISO
written by dina.hesham139
-
getting reference genome accession number
written by zizigolu
I usually extract reads mapped to exons from BAM files with featureCounts using a gtf file from ENSEMBL. So my guess is that you can also extract the reads that map to the introns only with featureCounts.
This question is not entirely clear if you are just looking get known
exonsas a BED file or are looking to find new exons that may be represented in your file file. If it is the former then these directions are still valid. Otherwise you are going to need to look for split-read mappings in your file.What is your actual end goal?