How do I create BSgenome Bos_taurus_UMD3.1?
1
0
Entering edit mode
mbk0asis • 0
@mbk0asis-7155
Last seen 9 months ago
Korea, Republic Of

Hi.

I have been trying to create BSgenome for Bos_taurus_UMD_3.1.1 over a week but haven't succeeded.

I keep the error message below.

> library(BSgenome)
> forgeBSgenomeDataPkg("/PATH/TO/genome.fa")
Error: Line starting '>1 ...' is malformed!

 

I am really frustrated now because I'm not fluent in both English and R language.

Can somebody explain what is the problem and/or how to create a custom BSgenome step by step?

 

 

bsgenome forge bovine • 989 views
ADD COMMENT
2
Entering edit mode
@herve-pages-1542
Last seen 19 hours ago
Seattle, WA, United States

Of course you cannot call forgeBSgenomeDataPkg() directly on a 2bit or FASTA file. You need to prepare a seed file first and then call forgeBSgenomeDataPkg() on the seed file. All this is explained in the BSgenomeForge vignette. Please see the thread you started on the bioc-devel mailing list for the follow-up on this:

https://stat.ethz.ch/pipermail/bioc-devel/2015-March/007223.html

Thanks,

H.

ADD COMMENT

Login before adding your answer.

Traffic: 364 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6