tools/dumpSpliceScoreInformation.py

   1 #!/usr/bin/env python
   2 # -*- coding: utf-8 -*-
   3
   4 import array
   5
   6 def dumpSpliceScoreInformation(in_fn):
   7
   8    pos_fn = in_fn + '.pos'
   9    scores_fn = in_fn + '.Conf'
  10
  11    positions   = array.array('I')
  12    positions.fromfile(open(pos_fn),60)
  13
  14    scores      = array.array('f')
  15    scores.fromfile(open(scores_fn),60)
  16
  17    print [p for p in positions if p < 560]
  18    #print scores[:10]
  19
  20
  21 if __name__ == '__main__':
  22    in_fn = '/fml/ag-raetsch/home/fabio/tmp/interval_query_files/acc/contig_1+'
  23    dumpSpliceScoreInformation(in_fn)
  24
  25    in_fn = '/fml/ag-raetsch/home/fabio/tmp/interval_query_files/don/contig_1+'
  26    dumpSpliceScoreInformation(in_fn)
  27
  28    in_fn = '/fml/ag-raetsch/home/fabio/tmp/interval_query_files/acc/contig_1-'
  29    dumpSpliceScoreInformation(in_fn)
  30
  31    in_fn = '/fml/ag-raetsch/home/fabio/tmp/interval_query_files/don/contig_1-'
  32    dumpSpliceScoreInformation(in_fn)