This is a test version of Biostars. For the public version, visit https://www.biostars.org.
New Hg18 Genome In R

Dear All,

Good morning! I would like to ask you a question......

I am doing a Hi-C experiment in a biological material at which at a specific position i have placed a 1.5 KB insertion.

So in the analysis of that experiment ( I am using different R packages such as IRanges) I need to load initially the library(BSgenome.Hsapiens.UCSC.hg18).

But in my experiment the normal hg18 (as you understand) is not the correct one because my material has this 1.5 KB insertion.

So I think that I have to make a new version of the Hg18 to load with this insertion at the specific position

How do you think I can do that? How can I make a new Hg18 library in order to load it?

I will wait for your answer

Best regards Eleni

r genome sequence
• 3,118 views • link
updated by Leonor Palmeira Leonor Palmeira • written by e.karasmani • e.karasmani

I don't really understand your question. If your sequence has an insertion, just add the insertion, and run your analyses on the new sequence.

What type of analyses will you be running? BLAST? We might be able to help you more if you give us more details on what exactly you are trying to run.

• 0 views • link
written by Leonor Palmeira Leonor Palmeira

Yes, please improve this question. What are you doing in R? Using a particular package? And why do you need the whole genome if your insertion is at only one specific location?

• 0 views • link
written by Neilfws Neilfws

1 answer

If you have your own pipelines using the BSgenome 'hg18' package and IRanges, and you do not want to change your scripts, then you might need your custom BSgenome package. Read here how to build your own.

However, depending on the analyses you are running, it might be simpler to adapt your pipelines to run on a modified fasta file containing your insert than to build you own custom BSgenome package.

• 0 views • link
written by Leonor Palmeira Leonor Palmeira

Log in to answer this question.

More posts like this