How to count number of reads aligned to rRNA and tRNA in bacterial genome

0

Entering edit mode

6.7 years ago

Javad ▴ 150

Hi everybody,

I am trying to find number of reads that are aligned to salmonella rRNA and tRNA genes. I use GenomicAlignment R package for counting reads. It works very well for counting the reads aligned to common genes but when it comes to rRNA and tRNA genes, it can not count them because in the related annotation in NCBI, these genes are not referred as genes. Do anybody have any solution for this or maybe a better tool that can handle this issue?

Best

next-gen • 1.8k views

ADD COMMENT • link updated 6.5 years ago by Biostar 20 • written 6.7 years ago by Javad ▴ 150

2

Entering edit mode

Change/add relevant entries for those genes in annotation file you are using.

On a different note: You do not appear to be validating answers for past questions you have posted. Accepting answers (green check mark against the answer, if answer helped you) or up-voting useful comments is the way you show appreciation/validation for effort people put in on this site.

ADD REPLY • link 6.7 years ago by GenoMax 141k

Login before adding your answer.