require 'biopieces'
require 'seq'
require 'biopieces'
require 'seq'
def disambiguate(adaptor)
adaptor_disamb = adaptor.dup
adaptor_disamb.gsub!('U', 'T')
def disambiguate(adaptor)
adaptor_disamb = adaptor.dup
adaptor_disamb.gsub!('U', 'T')
pos = self.length - adaptor.length
end
pos = self.length - adaptor.length
end
- #puts self.seq if VERBOSE
- #puts (" " * pos) + adaptor if VERBOSE
+ #puts (" " * pos) + adaptor
ed_max = (adaptor.length * ed_percent * 0.01).round
ed_max = (adaptor.length * ed_percent * 0.01).round
# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<