Professional Documents
Culture Documents
UniproUGENE UserManual
UniproUGENE UserManual
Version 1.9.4
Contents
1 About Unipro . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.1 Contacts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
9
9
About UGENE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.1 2.2 2.3 2.4 Key Features . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . User Interface . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . High Performance Computing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Cooperation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 10 11 11
Installation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
3.1 3.2 3.3 Installing UGENE on Windows . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Installing UGENE on Linux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Installing UGENE on Mac OS X . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12 13 14
Basic Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
4.1 4.2 4.3 4.4 UGENE Terminology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Main Menu Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Creating New Project . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Opening Document . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4.4.1 4.4.2 4.5 4.6 4.7 4.8 4.9 4.10 4.11 4.12 Opening for the First Time . . . . . . . . . . . . . . . . . . . . . . . . . . . . Opening Document Present in Project . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 16 17 18 18 19 20 21 22 23 24 24 26 26 27 28 29 29 30 30 31 31
Creating Document
Locked Documents . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Using Objects and Object Views . . . . . . . . . . . . . . . . . . . . . . . . . . . . Using Bookmarks . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Working with Projects . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Adding and Removing Plugins . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Fetching Data from Remote Database . . . . . . . . . . . . . . . . . . . . . . . . . UGENE Application Settings . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4.12.1 4.12.2 4.12.3 4.12.4 4.12.5 4.12.6 4.12.7 4.12.8 General . . . . . . Resources . . . . . Network . . . . . Logging . . . . . . Workow Designer OpenCL . . . . . CUDA . . . . . . External Tools . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.13
32
Sequence View . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
5.1 5.2 5.3 5.4 5.5 5.6 5.7 Sequence View Components . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Global Actions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Sequence Toolbar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Sequence Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Sequence Zoom View . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5.5.1 Managing Zoom View Rows . . . . . . . . . . . . . . . . . . . . . . . . . . . . Sequence Details View . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Manipulating Sequence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5.7.1 5.7.2 5.7.3 5.7.4 5.7.5 5.7.6 5.7.7 5.7.8 5.7.9 5.7.10 5.7.11 5.7.12 5.8 5.8.1 5.8.2 5.9 5.9.1 5.9.2 5.9.3 5.9.4 5.9.5 5.9.6 5.9.7 5.9.8 5.10 Going To Position . . . . . . . . . Toggling Views . . . . . . . . . . . Capturing Screenshot . . . . . . . Zooming Sequence . . . . . . . . . Creating New Ruler . . . . . . . . Selecting Amino Translation . . . . Showing and Hiding Translations . Selecting Sequence . . . . . . . . . Copying Sequence . . . . . . . . . Finding Pattern in Sequence . . . . Editing Sequence . . . . . . . . . . Locking and Synchronize Ranges of . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Several Sequences . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35 37 37 38 39 39 40 41 41 41 41 42 42 42 43 44 46 46 48 50 51 52 52 53 53 54 54 56 57 58 58 58 62
Annotations Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Automatic Annotations Highlighting . . . . . . . . . . . . . . . . . . . . . . . The "db_xref" Qualier . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Creating Annotation . . . . . . . . Editing Annotation . . . . . . . . . Highlighting Annotations . . . . . Creating and Editing Qualier . . . Adding Column for Qualier . . . . Copying Qualier Text . . . . . . . Deleting Annotations and Qualiers Importing Annotations from CSV . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Manipulating Annotations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Alignment Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
6.1 Overview 6.1.1 6.1.2 6.1.3 6.1.4 6.1.5 6.1.6 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Alignment Editor Features . . Alignment Editor Components Navigation . . . . . . . . . . Coloring Schemas . . . . . . Zooming and Fonts . . . . . Searching for Pattern . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65 65 66 67 67 67 68
6.1.7 6.2 6.2.1 6.2.2 6.2.3 6.2.4 6.2.5 6.3 6.3.1 6.3.2 6.3.3 6.4
Consensus . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Undo/Redo Framework Selecting Subalignment Editing Alignment . . . Aligning Sequences . . . Working with Sequences . . . . . . . . . . . . List . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
69 70 70 70 70 72 73 75 75 75 75 76
Assembly Browser . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
7.1 7.2 Import BAM File . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Browsing and Zooming Assembly . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.2.1 7.2.2 7.2.3 7.2.4 7.2.5 7.3 7.4 7.5 Opening Assembler Browser Window . . Assembly Browser Window Components Reads Area Description . . . . . . . . . Assembly Overview Description . . . . . Ruler Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78 79 79 80 80 81 81 82 83 84 84 84 84 85 85 85
Getting Information About Read . . . . . . . . . . . . . . . . . . . . . . . . . . . . Associating Reference Sequence . . . . . . . . . . . . . . . . . . . . . . . . . . . . Exporting . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7.5.1 7.5.2 7.5.3 Exporting Read . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Exporting Visible Reads . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Exporting Image . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Assembly Overview Hotkeys . . . . . . . . . . . . . . . . . . . . . . . . . . . . Reads Area Hotkeys . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
7.6
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
. . . .
90 90 91 91 91 91
Distributed Computing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
9.1 9.2 Remote Machines Monitor 9.2.1 9.2.2 9.2.3 9.2.4 9.2.5 9.3 9.4 9.5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92 94 94 94 94 95 96 97 98 99
Running Workows on Cloud . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Introduction . . . . . . . . . . . . Cloud Computing . . . . . . . . . Cloud Remote Machine . . . . . . Launching Workow . . . . . . . . Useful Tips and Recommendations
Running HMMER3 Search Task on Remote Machine . . . . . . . . . . . . . . . . . Running Smith-Waterman Search Task on Remote Machine . . . . . . . . . . . . . Running MUSCLE Align Task on Remote Machine . . . . . . . . . . . . . . . . . .
10
Plugins . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
10.1 3D Structure Viewer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100 10.1.1 10.1.2 10.1.3 10.1.4 10.1.5 10.1.6 10.1.7 10.2 10.3 Opening 3D Structure Viewer . . . . . . . . . Changing 3D Structure Appearance . . . . . . Moving, Zooming and Spinning 3D Structure . Selecting Sequence Region . . . . . . . . . . . Selecting Models to Display . . . . . . . . . . Exporting 3D Structure Image . . . . . . . . . Working with Several 3D Structures Views . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100 101 104 104 105 106 106
Workow Designer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108 Chromatogram Viewer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109 10.3.1 10.3.2 Exporting Chromatogram Data . . . . . . . . . . . . . . . . . . . . . . . . . . 110 Viewing Two Chromatograms Simultaneously . . . . . . . . . . . . . . . . . . . 110
Circular Viewer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113 DNA Annotator . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116 DNA Export . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 117 10.6.1 10.6.2 10.6.3 10.6.4 10.6.5 10.6.6 Exporting Sequences as Alignments . . . . . . . . . . . . . . Export Sequences to FASTA, Genbank and FASTQ Formats . Exporting Alignment to Sequence Format . . . . . . . . . . . Exporting Selected Sequence Region . . . . . . . . . . . . . Exporting Sequence of Selected Annotations . . . . . . . . . Exporting Annotations to CSV Format . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 117 117 119 120 120 121
10.7 10.8
ORF Marker . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 125 Remote BLAST . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 126 Repeat Finder . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130 10.11.1 10.11.2 Finding Repeats . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130 Finding Tandem Repeats . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 131 Selecting Restriction Enzymes . . Using Custom File with Enzymes Filtering by Number of Hits . . . Excluding Region . . . . . . . . . Circular Molecule . . . . . . . . . Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 135 135 135 135 135 136
10.12
10.13
Molecular Cloning in silico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 137 10.13.1 10.13.2 10.13.3 Digesting into Fragments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 137 Creating Fragment . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138 Constructing Molecule . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
Secondary Structure Prediction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 142 SITECON . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144 Smith-Waterman Search . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 146 HMM2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 148 10.17.1 10.17.2 10.17.3 Building HMM Model (HMM Build) . . . . . . . . . . . . . . . . . . . . . . . 148 Calibrating HMM Model (HMM Calibrate) . . . . . . . . . . . . . . . . . . . . 149 Searching Sequence Using HMM Prole (HMM Search) . . . . . . . . . . . . . 150 Building HMM Model (HMM3 Build) . . . . . . . . . . . . . . . . . . . . . . . 151 Searching Sequence Using HMM Prole (HMM3 Search) . . . . . . . . . . . . 152 Searching Sequence Against Sequence Database (Phmmer Search) . . . . . . . 154 Aligning with MUSCLE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 156 Aligning Prole to Prole with MUSCLE . . . . . . . . . . . . . . . . . . . . . 157 Aligning Sequences to Prole with MUSCLE . . . . . . . . . . . . . . . . . . . 157 Aligning Short Reads . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 159 Building Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 160 Aligning Short Reads . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 161 Building Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 162 Searching JASPAR Database . . . . . . . . . . . . . . . . . . . . . . . . . . . 166 Building New Matrix . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 167 Creating Dotplot . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 169 Navigating in Dotplot . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 172 Zooming to Selected Region . . . . . . . . . . . . . . . . . . . . . . . . . . . . 173
10.18
10.19
10.20
10.21
10.22
10.23
Selecting Repeat . . . . . . . . . . . . . . . . . . . . . . . . . Interpreting Dotplot: Identifying Matches, Mutations, Ivertions, Editing Parameters . . . . . . . . . . . . . . . . . . . . . . . . Saving Dotplot as Image . . . . . . . . . . . . . . . . . . . . . Saving and Loading Dotplot . . . . . . . . . . . . . . . . . . . Building Dotplot for Currently Opened Sequence . . . . . . . . Comparing Several Dotplots . . . . . . . . . . . . . . . . . . .
. . . etc. . . . . . . . . . . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
Primer3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 178 External Tools 10.25.1 10.25.2 10.25.3 10.25.4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 179 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 179 182 184 184 BLAST/BLAST+ ClustalW . . . . . MAFFT . . . . . . T-Coee . . . . .
10.26
11
12
APPENDICES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 204
12.1 Appendix A. Supported File Formats . . . . . . . . . . . . . . . . . . . . . . . . . . 204
Specic File Formats . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 204 UGENE Native File Formats . . . . . . . . . . . . . . . . . . . . . . . . . . . . 207 Other File Formats . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 207
1 About Unipro
Established in 1992 Unipro company has its headquarters located in Novosibirsk Akademgorodok (the home of Siberian Branch of Russian Academy of Sciences). The companys primary activity is IT outsourcing solutions. To learn more about the company, please, visit the company website.
1.1 Contacts
Company website: Address: http://unipro.ru UniPro, 6/1 Lavrentiev Avenue 630090, Novosibirsk, Russia Marketing department: Tel: +7 (383) 3326061 Fax: +7 (383) 3302960 Email: [email protected] UGENE website: UGENE technical support: http://ugene.unipro.ru Email: [email protected]
2 About UGENE
Unipro UGENE is a free cross-platform genome analysis suite. It is distributed under the terms of the GNU General Public License. To learn more about UGENE visit UGENE website. It works on Windows, Mac OS X or Linux and requires only a few clicks to install.
10
2.4 Cooperation
Can be used for education purposes in schools and universities Features to be included into the next release are initiated by users UGENE team is ready for collaboration in related projects, both free and commercial
11
3 Installation
Get the appropriate package from the UGENE download page http://ugene.unipro.ru/download.html. Follow the installation instructions on the same page to install UGENE on your system. Quick guides on how to install UGENE on Windows, Linux and Mac OS X are situated below.
Alternatively, to use UGENE without installing: 1. Download UGENE zip package. 2. Unpack it. 3. Launch the ugeneui.exe le.
12
The downloaded le has *.tar.gz extension. 2. Unpack the archive. You can use this command: tar -xf [name of the downloaded *.tar.gz file] 3. Change the working directory to the unpacked UGENE directory: cd [name of the unpacked directory] 4. Launch the UGENE GUI version using the command: ./ugene -ui or the command line version using the command: ./ugene Note: Several native packages for specic Linux distributions are also available. Find out details on the download page. Note: UGENE is a part of Ubuntu and Fedora Linux distributions.
13
2. Launch the *.dmg le and accept the GNU license agreement. The following window will appear:
3. To start UGENE click on the ugeneui icon. You can also copy UGENE to the Applications folder by dragging it.
14
Chapter 3. Installation
4 Basic Functions
4.1 UGENE Terminology
Project Storage for a set of data les and visualization options. Document A single le (can be stored on a local hard drive or be a remote web page). Each document contains a set of objects. Object A minimal and complete model of biological data. For example: a single sequence, a set of annotations, a multiple sequence alignment. Task A process, usually asynchronous, that works in background. For example: some computations, loading and writing les. Plugin A dynamically loaded module that adds new functionality to UGENE. Object view A graphical view for a single or a set of objects. Project view A visual component used to manage active project. Task view A visual component used to manage active tasks. Log view A visual component used to show logs. Plugin viewer A visual component used to manage plugins. Sequence view An Object view aimed to visualize DNA, RNA or protein sequences along with their properties like annotations, chromatograms, 3D models, statistical data, etc. Annotation Additional information about a sequence, identied by its name and the sequence region. Alignment editor An Object view used to visualize and edit DNA, RNA or protein multiple sequence alignments. In the image below you can see a typical UGENE window with a Project view and a single Object view window opened:
15
Actions
Settings Tools
Window Help
The menus can be dynamically populated with new actions added by plugins. Check the Plugins documentation to learn how each plugin aects global and context menus.
16
Here you need to specify the visual name for the project and the directory and le to store it. After you click the Create button the Project view window is opened. It shows an empty project without any documents and views.
17
18
The following parameters are available: Document location location of the document. It can be a local le, a shared network le or a web reference, for example: C:\store\mydocument.gb \\192.168.0.3\store\mydocument.gb http://someaddress.com/store/mydocument.gb Document format species how to interpret the data stored in the le. As specied above the format is detected automatically, but you can select it manually. Force read-only mode locks the document for editing. Save le to disk before opening this option becomes available if a web reference has been specied in the Document location. Saves the remote le to the local disk before opening it. Custom settings the button is available for Genbank, EMBL, FASTA and FASTQ document formats. The button opens the following dialog:
If there are several sequences in the document, then selecting the Separate sequences option will open several sequences in a Sequence view window. Contrariwise, selecting the Merge sequences option will merge the sequences into one sequence. The Gap length parameter species the length of the gaps inserted between the merged sequences. Your choice will be saved as default if you check the Save as default settings check box. Note that if you select to merge the sequences, then the annotations of the sequences if any are also relocated automatically.
19
You can input the created sequence to the Paste data here eld: The following Custom settings are available: Alphabet here you can select the alphabet:
Skip unknown symbols / Replace unknown symbols with you can select either to skip unknown input symbols or to replace them with the specied symbol. Document location location of the created document.
20
Document format format of the created document. Currently available formats are FASTA and Genbank. Sequence name name of the sequence in the created document. Save le immediately check this option if you want to save the document immediately after the Create button is pressed. The created document will be added to the current project and opened in the Sequence view .
For other le formats (for example, PDB) the Save a copy context menu item is not available:
You can use the built-in export utilities to create a copy of the document in this case. The list of the le formats supported by UGENE can be found here. The list also species whether it is possible to read / write a le of a certain format using UGENE.
21
Below is the list of object types supported by the current version of UGENE. Object types: Symbol [3d] [a] [as] [c] [i] [m] [s] [t] [tr] Icon Description A 3D model. Annotations for DNA sequence regions. An assembly. Chromatogram data. A le with index information for a set of other, usually large les. A multiple sequence alignment. A nucleic, protein or raw sequence. A plain text. A phylogenetic tree.
You can edit names of particular objects, such as sequence objects, by selecting them in the Project view and then pressing F2. To be able to do so, the document containing the target object must be unlocked. To see the list of all available views for a given object select the object and activate the context menu inside the Project view window and select the Open view submenu:
22
The picture above illustrates an option to visualize the selected DNA sequence object using the Sequence view a complex and extensible Object view that focuses on visualization of sequence objects in combination with dierent kinds of related data: sequence annotations, graphs, chromatograms, sequence analysis algorithms. Note, that the Sequence view is described in more details in the separate documentation section.
For every persistent view UGENE automatically saves the state of the view in the Auto saved bookmark when the view is closed.
23
Now, by activating bookmarks you can restore the original view state. For example for the Sequence view bookmarks you can store a visual position and zoom scale for the sequence region.
Use the F2 keyboard shortcut to rename a bookmark. To remove a bookmark press the Delete key. UGENE has limited set of built-in Object views. Extensions modules or plugins can be used to adjust the existing views or to add new views to the tool.
To load a saved project later, select File Open and specify the path to the project le.
24
The window shows the list of available plugins. To add or remove plugins use the Add plugin and the Remove plugin items available in the Plugin viewer context menu:
When you select the Remove plugin item for a plugin, the plugins status is changed to the to remove after restart value. The Remove plugin is no more available in the context menu of the plugin. Instead the Enable plugin item appears in the context menu:
If you select this item the plugin will be enabled again, i.e. it will not be removed after restart. Otherwise, the plugin will not be available after UGENE restart.
25
Here you need to enter unique id of the biological object and choose a database. Unique identiers are dierent for various databases. For example, for NCBI GenBank such unique id could be Accession Number or NCBI GI number. Optionally, you can browse for a directory to save the fetched le to. After you click the OK button, UGENE downloads the biological object (DNA sequence, protein sequence, 3d model, etc.) and adds it to the current project. If something goes wrong check the Log view , it will help you to diagnose the problem.
26
4.12.1 General
The following settings are available on the tab: Language of User Interface (applied after restart) here you can select UGENE localization. Currently available localizations are EN, RU and CS. The default value (Autodetection) species that UGENE should use the operating system regional options to select the localization. This setting is applied only after UGENE is reopened. Appearance denes the appearance of the application, for example here is a part of the same dialog when the Cleanlooks appearance has been applied:
Preferred Web browser you can use either System default browser or specify some other browser. 4.12. UGENE Application Settings 27
Open last project at startup if the option is checked, the last project is opened when UGENE is started. Path to downloaded data species the path where les downloaded from the remote databases will be stored. Enable statistical reports collecting collects information about UGENE usage and sends it to the UGENE team to help improve the application. Note: The collected information includes: 1. System info: UGENE version, OS name, Qt version, etc. 2. Counters info: number of launches of certain tasks (e.g. HMM search, MUSCLE align). The collected information DOESNT include any personal data.
4.12.2 Resources
On the Resources tab you can set resources that can be used by the application: Optimize for CPU count, Tasks memory limit and Threads limit.
28
4.12.3 Network
On the Network settings tab of the dialog you can specify Proxy server parameters, select SSL settings and congure the Remote request timeout.
4.12.4 Logging
29
On the Log settings tab you can select type of log information (ERROR, INFO, DETAILS, TRACE ) for each Category that will be output to the Log view . You can select format for each log message by checking the Show date, Show log level and Show log category options.
4.12.6 OpenCL
If you have a video card that supports OpenCL you can use it to speed up some calculations in UGENE. To do it install the latest video card driver and check the corresponding check box:
30
Now you can, for example, use OpenCL optimization for the Smith-Waterman algorithm.
4.12.7 CUDA
If you have a NVIDIA video card that supports Compute Unied Device Architecture (CUDA) you can use it to speed up some calculations in UGENE. To do it install the latest video driver and check the corresponding check box:
Now you can, for example, use OpenCL optimization for the Smith-Waterman algorithm.
31
32
The Index location value will appear automatically after you choose the document location, but you can choose some other location. Optionally you can check the Compress option to archive the result index le with the gzip utility (it is checked by default). Open the created index le. For example, the screenshot below demonstrates an index le created from a FASTA le with several sequences:
As you can see there is an index for each sequence. You can export the indexed sequence by selecting the corresponding index and using the Export Save selection to a new le context menu item.
33
You can also lter the indexes using the green arrow pop-up menu:
34
5 Sequence View
5.1 Sequence View Components
The Sequence view is one of the major Object views in UGENE aimed to visualize and edit DNA, RNA or protein sequences along with their properties like annotations, chromatograms, 3D models, statistical data, etc. For each le UGENE analyzes the le content and automatically opens the most appropriate view. To activate the Sequence view open any le with at least one sequence. For example you can use the $UGENE/data/samples/EMBL/AF177870.emb le provided with UGENE. After opening the le in UGENE the Sequence view window appears:
After the view is opened you can see a set of new buttons in the toolbar area. The actions provided by these buttons are available for all sequences opened in the view. In the picture below these buttons are pointed by the "Global actions arrow. Below the toolbar there is an area for a single or several sequences. For each sequence a smaller toolbar with actions for the sequence and the following areas are available: Sequence overview Shows the sequence in whole and provides handy navigation in the Sequence zoom view and the Sequence details view . Sequence zoom view Provide exible tools for navigation in large annotated sequence regions. Sequence details view A supplementary component of the Sequence overview . It is used to show sequence content without zooming. Annotations editor Contains tools to manipulate annotations for a sequence.
35
You can change the focus by clicking on the corresponding sequence area. All sequences that are not in focus have the sequence name and icon disabled. The bottom area of the Sequence view is the Annotations editor . It contains a tree-like structure of all annotations available for all sequences shown in the Sequence view and can be used to perform various actions on annotations: create a new annotation, modify the existing one, group, sort, etc. 36 Chapter 5. Sequence View
The global action toolbar provides possibility to go to the specied position (in all sequences at the same time). Also it allows to lock or adjust ranges of sequences in the same Sequence view . See this paragraph for details.
See also: Toggling Views Capturing Screenshot Zooming Sequence Showing and Hiding Translations Selecting Sequence
37
When the sigma button (in the right part of the Sequence overview ) is pressed, density of annotations in the sequence is shown. For example in the picture below there are annotations in the parts of the sequence that are marked with dark grey color:
38
Below the annotation rows there is a ruler to show coordinates in the sequence.
When the Show All Rows item is checked all available annotations are always shown. You can also add rows by selecting the +5 Rows and +1 Row items and remove rows by selecting the -5 Rows and -1 Row items. To restore the default number of rows select the Reset Rows Number item. See also: Navigating Sequence zoom view using Sequence overview Zooming Sequence Creating New Ruler Manipulating Annotations
39
See also: Navigating the Sequence details view using the Sequence overview Selecting Amino Translation Showing and Hiding Translations
40
Or use the Go to position context menu or the Actions main menu item.
The sequence can be removed from the view using the same menu. Once you remove the last sequence in the view, the view is automatically closed.
41
There are standard Zoom In and Zoom Out buttons. Additionally you can zoom to a selected region using the Zoom to Selection button. To restore the default view of the Sequence zoom view (when the sequence is not zoomed) use the Zoom to Whole Sequence button.
The new ruler will be shown right above the default one:
42
43
Selecting the Sequence region context menu item opens the Select range dialog:
Here you can specify the sequence range you would like to select. You can open the same dialog using the Select sequence region button on a sequence toolbar or using the Ctrl-A key sequence. To use the Sequence between selected annotations item, select two annotations in the Annotations editor (holding the Ctrl key at the same time):
44
And select the Select Sequence between selected annotations item in the context menu. The Sequence around selected annotations item selects the selected annotations and the sequences between these annotations.
45
Another way to select a sequence around annotations is to hold Shift and Ctrl keys while clicking on the annotations either in the Sequence details view or in the Sequence zoom view .
2. Using the following shortcuts: Ctrl-C copies direct sequence strand Ctrl-T copies direct amino translation Ctrl-Shift-C copies reverse-complement sequence Ctrl-Shift-T copies reverse-complement amino translation 3. Using the Copy submenu of the context menu:
46
The following parameters are available: Pattern species the pattern to search in the sequence. This parameter is mandatory. Search in where to search the pattern: in the Sequence or in its Translation. Strand where to search the pattern: in the Direct strand, in the Complement strand or in Both (direct and complement) strands. Algorithm if the Match percent is 100% (by default) the pattern is searched exactly in the sequence, otherwise the partial match is accepted. The Substitute option means that subsequences are searched for, where symbols can be substituted by other symbols. Example: If the pattern is AAAAB, the Match percent is 80% and the sequence is ...AAAAC..., then the AAAAC subsequence will be found by this algorithm as: 1. The last symbol B in the pattern is substituted by the C symbol in the subsequence. 2. The subsequence and the pattern dier by 20% only. The InsDel option means that subsequences are search for, where symbols can be inserted or deleted. Example: If the pattern is AAAAB, the Match percent is 80% and the sequence is ...AAAA..., then the AAAA subsequence will be found by this algorithm as: 1. The last symbol B in the pattern is deleted in the subsequence. 2. The subsequence and the pattern dier by 20% only. Find next nds the next result only. Find all nds all suitable results. Range start / Range end species the sequence range where to search for. 5.7. Manipulating Sequence 47
Position can be used with the Find next button, species from which position to search the next result. Has no meaning when the Find all button is pressed. Sequence range, Selection range appropriately change the Range start and the Range end parameters. Clear results becomes available when some results have been found, clears these results. Remove overlaps becomes available only if there are results that are fully contained in some other results, removes these results. Save as annotations becomes available when some results have been found, opens the Create annotation dialog to save the result as annotations to the sequence.
The Edit sequence submenu is available in the Actions main menu and in the Sequence view context menu. When you select the Insert subsequence item the following dialog is opened:
48
Description of the dialog parameters: Paste data here you must input the inserted subsequence. This parameter is mandatory. Annotated regions resolving denes either to Resize, Remove or Split an annotation (into two annotations) in case when the subsequence is inserted to the sequence position where some annotations are presented. Start position the sequence position where to insert the subsequence. Save resulted document to a new le the result sequence can be saves to a new le instead of modifying the current le. You must select the Document location. FASTA and Genbank le formats are available when you do not include annotations to the result le. If you check the Merge annotations to this le item, the annotations will also be saved to the result le (Genbank le format is only available in this case). In case a subsequence has been selected, the rst item in the Edit sequence submenu is called Replace subsequence instead of Insert subsequence. The dialog opened in this case is similar to the dialog described above, except it already contains the sequence to be edited an doesnt allow to input the start position.
Also it is possible to remove selected subsequence from a sequence. When you select corresponding item (in the context menu or in the Actions menu), the Remove subsequence dialog appears:
Description of the parameters: Region to remove species the region of the sequence that will be removed in the form. This parameter is mandatory. Annotated regions resolving species what to do with annotations that overlap with the region that is removed. You can select either Resize such annotations (i.e. make it smaller) or Remove them. Save resulted document to a new le similar to the same parameter in the Insert subsequence dialog (described above).
49
To unlock the scales click the same button again. You may use the Adjust scales button to synchronize scales without locking them. Note, that if you have a selected sequence region or a selected annotation the scales will be synchronized by the start position of the region or the annotation. If there are no active selection the regions are synchronized by the rst visible sequence position on the screen.
50
There are usually several objects with annotations in the Annotations editor . A special Auto-annotations object is always presented for each sequnce opened. It contains annotations automatically calculated for the sequence (see below for details). An object contains groups of annotations used by UGENE for logical organization of the annotations. An annotation must always belongs to some group.
51
For documents created not by UGENE annotations are grouped by their names. For annotations created in UGENE it is possible to use arbitrary group names. Groups can contain both annotations and other groups. The numbers in the brackets after a group name in the Annotations editor are the count of subgroups and annotations in the current group. A single annotation is allowed to be presented in several groups simultaneously. An annotation is physically removed from the document when it does not belong to any group.
To disable/enable the automatic annotations calculations use the Automatic Annotations Highlighting menu button on the Sequence view toolbar:
When you click on the value a web page is opened or a le is loaded specied in the reference. The loaded le is added to the current project.
52
The dialog asks where to save the annotation. It could be either an existing annotation table object or a new document (le). You can also specify the name of the group and the name of the annotation. If the group name is set to <auto> UGENE will use the annotation name as the name for the group. You can use the / characters in this eld as a group name separator to create subgroups. The Location eld contains annotation coordinates. The coordinates must be provided in the Genbank or EMBL le formats. If you want to annotate complement sequence strand surround the coordinates with the complement() word or press the last button in the Location row to do it automatically.
53
Note, that by default the Location eld contains the coordinates of the selected sequence region. Once the Create button is pressed the annotation is created and highlighted both in the Sequence overview and the Sequence details view areas:
54
To select a color of annotations in a group click on the corresponding color box in the dialog. The Select Color dialog will appear. The color selected will be applied to all the annotation in the group. To show annotation on the translation of the sequence check the Amino check box, for example:
The Highlight check box enables/disables the highlighting. You can also enable/disable this feature using the annotation context menu:
In the Visual Qualier column you can input a qualier name. The value of this qualier will be displayed on the annotation in the Sequence zoom view . For example if a note value is written for the misc_feature in the Visual Qualier column, then the value is displayed on the annotation (see the image below).
55
Note: If you input several qualiers names to the Visual Qualier (separated by comma), then the rst found qualier is taken into account.
56
Here you can specify the name and the value of the qualier. You can use the F2 key to rename a qualier:
To edit a qualier, select the qualier and press the F4 key or use the Edit qualier context menu item:
57
Basically you need to specify the le to read annotations table from (required):
58
And the format of and the path to the le to write the annotations table into (required):
Check Add result le to project to link the annotations to the currently opened sequence.
To use a separator to split the table, check the Column separator item and specify the separator symbols. Also you can press Guess to try to detect the separator from the input le.
Alternatively, you can press Edit and edit the script which will specify the separator for each parsed line. It is possible to use line number in the script.
Using the arrows, you exclude the necessary number of lines at the beginning of the document from parsing. You can also skip all lines that start with the specied text.
By pressing Preview one can bring up the view of the current annotations table (which is produced from the input le with the specied parameters values). The input le contents will also be shown at the bottom part of the dialog.
59
The preview table headline indicates the types of the information contained in the corresponding columns. By default the values are [ignored]. To specify a column role, click on the corresponding headline element:
The annotation start and end positions must be specied. It is possible to add an oset to every read start position by checking the Add oset checkbox, and to shorten annotations by one from the end by uncheking the Inclusive checkbox. When all the roles are specied, press Run. With the Add to project checkbox specied and a sequence view opened, on success you will see the sequence view with annotations linked:
60
61
The 3D Structure Viewer plugin adds 3D visualization for PDB and MMDB les:
62
The Chromatogram Viewer plugin adds support for chromatograms visualization and editing:
The DNA Graphs Package plugin shows various graphs for sequences:
63
A number of other plugins add graphical interface for popular sequence analysis methods:
64
6 Alignment Editor
6.1 Overview
This chapter gives an overview of the Alignment editor components and explains basic concepts of browsing an alignment.
65
The Alignment editor components: Sequence area This is the main component of the editor. It displays aligned sequences. The upper part of the Sequence area is the ruler, which shows the coordinates of the currently visible row sequences. Consensus area This component is situated above the Sequence area. It shows the consensus sequence for the current alignment calculated using currently selected algorithm. Sequence list This component is located in the left part of the Sequence area. It shows names of the corresponding sequences in the alignment. Editor toolbar The toolbar contains shortcuts for important editor actions, such as Undo/Redo, Zooming and others. Sequence offsets These are the osets for the rst and the last visible base for each alignment row. Note that the oset value doesnt include gaps. For example, lets assume that the coordinate of the rst visible base of the row is N, but the row contains K gaps before the position N. The starting oset value will be N-K. The same rule is true for the ending oset. You can turn o the Sequence osets by unchecking the Actions View Show osets main menu item or View Show osets context menu item. Global coordinates This component displays the coordinates of the upper left corner of the current selection. If no region is selected it shows the starting alignment point.
66
Alignment lock status As in the Sequence view this component shows whether the alignment is locked. Locked documents are not allowed to be modied.
6.1.3 Navigation
The Sequence area provides several exible ways to navigate through an alignment. The simplest way is to use the mouse and the scrollbars. Alternatively you can use arrow keys on the keyboard to navigate. The list of hot keys for quick navigation: PageUp to move one screen left. PageDown to move one screen right. Home to center the starting columns of the alignment. End to move to the trailing columns of the alignment Hint: if you use Shift key with the hot keys above you will navigate through the rows. For example, Shift-PageDown will move one screen down. Finally you can use the Go to position dialog from the Actions menu, the context menu or the editor toolbar.
Enter the column number (base coordinate) and the view will be centered to the corresponding base.
6.1. Overview
67
By default, the base characters are visible when zooming. But for rather long sequences there is another zoom mode available. In this mode the bases are not shown. This allows viewing very large sequence regions (up to 500 bp).
You can zoom to the selected region by clicking the Zoom to selection button. It is very convenient operation, when the alignment size is rather large. For example, you can zoom out to some percentage, select an interesting region and then zoom to the selection. You can change font by clicking the Change font button. To reset zoom and font click the Reset zoom button.
Press the right arrow to search in the direction "From left to right, from top to bottom". Press the left arrow to search in the direction "From right to left, from bottom to top". If the pattern is found, the result will be focused and highlighted in the Sequence area. You can continue the search in any direction from this position.
68
6.1.7 Consensus
Each base of a consensus sequence is calculated as a function of the corresponding column bases. There are dierent methods to calculate the consensus. Each method reveals unique biological properties of the aligned sequences. The Alignment editor allows switching between dierent consensus modes. To switch the consensus mode activate the context menu (using the right mouse button) or the Actions menu and select the Consensus mode item.
There are several modes: JalView (Default) it is based on the JalView algorithm. Returns + if there are 2 characters with high frequency. Returns symbol in lower case if the symbol content in a row is lower than the specied threshold. ClustalW emulates the ClustalW program and le format behavior. Levitsky this algorithm is proposed by Victor Levitsky to calculate consensus of DNA alignments. At rst, it collects global alignment frequencies for every symbol using extended (15 symbols) DNA alphabet. Then, for every column it selects the rarest symbol in the whole alignment with percentage in the column greater or equals to the threshold value. Strict the algorithm returns gap character () if symbol frequency in a column is lower than the threshold specied.
6.1. Overview
69
70
Extracting Selected as MSA It is possible to extract a subalignment and save it as new multiple sequence alignment (MSA). Select a subalignment and choose the Edit Extract selected as MSA item in the Actions main menu or in the context menu. The following dialog appears:
Specify the name of the new MSA le in the File name eld. The currently selected region is extracted by default when you press the Extract button. You can change the columns to be extracted using the From and to elds. And change the rows to be extracted by checking / unchecking required sequences in the Selected sequences list. Use buttons: Invert selection to invert the selection of the sequences. Select all to select all sequences. Clear selection to clear the selection of all sequences. The Add to project check box species to add the MSA le created from the subalignment to the active project. Removing All Gaps Use the Edit Remove all gaps item in the Actions main menu or in the context menu to remove all gaps from the alignment. Removing Selection To remove a subalignment select it and choose the Edit Remove selection item in the context menu or press the Delete key.
71
Removing Columns of Gaps To remove colums containg certain number of gaps select the Edit Remove columns of gaps item in the context menu. The dialog appears:
There are the following options: Remove columns with number of gaps removes columns with number of gaps greater than or equal to the specied value. Remove columns with percentage of gaps removes columns with percentage of gaps greater than or equal to the specied value. Remove all columns of gaps this option is selected by default. It species to remove columns from the alignment if they entirely consist of gaps. Select the option required and press the Remove button. Filling Selection with Gaps Select a region in the alignment and choose the Edit Fill selection with gaps item in the context menu or press the Spacebar. The region is lled with gaps shifting the subalignment from the region to the right.
72
You will see the Project view tree ltered to show only appropriate sequences. Select the items to add and press the Ok button. Copying Sequences To copy current selection click the Copy Copy selection item in the Actions main menu or the context menu. The hotkey for this action is Ctrl-C. To copy one or several sequences do the following: Select the sequences in the Sequence list area; Select the Copy Copy selection context menu item in the Sequence area or use hot key combination. Note, that if you activate context menu in the Sequence list area you will lose your current selection.
73
To copy consensus sequence use the Copy Copy consensus item. Sorting Sequences To sort sequences by name in the alphabetical order choose the View Sort sequences by name item from the Actions main menu or the context menu.
74
The le save dialog will appear where you should set name, location and format of the picture. UGENE supports export to the PNG, TIFF and JPEG image formats.
75
The following parameters are available: Tree building method method to build a tree. Currently only the neighbour-joining method is supported. The PHYLIP package implementation of the method is used under the hood. Distance matrix model model to compute a distance matrix. The following values are available for a nucleotide multiple sequence alignment: F84 Kimura Jukes-Cantor LogDet The following models are available for a protein alignment: Jones-Taylor-Thornton Heniko/Tillier PMB Dayho PAM Kimura 76 Chapter 6. Alignment Editor
Gamma distributed rates across sites species to take into account unequal rates of change at dierent sites. It is assumed that the distribution of the rates follows the Gamma distribution. Coecient of variation of substitution rate among sites becomes available if the Gamma distributed rates across sites parameter is checked. Species the coecient of the distribution of the rates. Transition/transversion ratio expected ratio of transitions to transversions. To enable bootstrapping check the Bootstrapping and Consensus Trees group check box. The following parameters are available: Number of replicates number of replicate date sets. Seed random number seed. By default, it is generated automatically. You can manually change this value in order to make results of dierent runs (of a tree building) reproducible. The should must be an integer greater than zero and less than 32767 and which is of the form 4n+1, that is, it leaves a remainder of 1 when divided by 4. Any odd number can also be used, but may result in a random number sequence that repeats itself after less than the full one billion numbers. Usually this is not a problem. Consensus type species the method to build the consensus tree. Select one of the following: Strict species that a set of species must appear in all input trees to be included in the strict consensus tree. Majority Rule (extended) species that any set of species that appears in more than 50% of the trees is included. The program then considers the other sets of species in order of the frequency with which they have appeared, adding to the consensus tree any which are compatible with it until the tree is fully resolved. This is the default setting. M1 includes in the consensus tree any sets of species that occur among the input trees more than a specied fraction of the time (see the Fraction parameter below). The Strict consensus and the Majority Rule consensus are extreme cases of the Ml consensus, being for fractions of 1 and 0.5 respectively. Majority Rule species that a set of species is included in the consensus tree if it is present in more than half of the input trees. Fraction becomes available when the Consensus type parameter is set to M1. Species the fraction. Save tree to le to save the tree built. Press the Build button to build a tree with the parameters selected.
77
7 Assembly Browser
The UGENE Assembly Browser project started in 2010 was inspired by Illumina iDEA Challenge 2011 and multiple requests from UGENE users. The main goal of the Assembly Browser is to let a user to visualize and eciently browse large next generation sequence assemblies. Currently supported format is BAM. It is a binary version of the SAM (Sequence Alignment/Map) format. The both formats are produced by SAMtools and described in the following specication: http://samtools.sourceforge.net/SAM1.pdf. The SAM format support is scheduled to the next UGENE iteration. Support of other formats is also planned, so please send us a request if youre interested in a certain format. To browse an assembly data in UGENE a BAM le should be imported to a UGENE database le. This approach has both advantages and disadvantages. The disadvantages are that the import may take time for a large BAM le and there should be enough disk space to store the database le. On the other hand this allows to overview the whole assembly and navigate in it rather rapidly. In addition during the import you can select contigs to be imported from the BAM le. So there is no need to import the whole le if youre going to work only with some contigs. Note that in future it is planned to support another approach as well when a BAM le is opened directly. The Assembly Browser has been tested on dierent BAM les from the 1000 Genomes Project and other sources. Read the documentation below to learn more about the Assembly Browser features.
The Source URL eld in the dialog species the le imported. The Info button nearby can be used to obtain additional information about the le.
78
There is a list of contigs below the Source URL. Check the contigs that you want to import to the database. You can use the Select All, Deselect All and Invert Selection buttons to manage selection. The Destination URL eld species the output database le. The Add to project check box adds the imported le to the current project. To start the import click the Import button in the dialog. You can see the progress of import in the Task view .
Each [as] object corresponds to a contig imported. When you double-click on an [as] object a new Assembly Browser window with the assembly data is opened. A window for the rst assembly object in the list is opened automatically after the import. The window appeared contains the list of well-covered regions of the assembly:
Note that for large assemblies it may take some time to calculate the overview and the well-covered regions.
79
To see the reads either select a region from the list or zoom in, for example by clicking the link above the well-covered regions or by rotating the mouse wheel.
80
By dragging the mouse while holding the left mouse button you can navigate in the Reads Area. To navigate long distances in the Reads Area use the Assembly Overview described below . Other ways to navigate in the assembly are: To use horizontal and vertical scroll bars of the Reads Area. Input a certain location on the toolbar and click the Go! button.
If you hold Shift and select a region on the overview, the overview is zoomed to the selection. Note that when the Assembly Overview is in focus and you use either the zoom buttons on the toolbar, the zoom items in the Actions main menu or a mouse wheel the Reads Area is resized appropriately. The Assembly Overview can also be resized. In this case it shows a coverage overview of some part of the assembly. To zoom in the overview select either the Zoom in or the Zoom in 100x item in the Assembly Overview context menu. You can scroll the resized overview by dragging the mouse while holding the mouse wheel. To zoom out the overview select the Zoom out item in the context menu. The Restore global overview item in the context menu restores the default overview size when the whole contig overview is shown. Notice that the Assembly Overview shows the coordiantes of the assembly areas visible in the Reads Area and in the Assembly Overview:
To scroll the resized overview drag the mouse while holding the mouse wheel. To learn about hotkeys available refer Assembly Browser Hotkeys paragraph.
81
To show/hide the coordinates on the ruler you can click the following button on the toolbar:
It shows the information about the read obtained from the: Read name Location Length Cigar Strand Read sequence Note: The hint may also contain the Unmapped label. This happens when a valid location was specied for a read in the imported le, but at the same time the read was marked as unmapped. The operations in the Cigar parameter are described as follows: M - Alignment match (can be a sequence match or mismatch). I - Insertion to the reference. Skipped when the read is aligned to the reference, i.e. it is not shown in the Reads Area, but is present in the read sequence. D - Deletion from the reference. Gaps are inserted to the read when read is aligned to the reference. For example:
N - Skipped region from the reference. Behaves as D, but has dierent biological meaning: for mRNA-togenome alignment represents an intron. S - Soft clipping (clipped sequences are present in the read sequence, i.e. behaves as I). 82 Chapter 7. Assembly Browser
H - Hard clipping (clipped sequences are not present in the read sequence). P - Padding (silent deletion from padded reference). = - Exact match to the reference. x - Reference sequence mismatch. To copy the information about the read to the clipboard select the Copy read information to clipboard item in the Reads Area context menu. Now you can paste it in any text editor.
To remove the association select the Unassociate item in the Reference Area context menu.
83
7.5 Exporting
7.5.1 Exporting Read
To export a read right-click on it in the Reads Area and select the Export Current read item in the context menu. The Export Reads dialog appears:
Select a le to export the read to and the le format. The read can be exported either to a FASTA or FASTQ le. When the parameters are set click the Export button. The read is exported to the le and if the Add to project check box has been checked it is added to the current project from where you can open it.
In the dialog you can select the image le name and its format (bmp, jpeg, png, etc.). For some le formats the Quality parameter also becomes available. When the parameters are set click the OK button. 84 Chapter 7. Assembly Browser
85
To load a tree from a le follow the instruction described in the Opening Document paragraph. For example, you may open the $UGENE\data\samples\Newick\COI.nwk sample le provided within UGENE package. To build a tree from a multiple sequence alignment see the Building Phylogenetic Tree paragraph. To learn what you can do with a tree using UGENE Phylogenetic Tree Viewer read the documentation below.
86
In the dialog you can tune the width of the tree. If the tree layout is set to rectangular you can tune the height of the tree also. And you can select the tree view: Phylogram Cladogram
Here you can select the color and the line width of the tree branches. Note that when a clade has been selected the branch settings are applied to the clade only.
87
88
Here you can select color, font, size and attributes (bold, italic, etc.) of the labels. Note that when a clade has been selected the labels formatting settings are applied to the clade only.
You can see that the corresponding branches are highlighted. To select several clades at the same time hold the Shift key and click on the root nodes of the clades.
To show the collapsed clade select the Expand item in the nodes context menu.
90
91
9 Distributed Computing
Distributed computing allows to notably increase the performance of computational tasks by distributing the task data among computational units. However the distributed computing assumes complex solutions: specialized versions of algorithms, network communication etc. Unipro UGENE project provides advanced distributed computing capabilities. Despite the complexity of the internal structure, for users running computational tasks on a remote machine is as easy as running it on a local machine. Starting with version 1.7.2 UGENE supports cloud computing. For example, computational workows can be launched on the Amazon EC2 cloud. Check for details the following documentation section: Running Workows on Cloud There are also several distributed algorithms that can be executed on a remote machine: HMMER3 search Smith-Waterman search Muscle3 align
The Remote machines monitor allows you to add, remove or modify remote machines.
92
To add a new remote machine, click the Add button. In the appeared dialog select the protocol and ll other required elds:
Modication of a remote machine is as simple as adding a new one. Just select the machine and click the Modify button. To remove a remote machine from the monitor select the machine in the table and click the Remove button. You can Ping a remote machine to check if its still alive and UGENE is still running there. Some network protocols (for example, direct socket protocol) can do scanning of local network. To search for running UGENEs through such protocols click the Scan button. Also, you can use one of the public UGENE machines to run your tasks on it. To add public machines to the monitor click the Get public machines button.
93
94
Check the cloud machine status in the Remote machines monitor. If the session has been initialized successfully, the green tick is highlighted in the Ping column. If there is no green tick, check the Log view for details about the occurred problem.
Run the schema, e.g. by clicking the corresponding Workow Designer toolbar button. The Remote machines monitor dialog will appear. Select the remote machine that represents the EC2 service and click the Run button.
95
96
In the appeared HMM3 search dialog ll required parameters and click the Remote run... button.
The Remote machines monitor dialog will appear. You can also add, remove or modify remote machines here. Select a machine to run and click the Run button. Note that only 1 machine can be selected in the current version of UGENE. Thats all. After the task is nished you will see the task report in the Task view .
97
You will see the Smith-Waterman search dialog. Fill required elds and click the Remote run button.
The Remote machines monitor dialog will appear. You can also add, remove or modify remote machines here. Select a machine to run and click the Run button. Note that only 1 machine can be selected in the current version of UGENE. Thats all. After the task is nished you will see the task report in the Task view .
98
You will see the Align with MUSCLE dialog. Fill required elds and click the Remote run button.
The Remote machines monitor dialog will appear. You can also add, remove or modify remote machines here. Select a machine to run and click the Run button. Note that only 1 machine can be selected in the current version of UGENE. Thats all. After the task is nished you will see the task report in the Task view .
99
10 Plugins
10.1 3D Structure Viewer
The 3D Structure Viewer is intended for visualization of 3D structures of biological molecules. Using the 3D Structure Viewer you can work with data from the Protein Data Bank (PDB) - a repository for the 3D structural data of large biological molecules, such as proteins and nucleic acids, maintained by the Worlwide Protein Data Bank (wwPDB). You can work as well with data from the NCBI Molecular Modeling DataBase (MMDB), also known as "Entrez Structure", a database of experimentally determined structures obtained from the RCSB Protein Data Bank. Find the description of the 3D Structure Viewer features below.
100
Notice the Links button on the toolbar. When you click the button the menu appears with quick links to online resources with detailed information about the molecule opened: PDB Wiki RSCB PDB PDBsum NCBI MMDB Note that if youre online, you can access the Protein Data Bank directly from UGENE and load a required le by its PDB ID (see Fetching Data from Remote Database for details). Hint: Dont forget to select the correct database (PDB) while fetching.
101
Selecting Coloring Scheme You can select one of the following coloring schemes: Chemical Elements Molecular Chains Secondary Structure To change the coloring scheme open the Coloring Scheme menu (available in the context menu and in the Display menu on the toolbar).
Calculating Molecular Surface To calculate the molecular surface of a molecule select the Molecular Surface item in the 3D Structure Viewer context menu or in the Display menu on the toolbar and check one of the following items: SAS (solvent-accessible surface) SES (solvent-excluded surface) vdWS (van der Waals surface) To remove the molecular surface that has already been calculated select the O item. You can also select the Molecular Surface Render Style to modify the calculated molecular surface appearance: Convex Map Dots
102
Selecting Background Color To change the background color open the Settings dialog (choose the Settings item in the 3D Structure Viewer context menu or in the Display menu on the toolbar), press the Set background color button and select a color in the dialog appeared. Selecting Detail Level To select the detail level of a 3D Structure representation open the Settings dialog of the 3D Structure Viewer and drag the Detail level slider. Enabling Anaglyph View UGENE allows you to view a molecule in the anaglyph mode. To enable the anaglyph view open the Settings dialog of the 3D Structure Viewer and check the Anaglyph view check box. You can modify the color settings: select one of the available Glasses colors or set custom colors, swap the colors. The oset of the color layers can be adjusted by dragging the Eyes shift slider.
103
You can also overview the whole structure by spinning it automatically. Select the Spin item either in the 3D Structure Viewer context menu or in the Display menu on the toolbar to do it. To stop the spinning uncheck the Spin item.
To adjust the shading drag the Unselected regions shading slider in the Settings dialog.
To show all the models check the Select All item. To show only one model check the Exclusive item and then check the model you want to display. To show several models uncheck both the Select All and the Exclusive items and check the models you would like to display.
105
Here you can browse for the le name, select the width and height of the image as well as its format: svg, png, ps, jpg or ti.
2. Press the Add button on the toolbar. The Select Item dialog will appear. Select [3d] objects to add. Hint: Use the Ctrl keyboard button to select several objects.
106
Below you can see the 3D Structure Viewer with two views:
To select an active view click on the view area or select an appropriate value in the Active view combo box on the toolbar. To synchronize the views press the Synchronize 3D Structure Views sticky button on the toolbar (see the image above). When the button has been pressed the 3D structures are moved, zoomed and spinned synchronously. Press the button again to stop the views synchronization. The views that are no more required can be closed by selecting the Close button in the 3D Structure Viewer context menu. Also you can hide/show views for a while. Use the menu of the green arrow button on the toolbar to do it:
Notice that the 3D Structure Viewer can be closed from this menu.
107
To learn more about the Workow Designer read the Workow Designer Manual (follow the link on the UGENE documentation page).
108
To edit a sequence data, right-click on the chromatogram view and select the Edit new sequence item in the appeared context menu. The original DNA sequence is not allowed to be changed; however you can add and modify a new sequence stored in a separate le. The sequence being edited is displayed right above the original one. Symbols can be changed by clicking on interesting value, modications are shown in bold.
109
After clicking on the item, the Export chromatogram le dialog will appear:
Check the Reversed and Complemented options if you want to create a reverse and complement chromatogram. Press the Export button. The exported le will be opened in the Sequence view .
110
You can also use the Lock scales and Adjust scales global actions for the chromatograms. For example if you lock the scales you are able to scroll the sequences simultaneously. Also when you select a 10.3. Chromatogram Viewer 111
sequence region in one sequence, the same region is selected in the second sequence.
112
Pressing the button will show the circular view of the sequence:
Note: The Circular Viewer is opened automatically when the Sequence view is opened for a plasmid. The inner circle represents the sequence clockwise and the scale marks show the corresponding sequence positions. The sequence annotations are represented as curved colored regions at the outer side of the circle. The Circular Viewer helps to navigate within the sequence. You can select an annotation on the circular view and the annotation will also be focused and highlighted in all Sequence view areas: Sequence overview , Sequence zoom view , Sequence details view and Annotations editor . 10.4. Circular Viewer 113
This will also aect the Sequence view . Note that the circular view is zoomed automatically when the Circular Viewer area is resized:
114
So you can adjust it to an appropriate size. It is possible to rotate the circular view using the mouse wheel. Use the Export Save circular view as image context menu or the Actions main menu item to save the image of the circular view.
Dierent le formats are available, including *.png, *.bmp, *.jpg, *.svg and *.pdf. Note, that if a sequence le contains several sequences it is possible to view the circular views of the sequences in the same Circular Viewer area.
You can work with these circular views at the same time.
115
Using this dialog you can search for DNA sequence regions that contain every annotation from the list on the left side. The found regions are displayed on the right side of the dialog. Use the Save regions as annotations... button to store the regions as new annotations to the sequence.
116
The Export sequences as alignment dialog will appear where you can point the result alignment le location, select a multiple alignment le format and optionally add the created document to the current project:
Here you can select the location of the result le and a sequence le format (FASTA, Genbank or FASTQ). Also it is possible to specify whether to merge the exported sequences into a single sequence or store them as separate sequences. You can choose to add newly created document to the current project. If you merge the sequences, youre allowed to select the Gap length. This is the length of the insertion region between sequences that contain N symbols for nucleic or X for protein sequences. 118 Chapter 10. Plugins
Here it is possible to specify the result le location, to select a sequence le format, to dene whether to keep or remove gaps ( chars) in the aligned sequences and optionally add the created document to the current project.
119
The Export selected sequence region dialog will appear which is similar to the Export sequences dialog described above.
120
The Export sequence of selected annotations will appear which is similar to the Export sequences dialog described above.
Here you can set the path to the CSV le and optionally save the sequence along with annotations.
121
A new area with a graph will appear right above the Sequence zoom view :
All graphs are always aligned to the range shown in the Sequence zoom view. It means that if you change the visible range in overview (either by zooming or scrolling) the graph will also be updated. The most of the graphs use window averaging algorithms to calculate point values. Use the Graph settings context menu item to change the Window and Steps per window settings. You can enable multiple graphs by clicking corresponding Graphs menu items. To close the graph, uncheck its item in the Graphs menu. Note that some graphs like GC content (%) are available only for nucleic sequences.
122
Here is a brief description of the options that can be set in the dialog: Prole mode: Counts/Percents select the Percents to have scores shown as percents in the report. Show scores for gaps check this item if you want gap characters () statistics to be shown in the report. 10.8. DNA Statistics 123
Show scores for symbols not used in alignment if a symbol is not used in the alignment at all it wont be shown in the report. Check this item to make all symbols of alignment alphabet reported. Skip gaps in consensus position increments consensus ruler conguration. If checked the gaps in consensus will not lead to ruler increments. Save prole to le allows to save prole to a le in the HTML or CSV format. The CSV format is convenient for further processing in worksheets editors like Excel. The result prole in the HTML mode:
124
Notes: DNA-to-Amino translation table denes the way start and stop codons are encoded. ORFs with length lower than Min length value will not be shown in the report. Checking the Must start with init codon item switches ORF marker algorithm to the mode when any non-stop amino acid code is interpreted as region start position. After the search is nished you can browse the results, sort them by length, strand or start position and save as annotations to the original sequence in the Genbank format. 1.. only:: latex
125
The following dialog will appear where you can choose the search options:
126
General options are: Select search in the remote databases the blastn search is used for nucleotide sequences, blastp and cdd searches are used for amino sequences. UGENE also provides a way to use blastp and cdd searches for nucleotide sequences. This is achieved by translating the nucleotide sequence into the amino sequences. When a sequence is translated the translation table from the active Sequence view is used. Finally, all 6 translations are used to query the remote database with the selected blastp or cdd search. Expectation value this option species the statistical signicance threshold for reporting matches against database sequences. Lower expect thresholds are more stringent, leading to fewer chance matches being reported. Max hits the maximum number of hits that will be shown (not equal to number of annotations). Database the target database. Search for short, nearly exact matches automatically adjusts the word size and other parameters to improve results for short queries. Megablast select this option to compare query with closely related sequences. It works best if the target percent identity is 95% or more, but it is very fast. Search timeout sometimes a database doesnt respond, therefore you need to re-wait for the response. This option sets the time that will be spent for re-appeal to the database. Note that in case of long sequences time for request preparation increases and the search takes several minutes. Also there is Advanced options tab:
127
The view of the Advanced options tab depends on the selected search. For the blastn search it looks like on the picture above. Word size the size of the subsequence parameter for the initiated search. Gap costs costs to create and extend a gap in an alignment. Increasing the Gap costs will result in alignments which decrease the number of Gaps introduced. Match/Mismatch scores reward and penalty for matching and mismatching bases. Filters lters for regions of low compositional complexity and repeat elements of the humans genome. Masks for lookup table only this option masks only for purposes of constructing the lookup table used by BLAST so that no hits are found based upon low-complexity sequence or repeats (if repeat lter is checked). Mask lower case letters with this option selected you can cut and paste a FASTA sequence in upper case characters and denote areas you would like ltered with lower case. When the blastp search is selected in the general options, the view of the Advanced options tab is the following:
As you can see there is no Match scores option, but there are Matrix and Service options. Matrix key element in evaluating the quality of a pair-wise sequence alignment is the "substitution matrix", which assigns a score for aligning any possible pair of residues. Service blastp service which needs to be performed: plain, psi or phi.
128
The Advanced options tab is not available when the cdd search is selected.
129
The dialog will appear that allows specifying repeat parameters and the annotations table document to save the results into:
The dialogues status line displays approximate repeats number that will be found with the current settings. The Advanced tab provides additional repeats nding options:
130
The found repeats are saved and displayed as annotations to the DNA sequence:
131
The dialog parameters: Tandem preset specify the tandem repeats parameters with predened values by selecting the available preset:
Min period, Max period the minimum and maximum acceptable repeat length measured in base symbols. Region to process specify the region to search in the whole sequence, a custom region or the region of the current selection (if any). Save annotation(s) to specify the existing or new annotations table le to save the resulting annotations into. Annotation parameters you can change the default group name and annotation(s) name values of the resulting annotation(s).
132
Additional search options can be found in the Advanced tab: Algorithm the algorithm parameter allows to select the search algorithm. The default and a fast one is optimized sux array algorithm. Minimum tandem size the minimum tandem size sets the limit on minimum acceptable length of the tandem, i.e. the minimum total repeats length of the searched tandem. Minimum repeat count the minimum number of repeats of a searched tandem. Show overlapped tandems check if the plugin should search for the overlapped tandems, otherwise keep unchecked. Tandem Repeats Search Result An example of the search results for the micro-satellite preset:
133
Alternatively, select either the Actions Analyze Find restriction sites item in the main menu or the Analyze Find restriction sites item in the context menu. The Find restriction sites dialog appears:
You can see the list of restriction enzymes that can be used to search for restriction sites. The information about enzymes was obtained from the REBASE database. For each enzyme in the list a brief description is available (the accession ID in the database, the recognition sequence, etc.). If youre online you can get more detailed information about an enzyme selected by clicking the REBASE Info button.
134
135
10.12.6 Results
When at least one enzyme has been selected and the OK button has been pressed in the dialog, the autoannotating becomes enabled. In the Annotations editor the Restriction Sites annotations can be found in the Auto-annotations\enzyme group. The direct and complement cut site positions are visualized as triangles on an annotation in the Sequence details view :
136
On the Restriction Sites tab of the dialog you can see the name of the molecule, the list of restriction enzymes found during the restriction analysis that can cut the molecule and the list of enzymes selected to perform the digestion. To digest the sequence into fragments you should select at least one enzyme. To move an enzyme to the Selected enzymes list click on it in the Available enzymes list and press the Add button. Note that you can select several items in a list by holding the Ctrl key while clicking on the items. To select all available enzymes press the Add All button. To remove enzymes from the Selected enzymes list select them in the list and press the Remove button. To remove all items from the Selected enzymes list press the Clear Selection button. 10.13. Molecular Cloning in silico 137
As soon as the required restriction enzymes are selected press the OK button. The fragments will be saved as annotations. When the annotations are being saved the parameters from the Output tab of the dialog are applied.
If a region has been selected you can choose to create the fragment from this region. Otherwise you can either choose to create the fragment from the whole sequence or choose the Custom item and input the custom region. To add a 5 overhang to the direct strand check the Include 5 End Overhang check box and input the required nucleotides. To add a 5 overhang to the reverse strand in addition to the described steps select the Reversecomplement item in the same group box. Similarly, to add a 3 overhang check the Include 3 End Overhang check box, input the required overhang and select either the direct or the reverse-complement strand. On the Output tab of the dialog you can optionally modify the annotations output settings. Finally, press the OK button to create the fragment. The fragment will be saved as an annotation.
138
Available Fragments All the fragments available in the current project are shown in the Available fragments list. You can automatically create a fragment from a DNA molecule from the current UGENE project. Click the From Project button to do so. The Select Item dialog appears with the sequence objects available. Select a sequence and press the OK button. After that create a fragment in the appeared Create DNA Fragment dialog as described in the Creating Fragment paragraph. The fragment created from the sequence appears in the list of available fragments. Fragments of the New Molecule The next step is to add required fragments to the new molecule contents. To add fragments select them in the list of available fragments and click the Add button. To add all the fragments click the Add All button. Changing Fragments Order in the New Molecule To change the order of fragments in the new molecule select a fragment in the new molecule contents list and click either the Up or the Down button to move the fragment in the corresponding direction.
139
Removing Fragment from the New Molecule To remove a fragment from the new molecule select it in the new molecule contents list and click the Remove button. To remove all the fragments click the Clear All button. Editing Fragment Overhangs To edit a fragments overhangs select the fragment in the new molecule contents list and click the Edit Overhangs button. The Edit Molecule Fragment dialog appears:
Here you can select the type of each DNA end and even input a custom overhang. The changes youve made are shown in the Preview area of the dialog. To conrm the changes and close the dialog click the OK button. Reverse Complement a Fragment To reverse complement a fragment check the Inverted check box for the fragment in the new molecule contents list.
140
Other Constuction Options To save the fragments of the new molecule as annotations check the Annotate fragments in new molecule check box. To make all DNA ends blunt check the Force "blunt" and omit all overhangs check box. All overhangs would be cut in this case. Check the Make circular check box to make the new molecule circular. Output On the Output tab of the dialog you can select the le to save the new molecule to. The molecule is opened by default as soon as it is created. To modify this behavior uncheck the Open view for new molecule check box on the same tab. To save the molecule le to the hard disk immediately after it is created check the Save immediately check box. Otherwise it would be stored in memory until you save or remove it.
141
It supports the following options: Algorithm you can choose the preferred algorithm. Currently, GORIV and PsiPred algorithms are available. Range start / Range end select the sequence range for prediction. Results visual representation of the prediction results, for example:
142
Save as annotation select this button to save the results as annotations of the current protein sequence.
143
10.15 SITECON
SITECON is a program package for recognition of potential transcription factor binding sites basing on the data about conservative conformational and physicochemical properties revealed on the basis of the binding sites sets analysis. To cite SITECON use the following article: "Oshchepkov D.Y., Vityaev E.E., Grigorovich D.A., Ignatieva E.V., Khlebodarova T.M.SITECON:a tool for detecting conservative conformational and physicochemicalproperties in transcription factor binding site alignments and for siterecognition. // Nucleic Acids Res. 2004 Jul 1;32(Web Server issue):W208-12." UGENE version of SITECON provides a tool for recognition of potential binding sites for over 90 types of transcription factors. Also UGENE version of SITECON provides a tool for recognition of potential binding sites basing site alignment proposed by user. For the detailed method description see the original SITECON site. Data about used context-dependent conformational and physicochemical properties are available in the PROPERTY Database. To search transcription factor binding sites in a DNA sequence select the Analyze Search TFBS with SITECON... context menu item. In the appeared search dialog you must select a le with TFBS prole. The proles supplied with UGENE are placed in the $UGENE/data/sitecon_models folder. After the prole is loaded the threshold-lter is populated with values read from prole. You can use the lter to remove low-scoring regions from the result.
The regions found by SITECON algorithm can be saved as annotations to the DNA sequence in the Genbank format. 144 Chapter 10. Plugins
Every SITECON prole supplied with UGENE contains complete information about calibration settings provided to UGENE team by the author of SITECON. The original TFBS alignments used to calculate proles can be requested directly from the author of SITECON.
10.15. SITECON
145
First of all you need to specify the pattern to search for. The rest parameters are optional: Search in select either to search in the sequence or in its translation. Strand select the strand to search in: direct, complementary or both strands. Range species the region of the sequence that will be used to search for the pattern. By default, if a subsequence has been selected when the dialog has been opened, then the selected subsequence is searched for the pattern. Otherwise, the whole sequence is used. You can also input a custom range. Algorithm version version of the algorithm implementation. Non-classic versions produce the same results as classic but much faster. To use these optimizations your system must support these capabilities. Classic 2 SSE2 CUDA OPENCL Scoring matrix can be chosen from a bunch of matrices supplied with UGENE. To view a matrix selected click the View button. Gap open penalty for opening a gap. Gap extension penalty for extending a gap 146 Chapter 10. Plugins
Report results simple heuristic which allows to lter intersected hits. If it is set to none, the algorithm may report large set of almost identical results in the same region. Minimal score another simple heuristic which measures sequences similarity. It is more convenient than using some abstract scores. If set to 100%, the algorithm will search for exact substring match. The results of the search are saved as annotations. To set the annotations parameters (Annotation name, Group name, a le to save the annotation to) go to the Input and output tab of the dialog.
147
10.17 HMM2
The HMM2 plugin is a toolkit based on the Sean Eddys HMMER2 package. While working on this plugin we were guided by the following principles: Make the HMMER2 tools accessible to a wider user audience by providing graphical interface for all supported utilities for most of the platforms. Be compatible with the original HMMER2 package. Create the high-performance solution utilizing modern multi-core processors and SIMD instructions. The current version of UGENE provides user interface for three HMM2 tools: HMM build , HMM calibrate and HMM search. In the original program the corresponding commands are: hmmbuild, hmmcalibrate and hmmsearch. To access these tools select the Tools HMMER2 tools submenu of the program main menu:
We highly recommend reading the original HMMER2 documentation to learn how to use utilities provided by the plugin. Note: SSE2 algorithm is implemented by Leonid Konyaev, Novosibirsk State University. Use of the SSE2 optimized version of the HMM search algorithm with quad-core CPU gives >30x performance boost when compared with the original single-threaded algorithm (single sequence mode).
148
Note: The HMM build tool does not automatically calibrate a prole. Use the HMM calibrate tool to calibrate the prole.
10.17. HMM2
149
The search results are stored as sequence annotations in the Genbank le format. Warning: All HMM2 UGENE tools work only with les that contain a single HMM model.
150
10.18 HMM3
The HMM3 plugin is a toolkit based on the Sean Eddys HMMER3 package. While working on this plugin we were guided by the following principles: Make the HMMER3 tools accessible to a wider user audience by providing graphical interface for all supported utilities for most of the platforms. Be compatible with the original HMMER3 package. Create the high-performance solution utilizing modern multi-core processors. The current version of UGENE provides user interface for three HMM3 tools: HMM3 build , HMM3 search and Phmmer search. In the original program the corresponding commands are: hmmbuild, hmmsearch and phmmer. To access these tools select the Tools HMMER3 tools submenu of the program main menu:
We highly recommend reading the original HMMER3 documentation to learn how to use utilities provided by the plugin.
10.18. HMM3
151
The HMM3 conguration dialog provides an easy way to set appropriate search parameters. Here you can see eective weighting strategies options:
152
For example, reporting thresholds options can be congured using the dialog:
The search results are stored as sequence annotations in the Genbank le format.
Warning: The HMM3 search works only with les that contain a single HMM model.
10.18. HMM3
153
You can set options of the Phmmer search by choosing the needed dialog tab. Here you can see the e-value calibration options:
154
10.18. HMM3
155
10.19 uMUSCLE
UGENE contains graphical ports of the Robert C. Edgars MUSCLE tool for multiple alignment. Note: MUSCLE4 is not supported since UGENE version 1.7.2. The package is integrated completely, so there is no need in extra les for using it. It is possible to run several multiple alignment tasks in parallel, check the progress and cancel the running tasks safely. Note: The k-mer clustering part of the MUSCLE algorithm was optimized for multicore systems by Timur Tleukenov, Novosibirsk State Technical University.
The dialog contains the list of MUSCLE modes: MUSCLE default, Large alignment, Rene only.
Warning: By default UGENE does not rearrange sequence order in an alignment, but the original MUSCLE package does. To enable sequence rearrangement uncheck the Do not re-arrange sequences (-stable) option in the dialog.
156
One of the improvements to the original MUSCLE package is the ability to align only a part of the model. When the Column range item is selected the region of the specied columns is only passed to the MUSCLE alignment engine. The resulted alignment is inserted into the original one with gaps added or removed on the region boundaries. Note: To visually select the column range to align, make a selection in the alignment editor rst. Then invoke the MUSCLE plugin. Its column range boundary values will automatically match the given selection.
There are two gap columns inserted into the source prole, and two gap columns inserted into the added one. Therefore the proles columns kept intact and the alignments havent been changed. Note: Aligning a prole to the active alignment you will modify the original alignment le, since it will contain 2 proles after the operation is completed.
10.19. uMUSCLE
157
The original alignment is not modied, only columns with gap () character can be inserted. The second prole was considered as a set of sequences and therefore is modied. Note that if a le with another alignment is used as a source of unaligned sequences, the gap characters are removed and each input sequence is processed independently. This method is quite fast, for example an alignment of 3000 sequences (1000 bases each) to the existing prole takes about 5 minutes on the usual Core2Duo computer.
158
10.20 Bowtie
The Bowtie plugin is a graphical port of the Bowtie short read aligner. The Bowtie packages are integrated completely, so there is no need in extra les to use these tools. The Bowtie functionality is available from the Tools DNA Assembly submenu of the main menu.
10.20. Bowtie
159
There are the following parameters: Reference sequence DNA sequence to align short reads to. This parameter is required. Result le name File in SAM format to write the result of the alignment into. This parameter is required. Short reads each added short read is a small DNA sequence le. At least one read should be added. Note: Short reads length for Bowtie cant be more than 1024. You can also congure other parameters. They are the same as in the original Bowtie (you can read their description on the Bowtie manual page). To start short reads aligning select these les and press the Start button. The Bowtie plugin can use a prebuilt index le. The index is a set of 6 les with suxes .1.ebwt, .2.ebwt, .3.ebwt, .4.ebwt, .rev.1.ebwt, and .rev.2.ebwt. The index is created during the aligning. Also you can build it manually . To use it check the Prebuilt index check box.
There are the following parameters: Reference sequence DNA sequence to which short reads would be aligned to. This parameter is required. Index le name le to save index to. This parameter is required.
160
161
The following parameters are available: Reference sequence DNA sequence to align short reads to. This parameter is required. Result le name le in SAM format to write the result of the alignment into. This parameter is required. Short reads each added short read is a small DNA sequence le. At least one read should be added. Note: UGENE Genome Aligner has no limitation on short reads length. Mismatches allowed check this box to allow mismatches between the reference sequence and a short read. Select one of the following: Mismatches number to set the number of mismatched nucleotides allowed. This parameter can take values: 1, 2 and 3. Percentage of mismatches to set the number of mismatches in percents. Note, that in the last case the absolute number of mismatches can vary for dierent reads. This parameter can take values: 1 - 10 %. Index use an appropriate index to optimize the aligning. Save index to value is selected by default, specifying to save index built during the aligning to the le. If this option is selected, the value is required. Use prebuilt index from option is used to optimize the aligning. Select an appropriate index built for the reference sequence. If this option is selected, the value is required. Note: An index is built for a reference sequence during the aligning or manually. It depends on the short reads length. Read Building Index for more information. Align reverse complement reads use both the read and its reverse complement during the aligning. Stop searching on the rst successful match search only for the rst occurence of each read in the reference sequence. Omit reads with qualities lower than omit all reads with qualities lower than the specied value. Reads that have no qualities are not omited.
162
The parameters are the following: Reference sequence DNA sequence to which short reads would be aligned to. This parameter is required. Index le name le to save index to. This parameter is required. Optimize for read length optimize the index created for reads with length equal to this value. Warning: The index can be used to align reads that are longer than the specied value (however it would be less optimized), but the index cant be used with reads shorter than this value. Note: An index is created by default when you align short reads to a reference sequence. The length of the shortest read is used in this case as the value to the Optimize for read length parameter. Mismatches allowed check this box to allow mismatches between the reference sequence and a short read. Select one of the following: Mismatches number to set the number of mismatched nucleotides allowed. This parameter can take values: 1, 2 and 3. Percentage of mismatches to set the number of mismatches in percents. Note, that in the last case the absolute number of mismatches can vary for dierent reads. This parameter can take values: 1 - 10 %.
163
164
In the search dialog you must specify a le with PWM or PFM. You can do so by pressing the browse button [1] and selecting the le. Also you can use the special interface to choose a JASPAR matrix by pressing the Search JASPAR database button [2]. Alternative way to specify the position weight/frequency matrix is to create a specic one from an alignment or a le with several sequences with the build a new matrix tool. After the prole (the matrix) is loaded, you can adjust the threshold value [3]. The threshold sets the minimal identity score for a result to pass. The more the result score is, the more it is homologically related to the aligned region. By changing the threshold you can lter low- scoring results. If the loaded matrix is a position frequency matrix, you must also specify the algorithm to build the corresponding position weight matrix which will represent the transcription factor. There are four algorithms available.
Also you can add a selected matrix with the specied Minimal score and the Algorithm to the matrices list. To do it, select the matrix and other options and press the Add to queue button. The plugin will search with all matrices specied in the list. You can use the Save list... button to export the list of matrices to a *.csv le. Later the list can be loaded from the le using the Load list... button. The rest options are standard sequence search options: the strand and the sequence region where to search for matches. After specifying the necessary options press the Search button. The found results will appear in the dialog table. The corresponding results identity scores are in the Score column.
The regions found by the weight matrix algorithm can be saved as annotations to the DNA sequence in the Genbank format by pressing the Save as annotations button. After saving, the le with resulting annotations will be automatically added to the current project, and the annotations will be added to the original sequence. Note that in case of selecting JASPAR or UNIPROBE matrix, the resulting annotations will contain the given matrix properties.
165
166
Here the matrices are divided into categories and you can read detailed information of a matrix which is represented by its properties. It could help you to choose the matrix properly. Note: The matrices provided with UGENE are located in the $UGENE/data/position_weight_matrix folder.
The following parameters are available: Input le an alignment or a le with several sequences to build the matrix from. The parameter is mandatory. Output le the resulting matrix will be saved in this le. The parameter is mandatory. Statistic type denes the way in which the statistics will be collected. The Mononucleic option is basically good for small alignments, and the Dinucleic option must give more appropriate results for big alignments. Matrix type denes the type of the resulting matrix.
167
If the Frequency matrix option is selected then the frequency matrix will be created and saved into the resulting le. If the Weight matrix option is selected then the intermediate frequency matrix will be created and then transformed into a weight matrix on basis of the selected Weight algorithm. Then the weight matrix will be saved into the resulting le. For some input les the colored Alignment Logo appears at the bottom of the dialog. It gives the representation of the selected alignment.
Note: The Alignment logo appears when: The input le format is *.pfm, *.aln or it is a le with several sequences; The size of the input le is small enough. To start the operation, press the Start button. The matrix will be created and saved. If the Build weight or frequency matrix dialog was invoked from the Weight matrix search dialog, then the matrix also will be chosen as the current prole.
168
10.23 Dotplot
The Dotplot plugin provides a tool to build dotplots for DNA or RNA sequences. This allows to compare these sequences graphically. Using a dotplot graphic, you can easily identify such dierences between sequences as mutations, invertions, insertions, deletions and low-complexity regions. Also the plugin provides advanced features: comparing multiple dotplots, navigation in a dotplot, dotplots synchronization, saving and loading a dotplot, etc. An example of a dotplot view:
Note: The Dotplot plugin uses the Repeat Finder plugin to build a dotplot, make sure you have the Repeat Finder plugin installed. The Dotplot features are described in more details below.
10.23. Dotplot
169
Here you should specify the File with rst sequence. Also you should either check the Compare sequence against itself option or select the File with second sequence. Optionally you can select to Join all sequences found in the le (for the rst and/or for the second le). If you select to join the sequences you can also select the Gap size. The gap of the specied size will be inserted between the joined sequences. After you press the Next button, the dialog to congure the dotplot parameters will appear:
The following parameters are available: X axis sequence the sequence for the X dotplot axis. Y axis sequence the sequence for the Y dotplot axis. If there are several sequences in the specied (the rst or the second) le and you havent selected to join the sequences in the previous dialog, then you can select a sequence in these elds. If you have selected to Join all sequences found in the le, then you cant select a separate sequence from the le, the joined Sequence can be selected instead. Search direct repeats check this option to search for direct repeats in the specied sequences. You can also select the color with which the repeats will be displayed in the picture. The default button sets the default color. Search inverted repeats check this option to search for inverted repeats in the specied sequences.
170
You can also select the color with which the repeats will be displayed in the picture. The default button sets the default color. Custom algorithm optionally you can select an algorithm to calculate the repeats: Auto Sux index Diagonals Note: The specied algorithm is provided to the Repeat Finder plugin as an input parameter. In most cases the Auto value is appropriate. Minimum repeat length allows to draw only such matches between the sequences that are continuous and long enough. For example if it equals to 3bp, then only repeats will be found that contain 3 and more base symbols. Press the 1k button to automatically adjust the Minimum repeat length value. Such value will be set, that there will be about 1000 repeats found. Repeats identity species the percents of the repeats identity. Press the 100 button to set the 100% identity. After the parameters are set, press the OK button. The dotplot will appear in the Sequence view :
Each dot on the plot corresponds to a matched base symbol at the "x" position of the horizontal sequence and the "y" position of the vertical sequence. Visible diagonal lines indicate matches between sequences in the given particular region. See also: Interpreting Dotplot: Identifying Matches, Mutations, Ivertions, etc. Building Dotplot for Currently Opened Sequence
Hold the middle mouse button and move the mouse cursor over the zoomed region of the doplot. Click on the desired region of the minimap in the right bottom corner. Activate the Scroll tool, hold the left mouse button and move the mouse cursor over the zoomed region:
172
hold down the left mouse button and drag the mouse cursor over the dotplot. When you select a region on a dotplot the corresponding region is also selected in other Sequence view areas (Sequence details view , Sequence zoom view , etc.). The opposite is true as well: if you select a region in a Sequence view area, the corresponding region is also selected in the dotplot view. To zoom to the region selected click the Zoom in on the left.
To deselect the repeat either click on other repeat or hold Ctrl and click somewhere on the dotplot.
10.23. Dotplot
173
2. Frame shifts a. Mutations Mutations are distinctions between sequences. On the graphic they are represented by gaps in diagonal lines. They interrupt matches. b. Insertions Insertions are parts of one sequence that are missed in the another, while the surrounding parts match. In other words, an insertion is a subsequence that was inserted into a sequence. Graphically, insertions are represented by gaps which lie only on one axis. A little shift towards the other axis indicates a mutation involved. c. Deletions A deletion is a subsequence that was deleted from a sequence. A deletion from sequence A found in sequence B can be considered as an insertion into sequence B and contained in sequence A.
3. Inverted repeats The Dotplot plugin allows to search for inverted repeats as well. Inverted repeats are shown contrary to the direct repeats. Use the Search direct repeats and Search inverted repeats options of the Dotplot parameters dialog to select which repeats to draw (the dialog is described here).
174
4. Low-complexity regions A low-complexity region is a region produced by redundancy in a particular part of the sequence. It is represented on a plot as a rectangular area lled with the matches.
The parameters dialog will be re-opened. See description of the available parameters here.
10.23. Dotplot
175
The Save Dotplot dialog will appear. A dotplot is saved in a le with the *.dpt extension. Later the dotplot can be loaded using the Dotplot Save/Load Load context menu item.
176
10.23. Dotplot
177
10.24 Primer3
The Primer3 plugin is a port of the Primer3 tool. It is intended to pick primers from a DNA sequence. To use the Primer3 , open a DNA sequence and select the Analyze Primer3... context menu item. The dialog will appear:
178
10.25.1 BLAST/BLAST+
The Basic Local Alignment Search Tool (BLAST) nds regions of local similarity between sequences. The program compares nucleotide or protein sequences to sequence databases and calculates the statistical signicance of matches. BLAST can be used to infer functional and evolutionary relationships between sequences as well as help identify members of gene families. BLAST+ is a new version of the BLAST package from the NCBI. From UGENE you can use the following tools of the old BLAST package: blastall the old program developed and distributed by the NCBI for running BLAST searches. formatdb formats protein or nucleotide source databases before these databases can be searched by blastall. And the following tools of the new BLAST+ package: blastn searches a nucleotide database using a nucleotide query. blastp searches a protein database using a protein query. blastx searches a protein database using a translated nucleotide query. 10.25. External Tools 179
tblastn compares a protein query against a translated nucleotide database (the all six reading frames). tblastx translates the query nucleotide sequence in all six possible frames and compares it against the six-frame translations of a nucleotide sequence database. makeblastdb formats protein or nucleotide source databases before these databases can be searched by other BLAST+ tools. BLAST home page: http://blast.ncbi.nlm.nih.gov/Blast.cgi?CMD=Web&PAGE_TYPE=BlastHome To make BLAST (or BLAST+) tools available from UGENE: 1. Install the required verion of BLAST (or BLAST+) on your system. 2. Set the paths to the executables, you are going to use, on the External tools tab of UGENE Application Settings dialog. After youve nished this conguration you can access the tools from the Tools BLAST submenu of the main menu. Creating Database To format a BLAST database do the following: If youre using BLAST open Tools BLAST FormatDB. If youre using BLAST+ open Tools BLAST BLAST+ make DB. The Format database dialog appears:
180
Here you must select the input les. If all the les you want to use are located in one directory, you can simply select the directory with the les. By default only the les are taken into account with *.fa and *.fasta extensions. You can change this by specifying either Include les lter or Exclude les lter. You can choose either protein or nucleotide type of the les. Then you must select the path to save the database le and specify a Base name for BLAST les and a Title for database le. Making Request to Database To make a request to a local BLAST database do the following: If youre using BLAST open Tools BLAST BLAST Search. If youre using BLAST+ open Open Tools BLAST BLAST+ Search. If there is a sequence opened you can also initiate the request to a local BLAST database from the Sequence view: If youre using BLAST select the Analyze Query with BLAST item in the context menu or in the Actions main menu. If youre using BLAST+ select the Analyze Query with BLAST+ item in the context menu or in the Actions main menu. The Request to local BLAST database dialog appears:
181
The dialog is very similar to the dialog described in the Remote BLAST chapter, except the following parameters: Select input le this parameter is only presented if the dialog has been opened from the Tools main menu. Here you must input a query sequence le that would be used to search the BLAST database. If the dialog has been opened e.g. using the Sequence view context menu, then the currently active sequence is used as a query sequence. Search type here you should select the tool you would like to use. If the query sequence is a nucleotide sequence then blastn, blastx and tblastx items are available. For a protein sequence the items are blastp and tblastn. Select database path path to the database les. Base name for BLAST DB les base name for the BLAST database les. Number of CPUs being used number of processors to use. To learn about other parameters, please, refer to the Remote BLAST chapter.
10.25.2 ClustalW
Clustal is a widely used multiple sequence alignment program. It is used for both nucleotide and protein sequences. ClustalW is a command-line version of the program. Clustal home page: http://www.clustal.org If you are using Windows OS, there are no additional conguration steps required, as ClustalW executable le is included to the UGENE distribution package. Otherwise: 1. Install the Clustal program on your system. 2. Set the path to the ClustalW executable on the External tools tab of UGENE Application Settings dialog. Now you are able to use Clustal from UGENE. Open a multiple sequence alignment le and select the Align with ClustalW item in the context menu or in the Actions main menu. The Align with ClustalW dialog appears (see below), where you can adjust the following parameters: Gap opening penalty cost of opening up a new gap in the alignment. Increasing this value will make gaps less frequent. Gap extension penalty cost of every item in a gap. Increasing this value will make gaps shorter. Weight matrix species a single weight matrix for nucleotide sequences or series of matrices for protein sequences. For nucleotide sequences the weight matrix selected denes the scores assigned to matches and mismatches (including IUB ambiguity codes), it can take values: IUB default scoring matrix used by BESTFIT for the comparison of nucleic acid sequences. Xs and Ns are treated as matches to any IUB ambiguity symbol. All matches score 1.9; all mismatches for IUB symbols score 0. CLUSTALW previous system used by ClustalW, in which matches score 1.0 and mismatches score 0. All matches for IUB symbols also score 0. For protein sequences it describes the similarity of each amino acid to each other. The following values are available:
182
BLOSUM BLOcks of Amino Acid SUbstitution Matrices rst introduced in a paper by Heniko and Heniko. These matrices appear to be the best available for carrying out data base similarity (homology searches). PAM Point Accepted Mutation matrices introduced by Margaret Dayho. These have been extremely widely used since the late 70s. GONNET these matrices were derived using almost the same procedure as the Dayho one (above) but are much more up to date and are based on a far larger data set. They appear to be more sensitive than the Dayho series. ID identity matrix which gives a score of 1.0 to two identical amino acids and a score of zero otherwise. Iteration type species the iteration type to use. During the iteration step each sequence is removed in turn and realigned. It is kept if the resulting alignment is better than the one has been made before. This process is repeated until the score converges or until the maximum number of iterations is reached. Available values are: NONE species not to use iterations. TREE species to iterate at each step of the progressive alignment. ALIGNMENT species to iterate on the nal alignment. Max iterations maximum number of iterations.
The following parameters are only available for protein sequences: Gap separation distance tries to decrease the chances of gaps being too close to each other. Gaps that are less than this distance apart are penalized more than other gaps. This does not prevent close gaps; it makes them less frequent, promoting a block-like appearance of the alignment. Hydrophilic gaps o increases the chances of a gap within a run of hydrophilic amino acids. No end gap separation penalty treats end gaps just like internal gaps to avoid gaps that are too close. 10.25. External Tools 183
Residue-specic gaps o amino acid specic gap penalties that reduce or increase the gap opening penalties at each position in the alignment or sequence. For example, positions that are rich in glycine are more likely to have an adjacent gap than positions that are rich in valine.
10.25.3 MAFFT
Originally, MAFFT is a multiple sequence alignment program for unix-like operating systems. However, currently it is available for Mac OS X, Linux and Windows. It is used for both nucleotide and protein sequences. MAFFT home page: http://mat.cbrc.jp/alignment/software To make MAFFT available from UGENE: 1. Install the MAFFT program on your system. 2. Set the path to the MAFFT executable on the External tools tab of UGENE Application Settings dialog. For example, on Windows you need to specify the path to the mafft.bat le. To use MAFFT open a multiple sequence alignment le and select the Align with MAFFT item in the context menu or in the Actions main menu. The following dialog appears:
The following parameters are available: Gap opening penalty Gap opening penalty at group-to-group alignment. Oset (works like gap extension penalty) oset value, which works like gap extension penalty, for group-to-group alignment. Maximum number of iterative rene species the number of cycles of iterative renement to perform.
10.25.4 T-Coffee
T-Coee is a multiple sequence alignment package. T-Coee home page: http://www.tcoee.org/Projects_home_page/t_coee_home_page.html Note: Currently, T-Coee tool is only available on UNIX and UNIX-like platforms (Linux, cygwin and Mac OS X). To make T-Coee available from UGENE: 1. Install the T-Coee program on your system. 2. Set the path to the T-Coee executable on the External tools tab of UGENE Application Settings dialog.
184
To use T-Coee open a multiple sequence alignment le and select the Align with T-Coee item in the context menu or in the Actions main menu. The following dialog appears:
The following parameters are available: Gap opening penalty indicates the penalty applied for opening a gap. The penalty must be negative. Gap extension penalty indicates the penalty applied for extending a gap. Number of iterations species the number of iterations.
185
Alternatively, you can create / edit a schema using a text editor. When the schema has been created and all its parameters have been set you can run it for a nucleotide sequence. The results are saved as a set of annotations to the specied le in the Genbank format. To learn more about the Query Designer read the Query Designer Manual (follow the link on the UGENE documentation page).
186
Here: task_name task to execute, it can be one of the predened tasks or a task you have created . task_parameter parameter of the specied task. Some parameters of a task are required, like in and out parameters of some tasks. option one of the CLI options. See the example below:
ugene align --in=COI.aln -out result.aln -log-level-details
--task=<task_name> [<task_parameter>=value ...] Species the task to run. A user-dened UGENE workow schema can be used as a task name. For example:
ugene --task=align --in=COI.aln -out result.aln ugene --task=C:\myschema.uwl --in=COI.aln --out=res.aln
187
--log-no-task-progress A task progress is shown by default when a task is running. This option species not to show the progress. --log-level="[<category1>=]<level1> [, ...]" Sets the log level per category. If a category is not specied, the log level is applied to all categories. The following categories are available: "Algorithms" "Console" "Core Services" "Input/Output" "Performance" "Remote Service" "Scripts" "Tasks". The following log levels are available: TRACE, DETAILS, INFO, ERROR or NONE. By default, loglevel=ERROR. For example:
ugene --log-level=NONE ugene --log-level="Tasks=DETAILS, Console=DETAILS"
--log-format="<format_string>" Species the format of a log line. Use the following notations: L - level, C - category, YYYY or YY - year, MM - month, dd - day, hh hour, mm - minutes, ss - seconds, zzz - milliseconds. By default, logformat="[L][hh:mm]". --license Shows license information. --lang=language_code Species the language to use (e.g. for the log output). The following values are available: CS (Czech) EN (English) RU (Russian) --log-color-output If log output is enabled, this option make it colored: ERROR messages are displayed in red, DETAILS messages are displayed in green, TRACE messages are displayed in blue.
188
189
190
191
Example:
ugene find-orfs --in=human_T1.fa --out=result.gb --require-init-codon=false
192
pam250 etc. The matrices available are stored in the $UGENE\data\weight_matrix directory. filter results ltering strategy. [String, Optional, Default: "lter-intersections"] The following values are available: lter-intersections none Example:
ugene find-sw --in=human_T1.fa --out=sw.gb --ptrn=TGCT --filter=none
193
p type of the BLAST search. [String, Optional, Default: "blastn"] The following values are available: blastn blastp blastx tblastn tblastx e expectation value threshold. [Number, Optional, Default: 10] Example:
ugene local-blast --in=input.fa --dbpath=. --dbname=mydb --out=output.gb
194
195
Example:
ugene query --in=input.fa --out=result.gb --schema=RepeatsWithORF.uql
seed seed for pseudo-random number generator. [Number, Optional] seedlen number of bases on the high-quality end of the read to which the n ceiling applies. The lowest permitted setting is 5 and the default is 28. Bowtie is faster for larger values of seedlen. [Number, Optional] tryhard nds valid alignments when they exist, including paired-end alignments. [Boolean, Optional] chunkmbs number of megabytes a certain thread is given to store path descriptors in best mode. [Number, Optional] best guarantees that reported singleton alignments are "best" in terms of stratum (i.e. number of mismatches, or mismatches in the seed in the case of n mode) and in terms of the quality values at the mismatched position(s). Example:
ugene bowtie --reads=r1.fa;r2.fa;r3.fa --ebwt=refindex --out=result.aln
197
score score based ltering which is an alternative to e-value ltering to exclude low-probability hits from the result. [Number, Optional, Default: -1000000000] Example:
ugene hmm2-search --seq=CBS_seq.fa --hmm=CBS.hmm --out=CBS_hmm.gb
198
Parameters: toolpath path to the MAFFT executable. By default, the path specied in the Application Settings is applied. [String, Optional, Default: "default"] tmpdir directory for temporary les. [String, Optional] in semicolon-separated list of input les. [String, Required] out output le. [String, Required] format format of the output le. [String, Required] op penalty for opening a gap. [Number, Optional] ep penalty for extending a gap. [Number, Optional] maxiterate maximum number of cycles of iterative renement. [Number, Optional] Example:
ugene align-mafft --in=COI.aln --out=COI_aligned.aln
199
Example:
ugene pfm-build --in=COI.aln --out=result.pfm
200
201
Example:
ugene pwm-search --seq=input.fa --matrix=Aro80.pwm;Aft1.pwm --out=res.gb
max-err2 setting for ltering results, maximum value of Error type II. [Number, Optional, Default: 0.001] strand strands to search in. [Number, Optional, Default: 0] The following values are available: 0 (both strands) 1 (direct strand) -1 (complement strand) Example:
ugene sitecon-search --in=input.fa --inmodel=profile.sitecon --out=res.gb
203
12 APPENDICES
12.1 Appendix A. Supported File Formats
Note: UGENE is able to read and write les compressed with Unix/Linux gzip utility. You dont have to unpack the les.
+ + +
+ +
EBWT EMBL
+ +
+ -
FASTA
FASTQ
Genbank
204
GFF
*.g
The Gene Finding Format (GFF) format is used to store features and annotations. See also: Sequence View A le format to store HMM proles. See also: HMM2 , HMM3 ASN.1 format used by the Molecular Modeling Database (MMDB). See also: 3D Structure Viewer A multiple sequence alignments le format. See also: Alignment Editor A tree le format. See also: create_a_phy_tree, Phylogenetic Tree Viewer A multiple alignment and phylogenetic trees le format. See also: Alignment Editor , create_a_phy_tree, Phylogenetic Tree Viewer The Protein Data Bank (PDB) format allows to view the 3D structure of the sequence. See also: 3D Structure Viewer A sequence le format used by pDRAW32 software. See also: Sequence View A le format for a position frequency matrix. See also: Weight Matrix A le format for a position weight matrix. See also: Weight Matrix A raw sequence format. See also: Sequence View The Sequence Alignment/Map (SAM) format is a generic alignment format for storing read alignments against reference sequences. See also: Alignment Editor , Bowtie, UGENE Genome Aligner It is a Standard Chromatogram Format. See also: Chromatogram Viewer A le format to store TFBS prole. See also: SITECON A multiple sequence alignments le format. See also: Alignment Editor
HMM MMDB
*.hmm *.prt
+ +
+ -
MSF Newick
+ +
+ +
Nexus
*.nex *.nxs
PDB
*.pdb
pDRAW32
*.pdw
+ + + +
+ + + +
+ + +
205
Swiss-Prot
*.txt *.sw
An annotated protein sequence in format of the UniProtKB/Swiss-Prot database. See also: Sequence View
206
Reads
+ +
+ +
UGENE Workow Designer schema Workow ement command tool elfor line
*.etc
207