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?
1 answer
redid the index construction, resolved
Log in to answer this question.
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.gamor 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.vgHello 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!