ACGTAGATCGT: "RL9",
ACTACGTCTCT: "RL10",
ACTATACGAGT: "RL11",
- ACTCGCGTCGT: "RL12"
+ ACTCGCGTCGT: "RL12",
+ AGACTCGACGT: "RL13",
+ AGTACGAGAGT: "RL14",
+ AGTACTACTAT: "RL15",
+ AGTAGACGTCT: "RL16",
+ AGTCGTACACT: "RL17",
+ AGTGTAGTAGT: "RL18",
+ ATAGTATACGT: "RL19",
+ CAGTACGTACT: "RL20",
+ CGACGACGCGT: "RL21",
+ CGACGAGTACT: "RL22",
+ CGATACTACGT: "RL23",
+ CGTACGTCGAT: "RL24",
+ CTACTCGTAGT: "RL25",
+ GTACAGTACGT: "RL26",
+ GTCGTACGTAT: "RL27",
+ GTGTACGACGT: "RL28",
+ ACACAGTGAGT: "RL29",
+ ACACTCATACT: "RL30",
+ ACAGACAGCGT: "RL31",
+ ACAGACTATAT: "RL32",
+ ACAGAGACTCT: "RL33",
+ ACAGCTCGTGT: "RL34",
+ ACAGTGTCGAT: "RL35",
+ ACGAGCGCGCT: "RL36",
+ ACGATGAGTGT: "RL37",
+ ACGCGAGAGAT: "RL38",
+ ACGCTCTCTCT: "RL39",
+ ACGTCGCTGAT: "RL40",
+ ACGTCTAGCAT: "RL41",
+ ACTAGTGATAT: "RL42",
+ ACTCACACTGT: "RL43",
+ ACTCACTAGCT: "RL44",
+ ACTCTATATAT: "RL45",
+ ACTGATCTCGT: "RL46",
+ ACTGCTGTACT: "RL47",
+ ACTGTAGCGCT: "RL48",
+ AGACACTCACT: "RL49",
+ AGACATATAGT: "RL50"
}
class BarCodeFinderError < StandardError; end
hamming_dist = 0
barcode = seq[@pos ... @pos + size].upcase.to_sym
- if @barcode_hash.has_key? barcode
+ if @barcode_hash[barcode]
return BarCode.new(barcode, @barcode_hash[barcode], @pos, size, hamming_dist)
elsif @max_mismatches > 0
@barcode_hash.each_key do |key|
Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
input.each_record do |record|
- if record.has_key? :SEQ
+ if record[:SEQ]
if barcode = bc_finder.find_barcode(record[:SEQ])
record.merge!(barcode.to_hash)