[2] | 1 | #!/usr/bin/env python |
---|
| 2 | |
---|
| 3 | from __future__ import division |
---|
| 4 | |
---|
| 5 | import sys |
---|
| 6 | from galaxy import eggs |
---|
| 7 | import pkg_resources; pkg_resources.require( "pysam" ) |
---|
| 8 | |
---|
| 9 | from pysam import csamtools |
---|
| 10 | from galaxy.visualization.tracks.summary import * |
---|
| 11 | |
---|
| 12 | def main(): |
---|
| 13 | |
---|
| 14 | input_fname = sys.argv[1] |
---|
| 15 | index_fname = sys.argv[2] |
---|
| 16 | out_fname = sys.argv[3] |
---|
| 17 | |
---|
| 18 | bamfile = csamtools.Samfile( filename=input_fname, mode='rb', index_filename=index_fname ) |
---|
| 19 | |
---|
| 20 | st = SummaryTree(block_size=25, levels=6, draw_cutoff=150, detail_cutoff=30) |
---|
| 21 | for read in bamfile.fetch(): |
---|
| 22 | st.insert_range(bamfile.getrname(read.rname), read.pos, read.pos + read.rlen) |
---|
| 23 | |
---|
| 24 | st.write(out_fname) |
---|
| 25 | |
---|
| 26 | if __name__ == "__main__": |
---|
| 27 | main() |
---|