Whole Genome Haplotyping
0
0
Entering edit mode
4.6 years ago

Hello,

I'm currently looking for recommendations on a method for haplotyping. The input is a .sam file with reads mapped to a (known) reference. The reads come from a species with 2 haplotypes (human), rather than reconstructing the genome I want to reconstruct its both haplotypes.

The reads are long (~10,000 bp) with high error rate and lots of indels (3-15%, nucleotides are errors, a majority are indels).

I've found one in https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4073643/ but do not seem to be publicly available for download. Primarily I have problems to find what methods right now are considered the best, as well as being easy to install and use.

Assembly sequencing genome Haplotyping • 958 views
ADD COMMENT
0
Entering edit mode

Why don't you explicitly write which technology was used? We can guess on nanopore data, but there is no need to be obscure about it.

Have a look at WhatsHap, I haven't used it though.

ADD REPLY
0
Entering edit mode

Sorry, Pacific Biosciences Real-time Sequencer.

ADD REPLY
0
Entering edit mode

So Falcon Unzip maybe?

ADD REPLY
0
Entering edit mode

Thank you, will try that one.

ADD REPLY
0
Entering edit mode

Seems like the program is currently under development and very much not trivial to use. Do you know about something for the slightly less tech-savy? (Like Bwa-Mem and samtools, those I can handle)

ADD REPLY
0
Entering edit mode

Actually I want to do the same (in other species), I had work with illumina and PacBio reads before and I think that the best way to do that is sequencing chromosomes separately, at least you have an overall RDC higher than 500x. You can see disPades assembler, it is a good aproximation https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4449708/

ADD REPLY

Login before adding your answer.

Traffic: 2063 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6