How to make snp calling with GATK under the same depth come true? · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

How to make snp calling with GATK under the same depth come true?

Recently I have dealt with 10 samples of rice genome resequencing. To control the results for comparison, I think it's better to do SNP calling under the same practical average depth. However, I couldn't get the same depth though I used 50X sequencing to all the samples. Is it possible to come true? How?

snp

• 2,622 views • link

written 10.6 years ago by zhengyunchaosky •

2 answers

Perform the variant calling the usual way and then filter the VCF file for minimum and maximum read depth (and other filter of choice). See here, here and here for some discussions on appropriate filters, you can easily find many more.

• 0 views • link

written 10.6 years ago by h.mon

Many thanks for your respond!! really helpful!!

• 0 views • link

written 10.6 years ago by zhengyunchaosky •

Log in to answer this question.

More posts like this

low alignment rate for RNA-Seq data
written 1.5 years ago by Bushra Bibi

Dear scientists, I have RNA-Seq data of rice samples. I analyzed two rice datasets from two different labs (Fastp, Hisat2, Samtools, FeatureCounts). Alignment rate of …
CNVkit analysis with different average depths of control normals
written 2.8 years ago by luna •

Hello. I would like to use CNVkit to detect SCNAs in tumor samples. I have several matched normal samples that were sequenced in the same …
average depth across samples
written 6.3 years ago by ijlal.hyder2012 •

Hello everyone, I need some help/ideas in calculating average depth per SNP across samples. I have already generated depth files for each sample (around 18000 …
What is the definition of "read depth" vs "coverage"? (again...)
written 9.3 years ago by ariel.balter

There are a number of biostars posts on how to calculate coverage and read depth, and what they mean. I'm still confused. This is how …
Somatic mutations from RNA-Seq data of tumor vs. normal
written 10.7 years ago by gil.hornung

<p>Hi,</p> <p>I have RNA-Seq data from a "control" cell line and the same cell line in which a an oncogenic transformation was induced (3 samples …
Does read-depth refer to a fold?
written 11.8 years ago by mangfu100

<p>Hi.</p> <p>In next generation sequencing field, we can easily see the read-depth concept. And I already understand the concept of it. But what is a …
Do I need to sort the .bam file for SNP calling by samtools?
written 12.1 years ago by juncheng

<p>Hi,</p> <p>I'm running samtools for snp calling with human RNAseq data. Is it necessary to sort the tophat aligned bam file for snp calling?</p> <p>I …
Low Coverage Depth On Exome Sequencing
written 12.3 years ago by NB

<p>Hello,</p> <p>I have samples sequenced at BGI using NimbleGen target enrichment kit which shows 55% of targeted sequence was covered 20x and an average coverage …
Forum: Outlier-Based Cnv Calling Method
written 12.6 years ago by Vikas Bansal

We recently published a paper ['Outlier-Based Identification of Copy Number Variations Using Targeted Resequencing in a Small Cohort of Patients with Tetralogy of Fallot'][1]. I …
How To Build The Index Between Genome And Mrnas?
written 15.1 years ago by Haiping •

<p>Hi. I finished calling SNP from genome resequencing data by using samtools. Now I need to know where the SNP located(5'-UTR, exon, intron, 3'UTR, et …