raise SeqError, "Missing seq_name" if self.seq_name.nil?
raise SeqError, "Missing seq" if self.seq.nil?
record = {}
- record['SEQ_NAME'] = self.seq_name
- record['SEQ'] = self.seq
- record['SEQ_LEN'] = self.length
+ record[:SEQ_NAME] = self.seq_name
+ record[:SEQ] = self.seq
+ record[:SEQ_LEN] = self.length
record
end
">#{seq_name}\n#{seq}\n"
end
+ # Method to reverse complement sequence.
+ def reverse_complement
+ self.reverse
+ self.complement
+ end
+
+ alias revcomp reverse_complement
+
+ # Method to reverse the sequence.
+ def reverse
+ self.seq.reverse!
+ end
+
# Method that complements sequence including ambiguity codes.
def complement
raise SeqError, "Cannot complement 0 length sequence" if self.length == 0