How to change the output reference for vg giraffe BAM output

I constructed a pan genome graph with mingraph cactus with two reference genomes specified. I am trying to align reads to the graph and produce a BAM file that is specific to one of the references. How can I specify which reference is used during the surjection to the output BAM?

vg bam giraffe

• 1,265 views • link

updated 2.9 years ago by Ram

• written 2.9 years ago by idiaz026 •

1 answer

You need to use this options in vg giraffe

--ref-paths FILE              ordered list of paths in the graph, one per line or HTSlib .dict, for HTSLib @SQ headers

You can see the paths in your graph with vg paths -Lv <graph>.gbz. When specifying the ref paths, you don't need to worry about suffixes like [32333-4444444].

• 0 views • link

written 2.9 years ago by glenn.hickey

More posts like this

vg gbwt missing sample and haplotype metadata when created using minigraph-cactus with vg construct
written 12 months ago by cassiwatt •

Hello, I've found strange behavior surrounding the vg's gbwt index and wanted to report it and potentially figure out why it occurs. When I create …
Question about `vg giraffe`
written 1.4 years ago by zhengluo •

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 …
abra2 error java.lang.NumberFormatException
written 1.7 years ago by chiyong4783 •

Hello, I have constructed a graph-based pan-genome using MC, and then I used VG to align short-read data to the pan-genome. When I used abra2, …
Identifying non-surjectable sites in a minigraph-cactus graph
written 1.8 years ago by avril •

For a graph built with Minigraph-Cactus, is there any way to determine which portions of a non-reference path wouldn't be able to be surjected to …
Detection of chimeric HIV-Host transcripts
written 2.1 years ago by Mrist01 •

I am trying to detect the presences of chimeric HIV-host transcripts in my bulk RNA-seq datasets (single-end reads). The current protocols I have found online …
What path list should be used to surject GAM to GRCh38/CHM13 reference genome
written 2.1 years ago by Hanc •

I have constructed an MC pangenome graph with `--reference CHM13 GRCh38` and want to align short read to the graph with haplotype sampling. In the …
Substantial reduction in read mapping following VG surject
written 3.1 years ago by peterdfields •

I did a small comparison of total mapping of a small short-read dataset to: 1. a graph constructed using cactus-minigraph composed of eight whole genome …
vg for gene presence-absence analysis
written 4.3 years ago by liorglic

I would like to perform a gene presence-absence analysis in order to create a matrix indicating the presence/absence of genes in a set of genomes …
Is it advisable to include reference genome for pan-genome anlysis when working on specific type st…
written 8.3 years ago by ThulasiS

I constructed and analyzed pan-genome but now I got this doubt: Is it advisable to include reference genome for pan-genome analysis when working on specific …
STAR output processing
written 9.5 years ago by kintany •

Hi all, I'm running STAR to align RNA-seq reads separately to paternal and maternal genomes for subsequent analysis of allele-specific expression with GATK tools. So …