VCF files incompatible with signeR
1
0
Entering edit mode
Dan • 0
@dan-19424
Last seen 5.2 years ago

I'm trying to analyze VCF files generated by the Qiagen CLC Genomics Workbench using signeR to detect mutation signatures, and signeR does not like the format of the VCF files for some reason. This is a small panel of targets (12 loci across 4 chromosomes in human samples).

Here is the following code and error message that I get:

library(signeR) library(VariantAnnotation) library(BSgenome.Hsapiens.UCSC.hg19) vcfobj <- readVcf("D:\R\Sample17.NT.vcf", "hg19")

library(Rsamtools) mygenome <- FaFile("D:\R\Homosapienshg19_sequence.fa") mut <- genCountMatrixFromVcf(mygenome, vcfobj) Error in attributes(.Data) <- c(attributes(.Data), attrib) : dims [product 10] do not match the length of object [20]

I'm pretty new with using R, so I'm having trouble interpreting the error message. Does anyone know of any good ways to identify and resolve the issue?

Thanks

software error signeR • 891 views
ADD COMMENT
0
Entering edit mode
Dan • 0
@dan-19424
Last seen 5.2 years ago

I think I may have figured it out; signeR doesn't like VCF files that have variants greater than 1 base. I generated a VCF file that only included single base changes and that worked.

ADD COMMENT
0
Entering edit mode

Hello Dan, thats weird signeR uses VariantAnnotation::isSNV to filter out indels, if you could send me an example of a vcf file with the issue I could look into it better.

ADD REPLY

Login before adding your answer.

Traffic: 791 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