Hi, I have loaded the vcf of dbSNPs with the readVcf function. Is it possible to get the strand of the SNPs annotated in dbSNPs with variant annotation?

Thank you.

Riccardo

rowRanges(vcf) returns a GRanges with the location and strand of each record. To extract strand only, use strand(rowRanges(vcf)). See ?readVcf for details.

Valerie

Thank you. I am using mutect and all the mutation are in the plus strand. Should I use the option ignore.strand?

Thank you.

Riccardo