Reconstruct SNP haplotypes from Illumina short reads
I have whole genome sequencing data (Illumina paired-end 150bp) for a diverse set of highly polyploid accessions with various levels of ploidy (2n=6-16x). I would like to know if there is a method to get sequencing depths for combination of alleles that are shared on the same read.
For example, for three SNPs with two alleles A and G within 100bp, a given individual could have have the following SNP genotyping :
- SNP1: depth_A=40, depth_G=60
- SNP2: depth_A=10, depth_G=90
- SNP3: depth_A=50, depth_G=50
I would like to obtain the following haplotype genotyping:
- Haplotype: depth_AAA=5, depth_AAG=2,…,depth_GGG=40
I hope this was clear, thanks for your help.
• 17 views
Read more here: Source link