This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Using repeat region filter post QC

We have NGS whole exome data for 2000 cases and 2000 controls study. Sequencing and QC is done by third party, in theory we already have a "clean" data.

Planning to do gene based tests (skat, burden, etc, or any other analysis), would you advise additionally remove any (all?) of below repeat regions (files are available at UCSC goldenpath):

  • simpleRepeats.txt.gz
  • rmsk.txt.gz
  • genomicSuperDups.txt.gz

We are settled to use simpleRepeats, but no real rationale why. Is it common practice to remove repeat regions, if yes which ones?

We were also suggested to use "blacklists", which ones to use?

ngs exome filter repeats
• 1,764 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by zx8754 zx8754

1 answer

in the VCF I don't remove anything but I flag the variants in the column FILTER.

I use the following resources to flag the variants:

http://hgdownload.soe.ucsc.edu/goldenPath/hg19/encodeDCC/wgEncodeMapability/ ( mappability, ConsensusExcludable...)

the BED provided in this Heng Li's paper: https://doi.org/10.1093/bioinformatics/btu356

the GATK annotation HomopolymerRun

there is also: https://sites.google.com/site/anshulkundaje/projects/blacklists ( Blacklisted genomic regions for functional genomics analysis )

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this