From: Magnus M. <mm...@sa...> - 2009-10-28 09:38:05
|
Hi all, as part of LookSeq [1], I have written a BAM alignment image renderer in C++. It can render coverage plots or paired alignments based on apparent fragment size as PNG [2], and can be used as a stand-alone software. Rendering speed is about 1 million read pairs/sec (using Sanger FTP as data source), which is not far from the pure reading speed. Data sources, display mode, and options are passed as command line parameters (run without parameters to see help text). Source is at [3]. Needs to be compiled against samtools and libpng. Cheers, Magnus [1] http://lookseq.sf.net/ [2] http://www.sanger.ac.uk/Teams/Team112/lookseq/lookseq_cpp_renderer_examp le.png [3] http://lookseq.svn.sourceforge.net/viewvc/lookseq/bam_c/render_image.cpp ?view=log -- The Wellcome Trust Sanger Institute is operated by Genome Research Limited, a charity registered in England with number 1021457 and a company registered in England with number 2742969, whose registered office is 215 Euston Road, London, NW1 2BE. |