Hello,
I am using vg 1.48.0 to try and calculate SV frequencies from short-read poolseq data mapped to a pangenome using giraffe. I was wondering about how the read depth for the reference vs alternate allele is calculated in vg call - is it the (total bp mapped)/(length of path) or (total bp mapped)/(length of path mapped). I have filtered for mapping quality, so I am wondering if one path had a repetitive region present elsewhere in the genome, this path would get less coverage (and for single sample data, be less likely to be called)? I wasn't sure from Hickey et al. 2020 but apologies if I've missed this somewhere.
Best wishes,
Daniel Wood
1 answer
For SV's it's looking at the average (base by base) read coverage across the entirety of each allele (reference or alt treated the same) in the site. If two alleles share a node, its coverage is split evenly between the two alleles. For smaller variants (<50bp), it uses the minimum instead of average coverage.
Log in to answer this question.
-
Is vg giraffe appropriate for RADseq sequencing reads?
written by RianeZ •
-
vg giraffe successfully mapped?
written by kingcohn •
-
Does vg giraffe output secondary and supplimentary alignments?
written by lushjia •
-
vg call AD frequencies for certain sites
written by Danny •
-
VG giraffe autoindex and alignement mapping warning "Refusing to perform too-large rescue alignment"
written by s_135 •
-
vg view and vg surject both crashing after running giraffe
written by Jen •
-
How to validate the assembly completeness of unaligned reads ?
written by Sony Nguyen •
-
warning[vg::giraffe]: Cannot cluster reads with a fragment distance smaller than read distance
written by cassiwatt •
-
HiFi read alignment with vg giraffe
written by shehzad_99 •
-
vg: Mapping fewer reads to a genome graph than a linear reference
written by John Lloyd •
