From: martinahansen <martinahansen@74ccb610-7750-0410-82ae-013aeee3265d>
Date: Tue, 2 Nov 2010 12:43:04 +0000 (+0000)
Subject: added find_mids biopiece
X-Git-Url: https://git.donarmstrong.com/?a=commitdiff_plain;h=5f3987b6ce1038a55abd7e0021c8bab7009e9263;p=biopieces.git

added find_mids biopiece

git-svn-id: http://biopieces.googlecode.com/svn/trunk@1154 74ccb610-7750-0410-82ae-013aeee3265d
---

diff --git a/bp_bin/find_mids b/bp_bin/find_mids
new file mode 100755
index 0000000..de00e84
--- /dev/null
+++ b/bp_bin/find_mids
@@ -0,0 +1,104 @@
+#!/usr/bin/env ruby
+
+# Copyright (C) 2007-2010 Martin A. Hansen.
+
+# This program is free software; you can redistribute it and/or
+# modify it under the terms of the GNU General Public License
+# as published by the Free Software Foundation; either version 2
+# of the License, or (at your option) any later version.
+
+# This program is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+# GNU General Public License for more details.
+
+# You should have received a copy of the GNU General Public License
+# along with this program; if not, write to the Free Software
+# Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
+
+# http://www.gnu.org/copyleft/gpl.html
+
+# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
+
+# This program is part of the Biopieces framework (www.biopieces.org).
+
+# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> DESCRIPTION <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
+
+#
+
+# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
+
+
+require 'biopieces'
+require 'pp'
+
+mids = %w{ ACGAGTGCGT ACGCTCGACA AGACGCACTC AGCACTGTAG ATCAGACACG
+           ATATCGCGAG CGTGTCTCTA CTCGCGTGTC TAGTATCAGC TCTCTATGCG
+           TGATACGTCT TACTGAGCTA CATAGTAGTG CGAGAGATAC ATACGACGTA
+           TCACGTACTA CGTCTAGTAC TCTACGTAGC TGTACTACTC ACGACTACAG
+           CGTAGACTAG TACGAGTATG TACTCTCGTG TAGAGACGAG TCGTCGCTCG
+           ACATACGCGT ACGCGAGTAT ACTACTATGT ACTGTACAGT AGACTATACT
+           AGCGTCGTCT AGTACGCTAT ATAGAGTACT CACGCTACGT CAGTAGACGT
+           CGACGTGACT TACACACACT TACACGTGAT TACAGATCGT TACGCTGTCT
+           TAGTGTAGAT TCGATCACGT TCGCACTAGT TCTAGCGACT TCTATACTAT
+           TGACGTATGT TGTGAGTAGT ACAGTATATA ACGCGATCGA ACTAGCAGTA
+           AGCTCACGTA AGTATACATA AGTCGAGAGA AGTGCTACGA CGATCGTATA
+           CGCAGTACGA CGCGTATACA CGTACAGTCA CGTACTCAGA CTACGCTCTA
+           CTATAGCGTA TACGTCATCA TAGTCGCATA TATATATACA TATGCTAGTA
+           TCACGCGAGA TCGATAGTGA TCGCTGCGTA TCTGACGTCA TGAGTCAGTA
+           TGTAGTGTGA TGTCACACGA TGTCGTCGCA ACACATACGC ACAGTCGTGC
+           ACATGACGAC ACGACAGCTC ACGTCTCATC ACTCATCTAC ACTCGCGCAC
+           AGAGCGTCAC AGCGACTAGC AGTAGTGATC AGTGACACAC AGTGTATGTC
+           ATAGATAGAC ATATAGTCGC ATCTACTGAC CACGTAGATC CACGTGTCGC
+           CATACTCTAC CGACACTATC CGAGACGCGC CGTATGCGAC CGTCGATCTC
+           CTACGACTGC CTAGTCACTC CTCTACGCTC CTGTACATAC TAGACTGCAC
+           TAGCGCGCGC TAGCTCTATC TATAGACATC TATGATACGC TCACTCATAC
+           TCATCGAGTC TCGAGCTCTC TCGCAGACAC TCTGTCTCGC TGAGTGACGC
+           TGATGTGTAC TGCTATAGAC TGCTCGCTAC ACGTGCAGCG ACTCACAGAG
+           AGACTCAGCG AGAGAGTGTG AGCTATCGCG AGTCTGACTG AGTGAGCTCG
+           ATAGCTCTCG ATCACGTGCG ATCGTAGCAG ATCGTCTGTG ATGTACGATG
+           ATGTGTCTAG CACACGATAG CACTCGCACG CAGACGTCTG CAGTACTGCG
+           CGACAGCGAG CGATCTGTCG CGCGTGCTAG CGCTCGAGTG CGTGATGACG
+           CTATGTACAG CTCGATATAG CTCGCACGCG CTGCGTCACG CTGTGCGTCG
+           TAGCATACTG TATACATGTG TATCACTCAG TATCTGATAG TCGTGACATG
+           TCTGATCGAG TGACATCTCG TGAGCTAGAG TGATAGAGCG TGCGTGTGCG
+           TGCTAGTCAG TGTATCACAG TGTGCGCGTG
+}
+
+count_hash = Hash.new { |hash, key| hash[key] = 0 }
+mid_hash   = {}
+
+mids.each_with_index do |mid, i|
+  mid_hash[mid] = true
+end
+
+casts = []
+
+bp = Biopieces.new
+
+options = bp.parse(ARGV, casts)
+
+bp.each_record do |record|
+  if record.has_key? :SEQ
+    if mid_hash.has_key? record[:SEQ][4 ... 14]
+      count_hash[record[:SEQ][4 ... 14]] += 1
+    end
+  end
+  bp.puts record
+end
+
+mids.each_with_index do |mid, i|
+  if count_hash[mid] > 0
+    record = {}
+    record[:REC_TYPE]  = "MID"
+    record[:MID_NUM]   = i
+    record[:MID_COUNT] = count_hash[mid]
+    record[:MID_SEQ]   = mid
+    bp.puts record
+  end
+end
+
+# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
+
+
+__END__