Hi Jens, hi Greg,
I would like to provide a gsnap genome index to detect ribosomal RNAs via the 'detectRRNA.rrna' config parameter of your awesome HTSeqGenie package. I couldn't find any details on how to generate a suitable index of ribosomal RNAs in the package documentation, though. (The unit test "test.detectRRNA" uses the TP53Genome() but not real rRNA data.)
For example, do you provide a set of rRNA transcript sequences (e.g. extracted from a TxDb object or a GTF file with gene annotations) and build an index from those sequences? Or is there an even better source for rRNA sequences, given that not all of the rRNA sequences may be annotated in the current reference genome?
Thanks a lot for any pointers and for providing such great tools to the community!
Thomas