Question about `vg giraffe`
When I use vg giraffe to align reads to a graph and generate a BAM file, is there a way to determine the positions aligned to any genome (non-reference genomes) in the graph?
• 1,534 views
•
link
updated
by
GenoMax
•
written
by
zhengluo •
•
written
by
zhengluo •
1 answer
The easiest way is to align the reads to the graph in GAM or GAF format. Then you select sample sample as the reference using vg gbwt:
vg gbwt --set-reference sample -g output.gbz --gbz-format -Z input.gbz
And then you convert the mappings to BAM using vg surject with the output graph. If the sample is diploid, both haplotypes will be used as references.
If you want to consider only one haplotype (e.g. 1), you can select path names for that haplotype using:
vg paths -L -S sample -x output.gbz | grep "#1#" > paths.txt
Then you can provide the path name list to vg surject using option --into-paths paths.txt.
• 0 views
•
link
written
by
Jouni Sirén
Log in to answer this question.
More posts like this
-
Is vg giraffe appropriate for RADseq sequencing reads?
written by RianeZ •
-
GATK SNP calling from BAM
written by zhengluo •
-
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 •
-
Identifying non-surjectable sites in a minigraph-cactus graph
written by avril •
-
VG : No reference-sense paths available in the graph; falling back to generic paths.
written by Hang •
-
How should I regenotype VCFs using vg?
written by Ricky •
-
Challenges in Variant Calling and Genotyping with Short-Read Data Mapped to a Pangenome Graph: Seek…
written by Kyran •
-
Substantial reduction in read mapping following VG surject
written by peterdfields •
-
vg: Mapping fewer reads to a genome graph than a linear reference
written by John Lloyd •
