This is a test version of Biostars. For the public version, visit https://www.biostars.org.
calling 1000G SNPs/indels from ngs data (chipseq or other source)

I would like to call known (e.g. from 1000G project) SNPs and indels from sequence reads, in my case a chip-seq project with 20 samples. Is there an easy way to call known SNPs/Indels ? (How) should i take into account the Indels/SNPs during allignment (because an indel is could be quite different from the reference genome)?

chip-seq snp indel sequencing alignment
• 2,747 views • link
written by Tafelplankje Tafelplankje

2 answers

AFAIK you can incorporate such info in the variant calling step with GATK's Variant Recalibration 'known variants' option. Not sure if you can take variants into effect at during assembly.

EDIT: I meant to tell you of the "known variants" option available with GATK, not the Variant Recalibrator. I apologize for the confusion.

• 0 views • link
written by Ram Ram

Maybe this is what you are looking for: http://www.bioconductor.org/packages/release/bioc/html/FunciSNP.html

• 0 views • link
written by Ying W Ying W

Log in to answer this question.

More posts like this