Quantifying The Area Under A K-Mer Frequency Distribution Curve To Estimate Heterozygosity

Hi,

I have produced a k-mer frequency distribution curve using Jellyfish 1.1.11. It has a clear bi-modal distribution suggesting the sample used to produce it is heterozygous. Is it possible to use this software/graph or any other method to quantify the level of heterozygosity?

Thanks,

Jom.

• 2,575 views • link

written 12.9 years ago by jomaco

0 answers

No answers yet.

More posts like this

Blog: The Complete Guide to Metatranscriptome Sequencing
written 1.8 years ago by Novogene •

It is estimated that humans have more bacteria living in and on our bodies than human cells. Some bacteria play key roles in protection from …
k-mer distribution to estimate the heterozygosity of my assembly
written 5.7 years ago by pablo

Hello, I had PacBio CCS reads I assembled using hifiasm/v12 . I got an assembly of 1.9Gb whereas I expected a genome of 1.3Gb, which …
puzzled about "-kmer" options during de novo assembly
written 8.1 years ago by Yingzi Zhang

Hi all, I am puzzled about "-kmer" options during de novo assembly. **First**, I did k-mer frequency analysis. **Reported**: For P(x): Possible peaks including: 100 …
What is chipmix and freemix in verifybamid
written 8.2 years ago by GK1610

have question on FREEMIX parameter I have this from you Hyun Min Kang "The key idea of FREEMIX estimate is to use excessive heterozygosity to …
Running ABySS at k-mer > 97
written 9.0 years ago by Anand Rao

This post is related to, but different (a continuation) from my earlier post at https://www.biostars.org/p/268060/ seeking help with running the genome assembler software - ABySS. …
Estimate the genome heterozygosity from a PacBio dataset
written 9.1 years ago by gabri

Hi! I know that in the case of Illumina sequences it is possible to use jellyfish for the k-mers count, using short sequences and the …
batch effect in radseq
written 9.7 years ago by oselm •

Hi all, We sent 600 samples of a plant species (no reference genome) for DNA-extraction, library prep and ddRadSeq (paired end, Illumina Hiseq). Firstly, 200 …
How would you write script to count all kmers to estimate the genome size, similar to the one in km…
written 10.7 years ago by Tom •

I'm trying to figure out how KmerGenie is estimating which good k-mers to count in each k-mer frequency histogram. Supposedly the script that counts all …
K-Mer Counting From Single-End, Paired-End And Mate-Pair Reads
written 12.8 years ago by Alice

Hello, biostars! I had one small issue, which grew up in a big question. My issue is following: i'm going to count k-mers and it's …
K-Mer Correction In Rna-Seq Data For Transcriptome Assembly?
written 14.7 years ago by Ryan Thompson

In whole-genome high-throughput sequencing data, one expects a clear separation between high-frequency k-mers (signal) and low-frequency k-mers (noise arising from sequencing errors): <img src='http://www.homolog.us/blogs/wp-content/uploads/2011/09/Capture.png' alt='k-mer …