odgi 把GFA文件按染色體拆分
odgi explode -i joint-clip.P.og -g --prefix=subgraph
$ ls subgraph*
subgraph.0.gfa subgraph.1.gfa subgraph.2.gfa subgraph.3.gfa subgraph.4.gfa
odgi 提取PGGB泛基因組子圖
gfabase 提取MC泛基因組子圖
gfabase sub GRCh38-f1g-90-mc-aug11.gfab GRCh38.chr1:25240000-25460000 --range --connected --view --cutpoints 1 --guess-ranges -o RH_locus.walk.gfa
odgi extract -i chr1.pan.smooth.og -o chr1.pan.RH_locus.og -b chr1.RH_locus.bed -E -P
坐標(biāo)系轉(zhuǎn)換
gfa 有路徑信息和Node信息，有時候為了研究需要將路徑信息和Node信息相互轉(zhuǎn)換
從y路徑到 Node舆声，
odgi position -i k.og -p y,10,+ -v
# 查找 y 路徑上第10個方向為+ 的堿基所在Node的位置
output
#source.path.pos target.graph.pos
y,10,+ 6,0,+
從y路徑到 轉(zhuǎn)x路徑
odgi position -i k.og -p y,10,+ -r x
# 查找 y 路徑上第10個方向為+ 的離x路徑的位置
output
#source.path.pos target.path.pos dist.to.ref strand.vs.ref
y,10,+ x,10,+ 0 +
從Node到路徑
odgi position -i k.og -g 6,2
# 查找 Node 6 偏移量為2在路徑的位置
output
#target.graph.pos target.path.pos dist.to.path strand.vs.ref
6,2,+ x,12,+ 0 +
odgi position -i k.og -g 4 -r x
# 查找 Node 6 偏移量為2在路徑x的位置
給gfa文件添加顏色
gffread -T Sr.gff3 -o Sr.gtf
利用grep 提取一個gene保存為Sr.gff3.adj淮椰，查看
$ cat Sr.gff3.adj
Sr_t1 GETA transcript 5000 7561 . + . SrXXXXXXX0
# 根據(jù)gff文件注釋GFA
odgi position -i merge.og -E gene.gtf.adj.exon > color.csv
odgi view -i merge.og -g merge.gfa
導(dǎo)入到bandage 中
## odgi 不支持walk 信息，而且sort 會重新命名節(jié)點，很討厭主穗。
