# HG changeset patch # User sanbi-uwc # Date 1483967971 18000 # Node ID 2f3d4c867e1d29ff74cb8ad83d84a002891ba231 # Parent e35fbd921b6bba0f403d9d37563e8f1b4fbe266c planemo upload for repository https://github.com/zipho/novo_sort commit 4db125fcd73327813d7bd511852bbe9ffc68a25c diff -r e35fbd921b6b -r 2f3d4c867e1d novo_sort.py --- a/novo_sort.py Fri Apr 01 10:26:19 2016 -0400 +++ b/novo_sort.py Mon Jan 09 08:19:31 2017 -0500 @@ -7,6 +7,7 @@ import logging log = logging.getLogger( __name__ ) + def novo_sort( bam_filename, output_filename ): cmdline_str = "novosort -c 8 -m 8G -s -f {} -o {}".format( bam_filename, output_filename ) cmdline = newSplit(cmdline_str) @@ -15,6 +16,7 @@ except CalledProcessError: print("Error running the nova-sort", file=sys.stderr) + def newSplit(value): lex = shlex.shlex(value) lex.quotes = '"' @@ -22,12 +24,15 @@ lex.commenters = '' return list(lex) + def main(): parser = argparse.ArgumentParser(description="Re-sorting aligned files by read position") parser.add_argument('output_filename') parser.add_argument('--bam_filename') args = parser.parse_args() - + novo_sort(args.bam_filename, args.output_filename) -if __name__ == "__main__": main() + +if __name__ == "__main__": + main() diff -r e35fbd921b6b -r 2f3d4c867e1d novo_sort.xml --- a/novo_sort.xml Fri Apr 01 10:26:19 2016 -0400 +++ b/novo_sort.xml Mon Jan 09 08:19:31 2017 -0500 @@ -1,29 +1,32 @@ - + Tool to run the Novo Sort mapping tool samtools - novocraft + novoalign - - - - - - novo_sort.py '${out_file}' --bam_filename ${input1} + - + - Help! + + + diff -r e35fbd921b6b -r 2f3d4c867e1d test-data/input1.bam Binary file test-data/input1.bam has changed diff -r e35fbd921b6b -r 2f3d4c867e1d test-data/output1.bam Binary file test-data/output1.bam has changed