Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

odgi paths: private/core/accessory sequence classification #565

Merged
merged 5 commits into from
Mar 28, 2024

Conversation

AndreaGuarracino
Copy link
Member

@AndreaGuarracino AndreaGuarracino commented Mar 14, 2024

Node classification by coverage:

odgi paths -i test/chr6.C4.gfa --coverage-levels 2,50,90 | sort -k 1,1n | head | column -t
  #node.id  node.len  class
  1         816       c>=90
  2         1         c<2
  3         1         50<=c<90
  4         125       c>=90
  5         1         c<2
  6         1         50<=c<90
  7         122       c>=90
  8         1         c<2
  9         1         50<=c<90
 
# classes
odgi paths -i test/chr6.C4.gfa --coverage-levels 2,50,90 | grep '^#' -v | cut -f 3 | sort -r | uniq -c
  795 c>=90
  478 50<=c<90
  320 2<=c<50
  155 c<2

Path range classification by fraction:

odgi paths -i test/chr6.C4.gfa --fraction-levels 0.1,0.5,1.0 --path-range-class | head | column -t   
  #path.name                    start  end   class
  chm13#chr6:31825251-31908851  0      816   f>=1
  chm13#chr6:31825251-31908851  816    817   0.5<=f<1
  chm13#chr6:31825251-31908851  817    942   f>=1
  chm13#chr6:31825251-31908851  942    943   0.5<=f<1
  chm13#chr6:31825251-31908851  943    1065  f>=1
  chm13#chr6:31825251-31908851  1065   1066  0.5<=f<1
  chm13#chr6:31825251-31908851  1066   1547  f>=1
  chm13#chr6:31825251-31908851  1547   1548  0.5<=f<1
  chm13#chr6:31825251-31908851  1548   1566  f>=1
  
odgi paths -i test/chr6.C4.gfa --fraction-levels 0.1,0.5,1.0 --path-range-class --show-step-ranges | head | column -t   
  #path.name                    start  end   class     steps
  chm13#chr6:31825251-31908851  0      816   f>=1       1+
  chm13#chr6:31825251-31908851  816    817   0.5<=f<1   3+
  chm13#chr6:31825251-31908851  817    942   f>=1       4+
  chm13#chr6:31825251-31908851  942    943   0.5<=f<1   6+
  chm13#chr6:31825251-31908851  943    1065  f>=1       7+
  chm13#chr6:31825251-31908851  1065   1066  0.5<=f<1   9+
  chm13#chr6:31825251-31908851  1066   1547  f>=1       10+
  chm13#chr6:31825251-31908851  1547   1548  0.5<=f<1   12+
  chm13#chr6:31825251-31908851  1548   1566  f>=1       13+

@AndreaGuarracino AndreaGuarracino merged commit aa08eae into master Mar 28, 2024
1 check passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

1 participant