Question: Read Counts In Intronic Region
3
gravatar for Varun Gupta
7.5 years ago by
Varun Gupta1.1k
United States
Varun Gupta1.1k wrote:

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 intron counts • 4.2k views
ADD COMMENTlink modified 7.4 years ago by Malachi Griffith17k • written 7.5 years ago by Varun Gupta1.1k
4
gravatar for Rm
7.5 years ago by
Rm7.9k
Danville, PA
Rm7.9k wrote:

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

ADD COMMENTlink modified 7.5 years ago • written 7.5 years ago by Rm7.9k
1

look in PomBase

ADD REPLYlink written 7.5 years ago by Rm7.9k

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 REPLYlink written 7.5 years ago by Varun Gupta1.1k

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

ADD REPLYlink written 7.5 years ago by Rm7.9k

bed format: http://uswest.ensembl.org/info/website/upload/bed.html or http://genome.ucsc.edu/FAQ/FAQformat.html#format1

ADD REPLYlink written 7.5 years ago by Rm7.9k

Hey RM

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

ADD REPLYlink written 7.5 years ago by Varun Gupta1.1k
2
gravatar for Malachi Griffith
7.5 years ago by
Washington University School of Medicine, St. Louis, USA
Malachi Griffith17k wrote:

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 COMMENTlink modified 7.5 years ago • written 7.5 years ago by Malachi Griffith17k

Thanks malachig

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

ADD REPLYlink written 7.5 years ago by Varun Gupta1.1k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1392 users visited in the last hour