ShujiaHuang

"""
Test pysam
"""
import pysam
import numpy as np


def Test1 () :
    samfile = pysam.Samfile('1006-01.bam', 'rb')

ShujiaHuang

// C++ read fa-sequence
void ReadFaSeq(const char *file, map<string, string> &fa) {

    ifstream I(file);
    if (!I){
        cerr << "Cannot open file : " << file << endl;
        exit(1);
    }

    string tmp, refId;

ShujiaHuang

# Author : Shujia Huang
# Date   : 2013-01-10

#!/usr/bin/perl
use strict;
use warnings;
use Getopt::Long;

my ( $posFile, $regionFile ); 

ShujiaHuang

#!/usr/bin/perl
use strict;
use warnings;

#1. 纯区域的Merge
sub Merge {

    my ( $dis_delta, @region) = @_; # An array , format: [chrId, start, end]

    my ( %prePos, @data, $id, $start, $end );