Read Counts In Intronic Region
2
3
Entering edit mode
10.5 years ago
Varun Gupta ★ 1.2k

Hi everyone. I am working on some RNA-seq data. I am interested in finding what are the read counts within an intron for my interest of genes. I have a bam file to start with in which 3rd column of the bam file is basically represented as chromosome1 , chromosome2 and chromosome3. So suppose i have some of my genes of interest on chromosome1 , some on chromosome2 and some on chromosome3, how can I find the number of read counts within the introns of those genes.

Again all i have is a BAM file. Any command line code or anything would be helpful

Thanks in advance

read counts intron • 5.6k views
ADD COMMENT
4
Entering edit mode
10.5 years ago
Rm 8.2k

coverageBed -abam input.bam -b introns.bed > introns.bed.coverage

ADD COMMENT
1
Entering edit mode

look in PomBase

ADD REPLY
0
Entering edit mode

Hi RM Thank you for your reply. Since i am new to bedtools, i guess i have to make this introns.bed file i guess. So can you tell me how much columns is necessary in creating an introns.bed file. Is BED3 format ok where we have name in column1 and then start and end coordinates of intronic region. Regards

ADD REPLY
0
Entering edit mode

you can download the intron information in a bed format from ucsc http://genome.ucsc.edu/cgi-bin/hgTables?command=start#Help

ADD REPLY
0
Entering edit mode
ADD REPLY
0
Entering edit mode

Hey RM

Thanks for the help, but since i am working on S.pombe genome it is not their in UCSC

ADD REPLY
2
Entering edit mode
10.5 years ago

Picard has a tool for this that seems to work nicely and gives you this specific metric:

CollectRnaSeqMetrics

Program to collect metrics about the alignment of RNA to various functional classes of loci in the genome: coding, intronic, UTR, intergenic, ribosomal.

Also determines strand-specificity for strand-specific libraries. Full list of: RNA-Seq metric definitions.

To narrow it down to a list of regions/genes of interest, you can filter your BAM and create new BAM with reads mapping to those regions. BEDTools could be used for this step.

ADD COMMENT
0
Entering edit mode

Thanks malachig

I will look into it. BTW can you provide me with a quick code i can go through.. Regards

ADD REPLY

Login before adding your answer.

Traffic: 682 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6