Question: counting the reads mapped to specific coordinate in the genome
0
gravatar for alirezamomeni707
2.9 years ago by
alirezamomeni7070 wrote:

I have RNA-seq data and aligned them to the genome, so I also have bam files. I also have the sequence of tRNA genes in the genome and corresponding coordinates. so, the question is that how can I count the number of reads that map to each coordinate separately in the genome. in other word, I want to count the number of reads that map to each tRNA variants. do you guys know how to do that?

rna-seq genome • 1.2k views
ADD COMMENTlink modified 2.9 years ago • written 2.9 years ago by alirezamomeni7070

did you look for an answer ? A simple google returns the answer.

ADD REPLYlink written 2.9 years ago by Pierre Lindenbaum127k

if you mean to use for instance HTseq_count, I have tried to use it but it needs GTF file which has information for all genes. if you get it from gencode or other websites, the information for the tRNA variants are not included. if I get the specific GTF file from the website relate to the the gene names and structure are different from what is mentioned in the GTF from gencode, HTseq count does not work. that is why I asked this question.

ADD REPLYlink written 2.9 years ago by alirezamomeni7070

Then make your own GTF file using the information you have (chromosome names, tRNA coordinates).

ADD REPLYlink modified 2.9 years ago • written 2.9 years ago by genomax80k

You can likely use your tRNA annotations (if GFF format) directly with several tools, otherwise convert to BED format which you should be able to use as intervals. Counting the reads can be done with several tools

ADD REPLYlink modified 2.9 years ago • written 2.9 years ago by Jon330
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: 1841 users visited in the last hour