Shujia Huang ShujiaHuang
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| """ | |
| Test pysam | |
| """ | |
| import pysam | |
| import numpy as np | |
| def Test1 () : | |
| samfile = pysam.Samfile('1006-01.bam', 'rb') |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| // 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
/ Overlap.pl
Created
July 22, 2014 21:17
Find the Overlap region between in the two file
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| # Author : Shujia Huang | |
| # Date : 2013-01-10 | |
| #!/usr/bin/perl | |
| use strict; | |
| use warnings; | |
| use Getopt::Long; | |
| my ( $posFile, $regionFile ); |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| #!/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 ); |
NewerOlder