FASTA format: Difference between revisions

No edit summary
Line 414:
 
# read tmp fasta file and store to hash
ref = tmpfile
id = seq = ""
fasta = {} of String => String
File.each_line(ref) do |line|
if line.starts_with?(">")
fasta[id] = seq.sub(/\s/, "") if id != ""
fasta[ id] = seqline.subsplit(/\s/, )[0].lstrip(">")
seq = ""
end
else
id = line.split(/\s/)[0].lstrip(">")
seq += ""line
end
else
seq += line
end
end
fasta[id] = seq.sub(/\s/, "")
 
# show fasta component
fasta.each do{ |k,v| puts "#{k}: #{v}"}
puts "#{k}: #{v}"
end
 
</lang>
{{out}}