README.md in bio-blast_seq_extract-0.2.4 vs README.md in bio-blast_seq_extract-0.2.5
- old
+ new
@@ -24,18 +24,22 @@
require 'bio-blast_seq_extract'
a = BlastSeqExtract.new
Read a sequence in the format of fasta
+
a.readSeq(seq_infile)
Read the blast output (tabular format)
+
a.readBlast(blast_infile, evalue_cutoff, max_no_seqs)
Get the aligned region from the blast output
+
seqObjs = a.getRegion
Assumming the input is in FASTA format, f is an object of Bio::FastaFormat
+
seqObjs.each do |f|
puts '>'+f.definition
puts f.seq
end