Search
Question: Error while running SomaticSignatures
0
gravatar for poojitha.ojamies
15 months 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,  :
  sequence chr1 not found
In addition: Warning message:
In .Seqinfo.mergexy(x, y) :
  The 2 combined objects have no sequence levels in common. (Use
  suppressWarnings() to suppress this warning.)

ADD COMMENTlink modified 15 months ago by Andy Lynch80 • written 15 months ago by poojitha.ojamies0
1

What does

seqlevelsStyle(vr_new) 

report? Presumably something different to

seqlevelsStyle(BSgenome.Hsapiens.1000genomes.hs37d5)

?

ADD REPLYlink written 15 months ago by Andy Lynch80
2
gravatar for Andy Lynch
15 months ago by
Andy Lynch80
United Kingdom
Andy Lynch80 wrote:

I think that 

seqlevelsStyle(vr_new)<-"NCBI" 

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

Best wishes,

Andy

ADD COMMENTlink written 15 months ago by Andy Lynch80
0
gravatar for poojitha.ojamies
15 months 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"

ADD COMMENTlink written 15 months ago by poojitha.ojamies0
0
gravatar for poojitha.ojamies
15 months ago by
poojitha.ojamies0 wrote:

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

Thanks for the pointer!

 

ADD COMMENTlink written 15 months ago by poojitha.ojamies0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 135 users visited in the last hour