vgteam / sequenceTubeMap

displays multiple genomic sequences in the form of a tube map
MIT License
180 stars 25 forks source link

Add a GAM example that has read group and sample name set #305

Closed adamnovak closed 1 year ago

adamnovak commented 1 year ago

This should help with #294.

To make this, I took one of the existing GAMs, converted to JSON, set the fields from the schema with jq, and converted back.

vg view -aj cactus0_10.sorted.gam | jq '.sample_name="Sample001" | .read_group="SomeRandomReadGroup"' | vg view -JGa - >cactus_metadata.gam
vg gamsort -i cactus_metadata.sorted.gam.gai cactus_metadata.gam > cactus_metadata.sorted.gam