Download MGRRF Bioinformatics Manual - Microbial Gene Research
Transcript
SINA Sina: 1. Download SINA (v 1.2.11, 4 bit, ubuntu 12.04) from: http://www.arb-silva.de/aligner/sina-download/ 2. Unpack the archive $ tar -xvzf sina-1.2.10.tgz 3. Change into the unpacked directory and check that SINA is working $ cd sina-1.2.10 $ ./sina --version SINA v1.2.10 (svn-20432) 4. Get yourself a suitable reference alignment. For SSU and LSU, the SILVA NR datasets are a good starting point. Alternatively, if you have some other reference alignment you want to use in multi-FASTA format, convert it to ARB format like this: $ ./sina -i reference.fasta -o reference.arb –prealigned 5. Make sure it's valid multi-FASTA though. The first word of each header must be unique for each sequence! 6. Try aligning some sequences: $ ./sina -i mysequences.fasta -o aligned.fasta --ptdb reference.arb The first time you do that, the ARB PT server used by SINA to quickly find reference sequences for alignment will build it's index. This may take a while, but you will only have to repeat this if you change the reference alignment. The PT server will also continue to run. Use "killall arb_pt_server" to stop all your running PT servers. 7. If you used the SILVA NR dataset, you can classify your sequences like this: $ ./sina -i mysequences.fasta -o aligned.arb --ptdb reference.arb \ --search --search-db reference.arb \ --lca-fields tax_slv 8. Check the manual to find out about the rest of the options. $ ./sina manual 9. Have fun, find out something great and cite us when you publish (Elmar Pruesse; Jörg Peplies; Frank Oliver Glöckner (2012). SINA: accurate high throughput multiple sequence alignment of ribosomal RNA genes.Bioinformatics 2012; doi: 10.1093/bioinformatics/bts252)