This is a test version of Biostars. For the public version, visit https://www.biostars.org.
vg call on augmented graph blank

Hello again. I've generated a gam (for many samples) using vg: variation graph tool, version v1.27.1 "Deliceto"

./vg map -t 16 -s $xg -g $gcsa -i -f "$R1" -f "$R2" > "$vg"_CP"$numbers"_2.gam

This command works and everything seems to be writing correctly to the .gam, but when I try to generate a vcf from the augmented graph it's completely blank. Here are the vg files' and augmented files' info (is the increase in size for augmented usual? Is that contributing to this issue, if so how much should I try to run this with?

bash-4.2$ ./vg pack -x sxab/mod_F_KS_P_RNA_scaffold_11_gfa.xg -g mod_F_KS_P_RNA_scaffold_11_gfa_CP41.gam -Q 5 -o mod_F_KS_P_RNA_scaffold_11_gfa_CP41.pack

bash-4.2$ ./vg call sxab/mod_F_KS_P_RNA_scaffold_11_gfa.xg -k mod_F_KS_P_RNA_scaffold_11_gfa_CP41.pack

bash-4.2$ ./vg call sxab/mod_F_KS_P_RNA_scaffold_11_gfa.xg -k mod_F_KS_P_RNA_scaffold_11_gfa_CP41.pack -a >> Reference_F_KS_P_RNA_scaffold_11_gfa.vcf

1.6G Jan 24 16:01 mod_F_KS_P_RNA_scaffold_11_gfa_CP41.gam

14M Jan 19 22:22 sxab/pruned_mod_F_KS_P_RNA_scaffold_11_gfa.vg

39M Jan 20 23:22 sxab/mod_F_KS_P_RNA_scaffold_11_gfa.xg

25M Jan 24 23:37 mod_F_KS_P_RNA_scaffold_11_gfa_CP41.pack

#

2.7G Jan 25 00:17 aug_mod_F_KS_P_RNA_scaffold_11_gfa_CP41.gam

1.1G Jan 25 00:20 aug_mod_F_KS_P_RNA_scaffold_11_gfa_CP41.vg

691M Jan 25 00:29 aug_mod_F_KS_P_RNA_scaffold_11_gfa_CP41.vg.xg

49M Jan 25 00:44 aug_mod_F_KS_P_RNA_scaffold_11_gfa_CP41.pack

Also, is there a way to generate a vcf with invariant regions <.> in the alt column?

vg

I can see that there are no references in the augmented vg/xg I'm generating. Is this an issue with my mapping ./vg map -t 16 -x sxab/*vg -g sxab/pruned_mod_F_KS_P_RNA_scaffold_11_gfa.vg.gcsa -i -f 16_19/CPBWGS_16_TCCGGAGA-GTACTGAC_L003_R1_001.fastq.gz -f 16_19/CPBWGS_16_TCCGGAGA-GTACTGAC_L003_R2_001.fastq.gz > mod_F_KS_P_RNA_scaffold_11_gfa_CP16_2.gam

or my augmenting? $ ./vg augment sxab/*vg mod_F_KS_P_RNA_scaffold_11_gfa_CP16_2.gam -A aug_mod_F_KS_P_RNA_scaffold_11_gfa_CP16_2.gam > aug_pruned_mod_F_KS_P_RNA_scaffold_11_gfa.vg

Hello kingcohn!

We believe that this post does not fit the main topic of this site.

redid the index construction, resolved

For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.

If you disagree please tell us why in a reply below, we'll be happy to talk about it.

Cheers!

1 answer

redid the index construction, resolved

Log in to answer this question.