Bandage
Bandage copied to clipboard
feature request: add Scopes for 'length' and 'connectivity'
It would be useful to be able to filter by contig length, or by connections e.g. to avoid plotting singletons, or combined, to avoid plotting short singletons.
For others looking into this, you can do simple filtering outside of Bandage via perl. These work for FASTG files, but not GFA.
to output sequences >1000 bp/aa from a FASTG file: perl -0076 -ne 'chomp; unless ($.==1) {s/(.*)//; $name=$&; print ">$name$_" if tr/A-Z/A-Z/ > 1000}' assembly_graph.fastg
to print only non-singletons from a FASTG file: perl -0076 -ne 'chomp; unless ($.==1) {print ">$_" if (m/:/);}' assembly_graph.fastg
That is a good idea - I'll add it to my (ever-growing) future feature list. Thanks!