Hi!
I am trying to use vg giraffe to map short reads against my pan-genome graph of a genus of Grasshoppers (created from 30 different mitochondrial genomes). The graph is not the most conserved because the divergence between different sequences is quite high (at least compared to humans).
My problem is that I am mapping very few reads with giraffe. I am assuming that my parameter settings are off (I am using the default atm), but I am unsure of how to adjust them to fit my graph. My minimizer file is set with a k = 20 and a w=10.
Do you have any suggestions?
Thank you in advance for your help! -Nicola
•
written
by
Nicola •
1 answer
Interesting topic - can you be more specific ? What is your command line ? What do you mean by "very few reads" ? Is your short read dataset genomic (then you'd expect few reads to map, right? ).
Log in to answer this question.
-
Is vg giraffe appropriate for RADseq sequencing reads?
written by RianeZ •
-
VG Crash with HPRC V2.0 Graph and Indexes
written by Alexandria A. •
-
Removing PCR duplicates in gam file
written by gulin201909 •
-
Question about vg giraffe, vg surject and Visualization tool for mapping result
written by gulin •
-
Does vg giraffe output secondary and supplimentary alignments?
written by lushjia •
-
The number of variations in the pan-genome is reduced compared to the variations in the input VCF f…
written by Wenke •
-
vg: Mapping fewer reads to a genome graph than a linear reference
written by John Lloyd •
-
VG Giraffe Mapping to a Single Connected Component
written by jdru •
-
Combining gcsa indeces for vg
written by Nicola •
-
Using VG to genotype a GAM using a graph constructed from GFA
written by Michael •
