File: blast2gff.rb

package info (click to toggle)
genometools 1.6.1%2Bds-3
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 50,412 kB
  • sloc: ansic: 271,241; ruby: 30,339; python: 4,880; sh: 3,193; makefile: 1,194; perl: 219; pascal: 159; haskell: 37; sed: 5
file content (7 lines) | stat: -rwxr-xr-x 253 bytes parent folder | download | duplicates (8) 
1
2
3
4
5
6
7
 
#!/usr/bin/env ruby
# transform blast tabular sseqid sstart send sstrand to gff3

ARGF.each do |line|
  m = line.match /\S+\s+(\S+)\s+\S+\s+\d+\s+\d+\s+\d+\s+(\d+)\s+(\d+)/
  puts "#{m[1]}\tBlast\tmatch\t#{m[2]}\t#{m[3]}\t.\t?\t.\tName=Regular Hit"
end