I have whole genome sequencing data from different sources. And the sequencing depth varies greatly among different samples.And I joint-called SNPs with all samples.
When I do SNP filtration after SNP calling, I found that the genotype quality is influenced by sequencing depth.
For example, I set GQ >= 30 as my threshold, a genotype that does not meet this threshold is considered a missing genotype. Then after max-missing filtration, I found only very few SNPs are retained, which can not meet the requirements of subsequent analysis.Even if I lowered the threshold to 20, SNPs left are still inadequate.
So I calculated the proportion of variants that meet the threshold in each individual. I found that the higher the sequencing depth of the individual, the higher the proportion of variants that met the threshold .
How can I fix this problem? Can I just skip the GQ filteration? Does that affect the analysis later on, like population structure analysis, demographic analysis or selective sweep indentification?
By the way, the species I study is not a model species, so VQSR in GATK is not usable. And I already did hard-filtring before GQ filtration.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Is it necessary to do genotype quality filtration after snp calling with GATK
Last edited by IdaHao0921; 05-14-2024, 01:27 AM.
Latest Articles
Collapse
-
by seqadmin
Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.
Long-Read Sequencing
Long-read sequencing has seen remarkable advancements,...-
Channel: Articles
12-02-2024, 01:49 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 12-02-2024, 09:29 AM
|
0 responses
159 views
0 likes
|
Last Post
by seqadmin
12-02-2024, 09:29 AM
|
||
Started by seqadmin, 12-02-2024, 09:06 AM
|
0 responses
58 views
0 likes
|
Last Post
by seqadmin
12-02-2024, 09:06 AM
|
||
Started by seqadmin, 12-02-2024, 08:03 AM
|
0 responses
48 views
0 likes
|
Last Post
by seqadmin
12-02-2024, 08:03 AM
|
||
Started by seqadmin, 11-22-2024, 07:36 AM
|
0 responses
77 views
0 likes
|
Last Post
by seqadmin
11-22-2024, 07:36 AM
|