Question: Error while running SomaticSignatures
0
2.7 years ago by
poojitha.ojamies0 wrote:

Hi,

I am trying to run soamtic signature on a test data and enountered the following issue. Can someone explain what this problem is and how to fix it? Thanks!

> sca_motifs = mutationContext(vr_new, BSgenome.Hsapiens.1000genomes.hs37d5)
Error in .getOneSeqFromBSgenomeMultipleSequences(x, name, start, NA, width,  :
In .Seqinfo.mergexy(x, y) :
The 2 combined objects have no sequence levels in common. (Use
suppressWarnings() to suppress this warning.)

somaticsignatures • 884 views
modified 2.7 years ago by Andy Lynch100 • written 2.7 years ago by poojitha.ojamies0
1

What does

seqlevelsStyle(vr_new) 

report? Presumably something different to

seqlevelsStyle(BSgenome.Hsapiens.1000genomes.hs37d5)

?

2
2.7 years ago by
Andy Lynch100
United Kingdom
Andy Lynch100 wrote:

I think that

seqlevelsStyle(vr_new)<-"NCBI"

should allow you to use your first choice of reference genome.

Best wishes,

Andy

0
2.7 years ago by
poojitha.ojamies0 wrote:

Ah, you are right! They are two different builds. Should I use my own FA reference files or is there a better fix for this?

> seqlevelsStyle(vr_new)
[1] "UCSC"
> seqlevelsStyle(BSgenome.Hsapiens.1000genomes.hs37d5)
[1] "NCBI"    "Ensembl"

0
2.7 years ago by
poojitha.ojamies0 wrote:

Fixed! I used BSgenome.Hsapiens.UCSC.hg19 insead and it works :)

Thanks for the pointer!
{"extraction_info": {"found_math": true, "script_math_tex": 0, "script_math_asciimath": 0, "math_annotations": 0, "math_alttext": 0, "mathml": 0, "mathjax_tag": 0, "mathjax_inline_tex": 0, "mathjax_display_tex": 0, "mathjax_asciimath": 1, "img_math": 0, "codecogs_latex": 0, "wp_latex": 0, "mimetex.cgi": 0, "/images/math/codecogs": 0, "mathtex.cgi": 0, "katex": 0, "math-container": 0, "wp-katex-eq": 0, "align": 0, "equation": 0, "x-ck12": 0, "texerror": 0, "math_score": 0.2281656116247177, "perplexity": 13026.590332331507}, "config": {"markdown_headings": false, "markdown_code": true, "boilerplate_config": {"ratio_threshold": 0.18, "absolute_threshold": 10, "end_threshold": 15, "enable": true}, "remove_buttons": true, "remove_image_figures": true, "remove_link_clusters": true, "table_config": {"min_rows": 2, "min_cols": 3, "format": "plain"}, "remove_chinese": true, "remove_edit_buttons": true, "extract_latex": true}, "warc_path": "s3://commoncrawl/crawl-data/CC-MAIN-2019-18/segments/1555578577686.60/warc/CC-MAIN-20190422175312-20190422201312-00110.warc.gz"}