This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Bulk RNA seq

Hi I just started indexing with Rsubread with R packages 3.6.0. I would like to know a few things in the scripts.

library(Rsubread) 
ref <- system.file("extdata","reference.fa",package="Rsubread") 
buildindex(basename="reference_index",reference=ref)

  1. What does mean about extdata means, in the script? reference.fa--- do I can download this file from https://www.ncbi.nlm.nih.gov/assembly/GCF_000001635.20/ Here in this link, I have two option in the Download Assembly

  2. Source database(Genbank or RefSeq), what I have to choose from here?

  3. File type, which files types, I have to select as my reference?
rna-seq rna-seq
• 1,027 views • link
updated by GenoMax GenoMax • written by B- for bioinformatics • B- for bioinformatics

1 answer

I assume you are following Rsubread vignette.

extdata refers to location where additional data installed with R packages. In this case Rsubread installs a reference file that you can find the location of by just typing ref at command prompt after you execute the second command.

If you are going to work with mouse genome then you could get the reference and annotation from GENCODE/Ensembl/NCBI. You will need a fasta reference sequence and a GTF format annotation file.

• 0 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this