Guide
Examples
= Prerequisites = * Linux 64-bit environment [http://www.centos.org CentOS] * Python 2.X [http://www.python.org Python]<br> - Under centos: yum install python * MUMmer 3.22 or higher [http://mummer.sourceforge.net MUMmer] * Blast 2.2.25+ or higher [http://blast.ncbi.nlm.nih.gov Blast+] = Installation = Please unpack the tar file.<br> [[Media:CISA_20120425.tar|CISA - 2012.04.25]] Available commands: >tar xvf CISA.tar Execute access: >chmod 755 -R CICA1.0 == Commands and configuration == === Merge.py === Mergy.py can convert data format to fit CISA. '''This is an essentail pre-work'''. Available commands: >python Merge.py Merge.config The content of configuration file: count=3 data=assembly1.fasta,title=Contig_m1 data=assembly2.fasta,title=Contig_m2 data=assembly3.fasta,title=Contig_m3 Master_file=Contigs_m.fa min_length=100 (default:100) Gap=11 #count<br>The number of assemblies you would like to merge #min_length<br>Contigs of length smaller than 100 bp will be discarded. #Gap<br>It's an optional variable.<br>With this variable, we will split the assemblies into contigs at >10 Ns. Please note that '''dos2unix''' may be necessary if your data are in DOS/MAC format. >dos2unix assembly1.fasta assembly2.fasta assembly3.fasta === CISA.py === Available commands: >python CISA.py CISA.config The content of configuration file: genome=genome size infile=file outfile=file nucmer=your installed path/nucmer R2_Gap=0.95 (default:0.95) CISA=your installed path/CISA1.0 makeblastdb=your installed path/makeblastdb blastn=your installed path/blastn #genome<br/>We suggest to use the longest length which is between attended contigs as genome variable.<br/>The break point of CISA will be set to 1.1 * genome variable. #infile<br/>File name with input. #nucmer<br/>Executive file for nucmer. If nucmer has beed set into the path, nucmer variable can be skipped. #makeblastdb<br/>Executive file for makeblastdb. If makeblastdb has beed set into the path, makeblastdb variable can be skipped. #blastn<br/>Executive file for blastn. If blastn has beed set into the path, blastn variable can be skipped. #CISA<br/>Home directory of CISA. #R2_Gap<br/>Tolerant amount of gap during CISA2 step. == Example == <ul><li><strong>Prepair datasets:</strong> the contigs which generate from different softwares you would like to integrate</li></ul> <blockquote>[[Ecoli|Ecoli]]</blockquote> <ul><li><strong>Convert the format and ID of your datasets:</strong><br/></li></ul> <blockquote><p> The content of the configuration file (Merge.config):</p></blockquote> <table width="457" border="1"> <tr> <td height="152" align="left" bgcolor="#66CCFF"> count=5 data=Abyss_contigs.fa,title=Abyss data=Edena_contigs.fa,title=Edena data=SOAPdenovo_contigs.fa,title=SOAP data=CLC_contigs.fa,title=CLC data=Velvet_contigs.fa,title=Velvet Master_file=Contigs.fa </td> </tr> </table> [[Media:Abyss contigs.fa|Abyss_contigs.fa]], [[Media:Edena contigs.fa|Edena_contigs.fa]], [[Media:SOAPdenovo contigs.fa|SOAPdenovo_contigs.fa]], [[Media:CLC contigs.fa|CLC_contigs.fa]], [[Media:Velvetcontigs.fa|Velvet_contigs.fa]] <ul><li><strong>Command:</strong><br /></li></ul><blockquote><p>>python Merge.py Merge.config</p></blockquote><ul><li><strong>Start to integrate contigs:</strong><br /></li></ul><blockquote><p>The content of the configuration file (CISA.config):</p></blockquote> <table width="457" border="1"> <tr> <td height="152" align="left" bgcolor="#66CCFF"> genome=4626205 infile=Contigs.fa nucmer=path/nucmer R2_Gap=0.95 CISA=path/CISA1.0 makeblastdb=path/makeblastdb blastn=path/blastn </td> </tr> </table> <blockquote><p>4626205 which is the longest whole genome between different result from 5 assemblers is set into genome variable.<br /></p></blockquote> <ul><li><strong>Command:</strong><br /></li></ul><blockquote><p>>python CISA.py CISA.config</p></blockquote>
