Estimate the genome heterozygosity from a PacBio dataset

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 same length. I would like to know if there is a method / software to estimate the genome heterozygosity from a PacBio dataset. Is jellyfish equally correct?

Thank you!

heterozygosity pacbio genome

• 2,817 views • link

updated 9.1 years ago by Istvan Albert

• written 9.1 years ago by gabri

1 answer

I would say that considering the high error rates of PacBio you can only hope such a method to work if your PacBio run was set up to produce shorter and higher quality consensus sequences. Even then you would need fairly high coverage.

• 0 views • link

written 9.1 years ago by Istvan Albert

Thanks a lot for the answer! This is a very useful advice.

• 0 views • link

written 9.1 years ago by gabri

More posts like this

Questions about using Jellyfish and genomescope with paired end reads?
written 3.2 years ago by asilidae •

I recently used jellyfish with some paired-end reads to count k-mers, and I saw another post that mentioned you should concatenate the paired-end reads before …
jellyfish count discrepancy between results obtained with and without -C option
written 4.6 years ago by ashenflower •

I ran jellyfish on my Illumina reads, first using the -C options, and then without, on the same data. Those are the commands I used: …
How to choose k-mer for the genome size estimation?
written 5.6 years ago by 524730309 •

Hi, there. I am using Jellyfish-2 and Genomescope2 to estimate the genome size of my object. I choose different k-mers (e.g. 17 and 21) and …
How to properly measure the length of a reference genome
written 6.9 years ago by elcortegano

I'm working in a non-model organism for which we do have a reference genome. I was recently learning about the [k-mer sprectrum](https://genomebiology.biomedcentral.com/articles/10.1186/gb-2009-10-10-r108) and learned that …
PacBio long read error correction tools
written 8.4 years ago by freddiejung •

Dear all, in my previous post (https://www.biostars.org/p/303214/#303734), some people suggested that it's better to correct errors in log reads from PacBio sequel using short reads …
Jellyfish-Kmer, Genome Size Estimation
written 12.5 years ago by jackuser1979

I have filtered paired-end genomic DNA sequence of non-model plant, sequenced through Hiseq 2000.Now I want to do k-mer analysis and want to estimate genome …
Ngs-Qc & K-Mer
written 12.5 years ago by williamjohn360

I am new to NGS, I have just received sequenced raw files of non-model plant from Hiseq. I would like to know sequence quality check …
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 …
Quantifying The Area Under A K-Mer Frequency Distribution Curve To Estimate Heterozygosity
written 12.9 years ago by jomaco

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 …
Tool: Jellyfish - Fast, Parallel K-Mer Counting For Dna
written 14.2 years ago by Istvan Albert

JELLYFISH is a tool for fast, memory-efficient counting of k-mers in DNA. A k-mer is a substring of length k, and counting the occurrences of …