nickp60 / sis Goto Github PK
View Code? Open in Web Editor NEWa program to generate draft genome sequence scaffolds for prokaryotes
License: GNU General Public License v3.0
a program to generate draft genome sequence scaffolds for prokaryotes
License: GNU General Public License v3.0
## Note: This content was added on Nov. 10 2017 due to issues with the (original site)[http://marte.ic.unicamp.br:8747/]. ## Changes to Implementation - SIS is now available both as standalone scripts (to be install with pip), or as a python package (`scaffoldsis`, not to be confused with a different package called `sis`) to include in other projects - sis.py and multifasta.py return a return code 1 if there is an error in execution - as of 0.1.0, multifasta.py writes scaffolds to stdout as well as to the local directory ------------ | AUTHORS | ------------ - Zanoni Dias <[email protected]> - Ulisses Martins Dias <[email protected]> - Joao Carlos Setubal <[email protected]> Last updated: November 10, 2017 This file explains how to run the program SIS in standalone mode. ------------------- | REQUIREMENTS | ------------------- SIS requires an installed python distribution, version 2.5 or later. We have tested SIS on linux (ubuntu), windows 7, and mac OSX. SIS also requires a coordinate file generated by scripts nucmer or promer, which come with the MUMmer package. MUMmer can be obtained here: http://mummer.sourceforge.net/ ----------- | USAGE | ----------- SIS requires MUMmer coord files as the only input. To run the program, type this in a command line window: python sys.py <coord file> where <coord file> is a file containing coordinate data. This file should be generated by nucmer or promer. As a post-processing step after using nucmer or promer, we recommend the use of the delta-filter program, which is also available in the MUMmer package. As a general example, let "reference.fna" and "contigs.fna" be one genome and a set of contigs respectively. Both files have to be in FASTA format. The following command can generate a coordinate file: <nucmer/promer> reference.fna contigs.fna -p example delta-filter -1 example.delta > example.filter show-coords example.filter > example.coords Then you can use 'example.coords' to run SIS: python sys.py example.coords ------------ | OUTPUT | ------------ SIS generates scaffolds by listing arbitrarily numbered scaffolds followed the names of the contigs sorted by their position in the scaffold order: >scaffold_1 <contig_name_a> <orientation_a> <contig_name_b> <orientation_b <contig_name_c> <orientation_c> ... >scaffold_2 <contig_name_d> <orientation_d> <contig_name_e> <orientation_e <contig_name_f> <orientation_f> ... etc The contig orientation is represented by the numbers 0 or 1, where 0 means the contig should be used in the scaffold in its original orientation and 1 means the contig should be reverse complemented. MULTIFASTA output As part of the SIS package you will find program multifasta.py. This programs accepts 2 inputs: the output of sis.py and the draft genome for which the output was generated: python multifasta.py <sis_file> <draft_genome_in_fasta_format> This program generates a multifasta version of the scaffolds found by sys.py. Example of utilization: <nucmer/promer> reference.fna contigs.fna -p example delta-filter -1 example.delta > example.filter show-coords example.filter > example.coords python sys.py example.coords > example.sis python multifasta.py example.sis contigs.fna multifasta.py creates one multifasta file for each scaffold, using the numbering determined by sis.py to name output files, with .fna extension. ------------- | EXAMPLE | ------------- Some examples can be found in the folder "examples": "nucmer.coords", "promer.coords" and "NC_009525-NC_002945.coords". The first two files are toy examples that are useful for understanding the coord files generated by promer or nucmer. The third file is the result of running SIS on real data. The reference file comes from the organism "Mycobacterium bovis AF2122_97" and the contig set comes from the "Mycobacterium tuberculosis H37Ra". Both files can be acessed at NCBI: reference: ftp://ftp.ncbi.nih.gov/genomes/Bacteria/Mycobacterium_bovis_AF2122_97_uid57695 contigs: ftp://ftp.ncbi.nih.gov/genomes/Bacteria_DRAFT/Mycobacterium_tuberculosis_H37Ra_uid54595 In this example we used 'nucmer' and 'delta-filter' to find 219 conserved blocks and generate the coordinate file. SIS was able to find 204 correct adjacencies from those 219 blocks. ------------------------------ | HOW TO CITE THIS PROGRAM | ------------------------------ If you find SIS useful, please cite the following paper: Z. Dias, U. Dias, J.C. Setubal. SIS: a program to generate draft genome sequence scaffolds for prokaryotes. BMC Bioinformatics, 13(1):96, 2012
A declarative, efficient, and flexible JavaScript library for building user interfaces.
๐ Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. ๐๐๐
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google โค๏ธ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.