Question: Freebayes usage of --haplotype-length and --max-complex-gap
0
gravatar for jingjin2203
20 months ago by
jingjin220340
jingjin220340 wrote:

Hi All,

I am analyzing my ddRADseq data from a diploid organism using Freebayes. In my case, fragments 450-550bp were sequenced with 300bp paired-end reads. My goal is to figure out the two haplotypes of each sample.

I tried to follow the method below (From BioStars provided by Erik). But I did't quite understand what "--max-complex-gap" does. It seems "--max-complex-gap" is limited by the length of the reads, does that mean I can set it as high as 300 in my case? Also I was wondering what "--haplotype-length" does?

Since ddRADseq does not provide the information across whole genome, I also wonder how Freebayes calls a haplotype across whole genome or across different chromosomes?

Hope someone could help me out. I would really appreciate it. Thank you in advance!

freebayes --max-complex-gap 0 --ploidy 2 input.bam | vcffilter -f "QUAL > 20" >high_confidence.vcf bgzip high_confidence.vcf tabix -p vcf high_confidence.vcf.gz freebayes --max-complex-gap 200 --ploidy 2 --haplotype-basis-alleles high_confidence.vcf.gz input.bam >haplotypes.vcf

sequencing snp ddrad haplotype • 624 views
ADD COMMENTlink written 20 months ago by jingjin220340
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1925 users visited in the last hour