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.)

I think that 

seqlevelsStyle(vr_new)<-"NCBI" 

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

Best wishes,

Andy

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"

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

Thanks for the pointer!