heterozygosity rate

joxcargator73

Member

Join Date: Dec 2012

Posts: 28
- Share
- Tweet
#1

heterozygosity rate

10-08-2014, 03:41 PM

I am trying to have an idea of the level of heterozygosity in a genome of 50 Mb that I just assembled. Assembly with 250 PE MySeq
I tried botiew2 to re-map high quality pair end-reads to my contigs. Them using sam tools created a bam sorted and indexed file. Finally SNPs were called using freebayes. I filtered for Quality >20 and called only for SNPs.
My vcf file contains more than 100,000 SNPs. I believe that some how I am over estimating the number of SNPs. Maybe I am missing something.
Besides a couple of reports http://www.nature.com/ncomms/2013/13...ncomms3673.pdf) I do not have others to compare the level of SNPs I have obtained.
Any comments about this result will be appreciated
Tags: None
Brian Bushnell

Super Moderator

Join Date: Jan 2014

Posts: 2709
- Share
- Tweet
#2

10-08-2014, 04:11 PM

It depends on the organism, but they can have SNP rates of 1/1000 (like human) or less... or as many as 1/100 or more. 100k SNPs in a 50Mbp diploid genome is not surprising.
Comment
joxcargator73

Member

Join Date: Dec 2012

Posts: 28
- Share
- Tweet
#3

10-08-2014, 04:24 PM

Thanks. YES the organism is a oomycete. Usually highly heterozygous. Actually there the filtered VCF file contains 545,000 SNPS
Comment

Previous template Next

Reply to Nine Things a Sample Prep Scientist Thinks About Before Sequencing

by GATTACAT

Love this - good data definitely starts from good input, and poor input can only give relatively poor data. I particularly like the mention of Nanodrop/absorbance based methods for quantification. It's such a toss up if you'll get an accurate reading or what amounts to a randomly generated number, and a lot of library/sequencing related issues can be traced back to poor quant.
- Channel: Articles
Yesterday, 11:43 AM
Nine Things a Sample Prep Scientist Thinks About Before Sequencing

by SEQadmin2

I’m not a sequencing expert. I’m a purification scientist who uses NGS to evaluate workflows my group develops. With this perspective, we think about the sample first and the NGS workflow second. The sequencer is an exceptionally honest reporter, but it can only report on what you give it, so whether you get clean, interpretable data from an NGS workflow is largely determined before you begin.

Here are nine questions we think about, in roughly the order they matter, before...
- Channel: Articles
06-18-2026, 07:11 AM
From Collection to Sequencing: Why Sample Preparation and Preservation Define Sequencing Data

by SEQadmin2

Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.

The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
...
- Channel: Articles
06-02-2026, 10:05 AM

Topics	Statistics	Last Post
New AI Model Captures Long-Range Genomic Signals to Improve RNA Splice Site Prediction by SEQadmin2 Started by SEQadmin2, 06-30-2026, 05:37 AM	0 responses 10 views 0 reactions	Last Post by SEQadmin2 06-30-2026, 05:37 AM
Large-Scale Protein Screen Uncovers Hidden Regulators of Alternative Polyadenylation by SEQadmin2 Started by SEQadmin2, 06-26-2026, 11:10 AM	0 responses 18 views 0 reactions	Last Post by SEQadmin2 06-26-2026, 11:10 AM
Whole-Genome Sequencing Traces Faroe Islands Ancestry to a North Atlantic Founder Population by SEQadmin2 Started by SEQadmin2, 06-17-2026, 06:09 AM	0 responses 52 views 0 reactions	Last Post by SEQadmin2 06-17-2026, 06:09 AM
Sequencing the Two-Toed Sloth Genome Reveals Jumping Genes Tied to Its Extreme Metabolism by SEQadmin2 Started by SEQadmin2, 06-09-2026, 11:58 AM	0 responses 111 views 0 reactions	Last Post by SEQadmin2 06-09-2026, 11:58 AM

Unconfigured Ad

Comment

Comment

Latest Articles

ad_right_rmr

News