|
@@ -2,88 +2,90 @@
|
|
|
|
|
|
/**
|
|
/**
|
|
* @file
|
|
* @file
|
|
- * Contains all functions for the protein blast module
|
|
|
|
|
|
+ * Contains all functions for the Protein BLAST
|
|
*/
|
|
*/
|
|
|
|
|
|
|
|
|
|
/**
|
|
/**
|
|
- * Form constructor for the protein blast form.
|
|
|
|
|
|
+ * Protein BLAST Submission Form
|
|
*
|
|
*
|
|
* @see blast_protein_form_validate()
|
|
* @see blast_protein_form_validate()
|
|
* @see blast_protein_form_submit()
|
|
* @see blast_protein_form_submit()
|
|
- *
|
|
|
|
*/
|
|
*/
|
|
function blast_protein_form($form, &$form_state) {
|
|
function blast_protein_form($form, &$form_state) {
|
|
-
|
|
|
|
-// Protein BLAST
|
|
|
|
|
|
|
|
-$form['pBLAST'] = array(
|
|
|
|
|
|
+ // PROTEIN QUERY
|
|
|
|
+ //.........................
|
|
|
|
+
|
|
|
|
+ $form['query'] = array(
|
|
'#type' => 'fieldset',
|
|
'#type' => 'fieldset',
|
|
'#title' => t('Enter Query Sequence'),
|
|
'#title' => t('Enter Query Sequence'),
|
|
'#description' => t('Enter one or more queries in the top text box or use the browse button to upload a file from your local disk. The file may contain a single sequence or a list of sequences. In both cases, the data must be in FASTA format. <a href="http://www.ncbi.nlm.nih.gov/BLAST/blastcgihelp.shtml" target="_blank">More information..</a> '),
|
|
'#description' => t('Enter one or more queries in the top text box or use the browse button to upload a file from your local disk. The file may contain a single sequence or a list of sequences. In both cases, the data must be in FASTA format. <a href="http://www.ncbi.nlm.nih.gov/BLAST/blastcgihelp.shtml" target="_blank">More information..</a> '),
|
|
- '#collapsible' => TRUE,
|
|
|
|
- '#collapsed' => FALSE,
|
|
|
|
|
|
+ '#collapsible' => TRUE,
|
|
|
|
+ '#collapsed' => FALSE,
|
|
);
|
|
);
|
|
|
|
|
|
-$form['pBLAST']['FASTA'] = array(
|
|
|
|
|
|
+ $form['query']['FASTA'] = array(
|
|
'#type' => 'textarea',
|
|
'#type' => 'textarea',
|
|
'#title' => t('Enter FASTA sequence(s)'),
|
|
'#title' => t('Enter FASTA sequence(s)'),
|
|
'#description'=>t('Enter query sequence(s) in the text area.'),
|
|
'#description'=>t('Enter query sequence(s) in the text area.'),
|
|
);
|
|
);
|
|
|
|
|
|
-// Upload a file as an alternative to enter a query sequence
|
|
|
|
-$form['#attributes']['enctype'] = 'multipart/form-data';
|
|
|
|
-$form['pBLAST']['UPLOAD'] = array(
|
|
|
|
|
|
+ // Upload a file as an alternative to enter a query sequence
|
|
|
|
+ $form['#attributes']['enctype'] = 'multipart/form-data';
|
|
|
|
+ $form['query']['UPLOAD'] = array(
|
|
'#prefix' => 'Or upload your query files: ',
|
|
'#prefix' => 'Or upload your query files: ',
|
|
'#type' => 'file',
|
|
'#type' => 'file',
|
|
'#description' => t('The file should be a plain-text FASTA file and not a .doc, .docx, etc. It cannot be greater than 10 Mb in size.'),
|
|
'#description' => t('The file should be a plain-text FASTA file and not a .doc, .docx, etc. It cannot be greater than 10 Mb in size.'),
|
|
);
|
|
);
|
|
|
|
|
|
|
|
|
|
-//DATABASE//
|
|
|
|
-
|
|
|
|
-$form['DB'] = array(
|
|
|
|
- '#type' => 'fieldset',
|
|
|
|
- '#title' => t('Choose Search Set'),
|
|
|
|
- '#description' => t('Choose from one of the protein BLAST databases listed below. You can also use the browse button to upload a file from your local disk. The file may contain a single sequence or a list of sequences. '),
|
|
|
|
- '#collapsible' => TRUE,
|
|
|
|
- '#collapsed' => FALSE,
|
|
|
|
-);
|
|
|
|
-
|
|
|
|
-$options = DB_options();
|
|
|
|
-$form['DB']['SELECT_DB'] = array(
|
|
|
|
- '#type' => 'select',
|
|
|
|
- '#title' => t('Protein BLAST Databases:'),
|
|
|
|
- '#options' => $options,
|
|
|
|
- '#default_value' => t('Select a database'),
|
|
|
|
-);
|
|
|
|
-
|
|
|
|
-// Upload a file as an alternative to enter a query sequence
|
|
|
|
-$form['#attributes']['enctype'] = 'multipart/form-data';
|
|
|
|
-$form['DB']['DBUPLOAD'] = array(
|
|
|
|
|
|
+ // BLAST DATABASE
|
|
|
|
+ //.........................
|
|
|
|
+
|
|
|
|
+ $form['DB'] = array(
|
|
|
|
+ '#type' => 'fieldset',
|
|
|
|
+ '#title' => t('Choose Search Set'),
|
|
|
|
+ '#description' => t('Choose from one of the protein BLAST databases listed below. You can also use the browse button to upload a file from your local disk. The file may contain a single sequence or a list of sequences. '),
|
|
|
|
+ '#collapsible' => TRUE,
|
|
|
|
+ '#collapsed' => FALSE,
|
|
|
|
+ );
|
|
|
|
+
|
|
|
|
+ $options = DB_options();
|
|
|
|
+ $form['DB']['SELECT_DB'] = array(
|
|
|
|
+ '#type' => 'select',
|
|
|
|
+ '#title' => t('Protein BLAST Databases:'),
|
|
|
|
+ '#options' => $options,
|
|
|
|
+ '#default_value' => t('Select a database'),
|
|
|
|
+ );
|
|
|
|
+
|
|
|
|
+ // Upload a file as an alternative to enter a query sequence
|
|
|
|
+ $form['#attributes']['enctype'] = 'multipart/form-data';
|
|
|
|
+ $form['DB']['DBUPLOAD'] = array(
|
|
'#prefix' => 'Or upload your own dataset: ',
|
|
'#prefix' => 'Or upload your own dataset: ',
|
|
'#type' => 'file',
|
|
'#type' => 'file',
|
|
'#description' => t('The file should be a plain-text FASTA file and not a .doc, .docx, etc. It cannot be greater than 10 Mb in size.'),
|
|
'#description' => t('The file should be a plain-text FASTA file and not a .doc, .docx, etc. It cannot be greater than 10 Mb in size.'),
|
|
-);
|
|
|
|
-
|
|
|
|
-//Algorithm parameters
|
|
|
|
|
|
+ );
|
|
|
|
+
|
|
|
|
+ // ALGORITHM PARAMETERS
|
|
|
|
+ //.........................
|
|
|
|
|
|
$form['ALG'] = array(
|
|
$form['ALG'] = array(
|
|
'#type' => 'fieldset',
|
|
'#type' => 'fieldset',
|
|
'#title' => t('Algorithm parameters'),
|
|
'#title' => t('Algorithm parameters'),
|
|
- '#collapsible' => TRUE,
|
|
|
|
- '#collapsed' => TRUE,
|
|
|
|
-);
|
|
|
|
|
|
+ '#collapsible' => TRUE,
|
|
|
|
+ '#collapsed' => TRUE,
|
|
|
|
+ );
|
|
|
|
|
|
-//General parameters
|
|
|
|
|
|
+ //General parameters
|
|
|
|
|
|
- $form['ALG']['GParam'] = array(
|
|
|
|
|
|
+ $form['ALG']['GParam'] = array(
|
|
'#type' => 'fieldset',
|
|
'#type' => 'fieldset',
|
|
'#title' => t('General parameters'),
|
|
'#title' => t('General parameters'),
|
|
'#collapsible' => FALSE,
|
|
'#collapsible' => FALSE,
|
|
- );
|
|
|
|
-
|
|
|
|
-$form['ALG']['GParam']['maxTarget'] = array(
|
|
|
|
|
|
+ );
|
|
|
|
+
|
|
|
|
+ $form['ALG']['GParam']['maxTarget'] = array(
|
|
'#type' => 'select',
|
|
'#type' => 'select',
|
|
'#title' => t('Max target sequences:'),
|
|
'#title' => t('Max target sequences:'),
|
|
'#options' => array(
|
|
'#options' => array(
|
|
@@ -99,24 +101,24 @@ $form['ALG']['GParam']['maxTarget'] = array(
|
|
),
|
|
),
|
|
'#default_value' => 2,
|
|
'#default_value' => 2,
|
|
'#description' => t('Select the maximum number of aligned sequences to display'),
|
|
'#description' => t('Select the maximum number of aligned sequences to display'),
|
|
-);
|
|
|
|
-
|
|
|
|
-$form['ALG']['GParam']['shortQueries'] = array(
|
|
|
|
- '#type' => 'checkbox',
|
|
|
|
|
|
+ );
|
|
|
|
+
|
|
|
|
+ $form['ALG']['GParam']['shortQueries'] = array(
|
|
|
|
+ '#type' => 'checkbox',
|
|
'#title' => t('Automatically adjust parameters for short input sequences'),
|
|
'#title' => t('Automatically adjust parameters for short input sequences'),
|
|
'#default_value' => TRUE,
|
|
'#default_value' => TRUE,
|
|
-);
|
|
|
|
-
|
|
|
|
-$form['ALG']['GParam']['eVal'] = array(
|
|
|
|
- '#type' => 'textfield',
|
|
|
|
- '#title' => t('e-value(Expect threshold)'),
|
|
|
|
- '#default_value' => 10,
|
|
|
|
- '#size' => 12,
|
|
|
|
- '#maxlength' => 20,
|
|
|
|
|
|
+ );
|
|
|
|
+
|
|
|
|
+ $form['ALG']['GParam']['eVal'] = array(
|
|
|
|
+ '#type' => 'textfield',
|
|
|
|
+ '#title' => t('e-value(Expect threshold)'),
|
|
|
|
+ '#default_value' => 10,
|
|
|
|
+ '#size' => 12,
|
|
|
|
+ '#maxlength' => 20,
|
|
'#description' => t('Expected number of chance matches in a random model.'),
|
|
'#description' => t('Expected number of chance matches in a random model.'),
|
|
-);
|
|
|
|
-
|
|
|
|
-$form['ALG']['GParam']['wordSize'] = array(
|
|
|
|
|
|
+ );
|
|
|
|
+
|
|
|
|
+ $form['ALG']['GParam']['wordSize'] = array(
|
|
'#type' => 'select',
|
|
'#type' => 'select',
|
|
'#title' => t('Word size:'),
|
|
'#title' => t('Word size:'),
|
|
'#options' => array(
|
|
'#options' => array(
|
|
@@ -125,71 +127,70 @@ $form['ALG']['GParam']['wordSize'] = array(
|
|
),
|
|
),
|
|
'#default_value' => 1,
|
|
'#default_value' => 1,
|
|
'#description' => t('The length of the seed that initiates an alignment'),
|
|
'#description' => t('The length of the seed that initiates an alignment'),
|
|
-);
|
|
|
|
-
|
|
|
|
-$form['ALG']['GParam']['qRange'] = array(
|
|
|
|
- '#type' => 'textfield',
|
|
|
|
- '#title' => t('Max matches in a query range'),
|
|
|
|
- '#default_value' => 0,
|
|
|
|
- '#size' => 12,
|
|
|
|
- '#maxlength' => 20,
|
|
|
|
|
|
+ );
|
|
|
|
+
|
|
|
|
+ $form['ALG']['GParam']['qRange'] = array(
|
|
|
|
+ '#type' => 'textfield',
|
|
|
|
+ '#title' => t('Max matches in a query range'),
|
|
|
|
+ '#default_value' => 0,
|
|
|
|
+ '#size' => 12,
|
|
|
|
+ '#maxlength' => 20,
|
|
'#description' => t('Limit the number of matches to a query range. This option is useful if many strong matches to one part of a query may prevent BLAST from presenting weaker matches to another part of the query.'),
|
|
'#description' => t('Limit the number of matches to a query range. This option is useful if many strong matches to one part of a query may prevent BLAST from presenting weaker matches to another part of the query.'),
|
|
-);
|
|
|
|
-
|
|
|
|
-// Scoring parameters
|
|
|
|
|
|
+ );
|
|
|
|
|
|
-$form['ALG']['SParam'] = array(
|
|
|
|
|
|
+ // Scoring parameters
|
|
|
|
+
|
|
|
|
+ $form['ALG']['SParam'] = array(
|
|
'#type' => 'fieldset',
|
|
'#type' => 'fieldset',
|
|
'#title' => t('Scoring parameters'),
|
|
'#title' => t('Scoring parameters'),
|
|
- '#collapsible' => FALSE,
|
|
|
|
-);
|
|
|
|
-
|
|
|
|
-$options_first = _ajax_example_get_first_dropdown_options();
|
|
|
|
-$selected = isset($form_state['values']['MATRIX'] ) ? $form_state['values']['MATRIX'] : key($options_first);
|
|
|
|
-
|
|
|
|
-$form['ALG']['SParam']['MATRIX'] = array(
|
|
|
|
- '#type' => 'select',
|
|
|
|
- '#title' => 'Matrix',
|
|
|
|
- '#options' => $options_first,
|
|
|
|
- '#default_value' => $selected,
|
|
|
|
- '#description' => t('Assigns a score for aligning pairs of residues, and determines overall alignment score..'),
|
|
|
|
- '#ajax' => array(
|
|
|
|
|
|
+ '#collapsible' => FALSE,
|
|
|
|
+ );
|
|
|
|
+
|
|
|
|
+ $options_first = _ajax_example_get_first_dropdown_options();
|
|
|
|
+ $selected = isset($form_state['values']['MATRIX'] ) ? $form_state['values']['MATRIX'] : key($options_first);
|
|
|
|
+
|
|
|
|
+ $form['ALG']['SParam']['MATRIX'] = array(
|
|
|
|
+ '#type' => 'select',
|
|
|
|
+ '#title' => 'Matrix',
|
|
|
|
+ '#options' => $options_first,
|
|
|
|
+ '#default_value' => $selected,
|
|
|
|
+ '#description' => t('Assigns a score for aligning pairs of residues, and determines overall alignment score..'),
|
|
|
|
+ '#ajax' => array(
|
|
'callback' => 'ajax_example_dependent_dropdown_callback',
|
|
'callback' => 'ajax_example_dependent_dropdown_callback',
|
|
'wrapper' => 'dropdown-second-replace',
|
|
'wrapper' => 'dropdown-second-replace',
|
|
- ),
|
|
|
|
-);
|
|
|
|
-
|
|
|
|
-$form['ALG']['SParam']['gapCost'] = array(
|
|
|
|
- '#type' => 'select',
|
|
|
|
- '#title' => t('Gap Costs:'),
|
|
|
|
- '#prefix' => '<div id="dropdown-second-replace">',
|
|
|
|
- '#suffix' => '</div>',
|
|
|
|
- '#options' => _ajax_example_get_second_dropdown_options($selected),
|
|
|
|
- '#default_value' => 2,
|
|
|
|
- '#description' => t('Cost to create and extend a gap in an alignment.'),
|
|
|
|
-);
|
|
|
|
-
|
|
|
|
-$form['ALG']['SParam']['M&MScores'] = array(
|
|
|
|
- '#type' => 'select',
|
|
|
|
- '#title' => t('Match/Mismatch Scores:'),
|
|
|
|
- '#options' => array(
|
|
|
|
- 0 => t('No adjustment'),
|
|
|
|
- 1 => t('Composition-based statistics'),
|
|
|
|
- 2 => t('Conditional compositional score matrix adjustment'),
|
|
|
|
- 3 => t('Universal composition score matrix adjustment '),
|
|
|
|
- ),
|
|
|
|
- '#default_value' => 2,
|
|
|
|
- '#description' => t('Matrix adjustment method to compensate for amino acid composition of sequences'),
|
|
|
|
-);
|
|
|
|
-
|
|
|
|
-//Submit
|
|
|
|
-
|
|
|
|
-$form['submit'] = array(
|
|
|
|
- '#type' => 'submit',
|
|
|
|
- '#default_value' => ' BLAST ',
|
|
|
|
-);
|
|
|
|
-
|
|
|
|
-return $form;
|
|
|
|
|
|
+ ),
|
|
|
|
+ );
|
|
|
|
+
|
|
|
|
+ $form['ALG']['SParam']['gapCost'] = array(
|
|
|
|
+ '#type' => 'select',
|
|
|
|
+ '#title' => t('Gap Costs:'),
|
|
|
|
+ '#prefix' => '<div id="dropdown-second-replace">',
|
|
|
|
+ '#suffix' => '</div>',
|
|
|
|
+ '#options' => _ajax_example_get_second_dropdown_options($selected),
|
|
|
|
+ '#default_value' => 2,
|
|
|
|
+ '#description' => t('Cost to create and extend a gap in an alignment.'),
|
|
|
|
+ );
|
|
|
|
+
|
|
|
|
+ $form['ALG']['SParam']['M&MScores'] = array(
|
|
|
|
+ '#type' => 'select',
|
|
|
|
+ '#title' => t('Match/Mismatch Scores:'),
|
|
|
|
+ '#options' => array(
|
|
|
|
+ 0 => t('No adjustment'),
|
|
|
|
+ 1 => t('Composition-based statistics'),
|
|
|
|
+ 2 => t('Conditional compositional score matrix adjustment'),
|
|
|
|
+ 3 => t('Universal composition score matrix adjustment '),
|
|
|
|
+ ),
|
|
|
|
+ '#default_value' => 2,
|
|
|
|
+ '#description' => t('Matrix adjustment method to compensate for amino acid composition of sequences'),
|
|
|
|
+ );
|
|
|
|
+
|
|
|
|
+ //Submit
|
|
|
|
+ $form['submit'] = array(
|
|
|
|
+ '#type' => 'submit',
|
|
|
|
+ '#default_value' => ' BLAST ',
|
|
|
|
+ );
|
|
|
|
+
|
|
|
|
+ return $form;
|
|
}
|
|
}
|
|
|
|
|
|
/**
|
|
/**
|
|
@@ -198,40 +199,49 @@ return $form;
|
|
* @see blast_protein_form_validate()
|
|
* @see blast_protein_form_validate()
|
|
*/
|
|
*/
|
|
function blast_protein_form_validate($form, &$form_state) {
|
|
function blast_protein_form_validate($form, &$form_state) {
|
|
- $fastaSeq = $form_state['input']['FASTA'];
|
|
|
|
- if (isset($fastaSeq)) {
|
|
|
|
- if(validateFasta($fastaSeq)){
|
|
|
|
- form_set_error('pBLAST', t('Error: Failed to read the Blast query: Wrong format provided for FASTA protein sequence'));
|
|
|
|
- } else {
|
|
|
|
- $form_state['qFlag'] = 'seqQuery';
|
|
|
|
- }
|
|
|
|
- }
|
|
|
|
- $upQuery = file_save_upload('UPLOAD', array('file_validate_extensions' => array('txt fasta fa fna')), FILE_EXISTS_RENAME);
|
|
|
|
- if ($upQuery) {
|
|
|
|
- $upQuery_uri = $upQuery->uri;
|
|
|
|
- $form_state['upQuery_path'] = drupal_realpath($upQuery_uri);
|
|
|
|
- $upQuery_content = file_get_contents($form_state['upQuery_path']);
|
|
|
|
- if(validateFasta($upQuery_content)){
|
|
|
|
- form_set_error('pBLAST', t('Error: Failed to upload the Blast query: Wrong format provided for FASTA protein sequence'));
|
|
|
|
- } else {
|
|
|
|
- $form_state['qFlag'] = 'upQuery';
|
|
|
|
- }
|
|
|
|
- }
|
|
|
|
-
|
|
|
|
- $upDB = file_save_upload('DBUPLOAD', array('file_validate_extensions' => array('txt fasta fa fna')), FILE_EXISTS_RENAME);
|
|
|
|
- if ($upDB) {
|
|
|
|
- $upDB_uri = $upDB->uri;
|
|
|
|
- $form_state['upDB_path'] = drupal_realpath($upDB_uri);
|
|
|
|
- $upDB_content = file_get_contents($form_state['upDB_path']);
|
|
|
|
- if(validateFasta($upDB_content)){
|
|
|
|
- form_set_error('DB', t('Error: Failed to upload the Blast subject sequence file: Wrong format provided for FASTA protein sequence'));
|
|
|
|
- } else {
|
|
|
|
- $form_state['dbFlag'] = 'upQuery';
|
|
|
|
- }
|
|
|
|
- } else {
|
|
|
|
- $form_state['dbFlag'] = 'blastdb';
|
|
|
|
- }
|
|
|
|
-
|
|
|
|
|
|
+
|
|
|
|
+ // Validate query sequence
|
|
|
|
+ $fastaSeq = $form_state['input']['FASTA'];
|
|
|
|
+ if (isset($fastaSeq)) {
|
|
|
|
+ if (validateFasta($fastaSeq)){
|
|
|
|
+ form_set_error('pBLAST', t('Error: Failed to read the Blast query: Wrong format provided for FASTA protein sequence'));
|
|
|
|
+ }
|
|
|
|
+ else {
|
|
|
|
+ $form_state['qFlag'] = 'seqQuery';
|
|
|
|
+ }
|
|
|
|
+ }
|
|
|
|
+
|
|
|
|
+ // Validate Query Upload
|
|
|
|
+ $upQuery = file_save_upload('UPLOAD', array('file_validate_extensions' => array('txt fasta fa fna')), FILE_EXISTS_RENAME);
|
|
|
|
+ if ($upQuery) {
|
|
|
|
+ $upQuery_uri = $upQuery->uri;
|
|
|
|
+ $form_state['upQuery_path'] = drupal_realpath($upQuery_uri);
|
|
|
|
+ $upQuery_content = file_get_contents($form_state['upQuery_path']);
|
|
|
|
+ if(validateFasta($upQuery_content)){
|
|
|
|
+ form_set_error('pBLAST', t('Error: Failed to upload the Blast query: Wrong format provided for FASTA protein sequence'));
|
|
|
|
+ }
|
|
|
|
+ else {
|
|
|
|
+ $form_state['qFlag'] = 'upQuery';
|
|
|
|
+ }
|
|
|
|
+ }
|
|
|
|
+
|
|
|
|
+ // Validate uploaded database
|
|
|
|
+ $upDB = file_save_upload('DBUPLOAD', array('file_validate_extensions' => array('txt fasta fa fna')), FILE_EXISTS_RENAME);
|
|
|
|
+ if ($upDB) {
|
|
|
|
+ $upDB_uri = $upDB->uri;
|
|
|
|
+ $form_state['upDB_path'] = drupal_realpath($upDB_uri);
|
|
|
|
+ $upDB_content = file_get_contents($form_state['upDB_path']);
|
|
|
|
+ if(validateFasta($upDB_content)){
|
|
|
|
+ form_set_error('DB', t('Error: Failed to upload the Blast subject sequence file: Wrong format provided for FASTA protein sequence'));
|
|
|
|
+ }
|
|
|
|
+ else {
|
|
|
|
+ $form_state['dbFlag'] = 'upQuery';
|
|
|
|
+ }
|
|
|
|
+ }
|
|
|
|
+ else {
|
|
|
|
+ $form_state['dbFlag'] = 'blastdb';
|
|
|
|
+ }
|
|
|
|
+
|
|
}
|
|
}
|
|
|
|
|
|
/**
|
|
/**
|
|
@@ -239,19 +249,21 @@ function blast_protein_form_validate($form, &$form_state) {
|
|
*
|
|
*
|
|
* @see blast_protein_form_submit()
|
|
* @see blast_protein_form_submit()
|
|
*/
|
|
*/
|
|
-function blast_protein_form_submit($form, &$form_state) {
|
|
|
|
|
|
+function blast_protein_form_submit($form, &$form_state) {
|
|
|
|
+
|
|
$eVal = $form_state['values']['eVal'];
|
|
$eVal = $form_state['values']['eVal'];
|
|
-
|
|
|
|
|
|
+
|
|
$trgtKey = $form_state['values']['maxTarget'];
|
|
$trgtKey = $form_state['values']['maxTarget'];
|
|
$numAlign = $form['ALG']['GParam']['maxTarget']['#options'][$trgtKey];
|
|
$numAlign = $form['ALG']['GParam']['maxTarget']['#options'][$trgtKey];
|
|
-
|
|
|
|
|
|
+
|
|
$wsKey = $form_state['values']['wordSize'];
|
|
$wsKey = $form_state['values']['wordSize'];
|
|
$wordSize = $form['ALG']['GParam']['wordSize']['#options'][$wsKey];
|
|
$wordSize = $form['ALG']['GParam']['wordSize']['#options'][$wsKey];
|
|
-
|
|
|
|
|
|
+
|
|
|
|
+ // Expand Gap Cost key into open and extend penalties
|
|
$gapKey = $form_state['values']['MATRIX'];
|
|
$gapKey = $form_state['values']['MATRIX'];
|
|
switch ($gapKey) {
|
|
switch ($gapKey) {
|
|
case 0:
|
|
case 0:
|
|
- $matrix ="PAM30";
|
|
|
|
|
|
+ $matrix ="PAM30";
|
|
$gapKey = $form_state['values']['gapCost'];
|
|
$gapKey = $form_state['values']['gapCost'];
|
|
switch ($gapKey) {
|
|
switch ($gapKey) {
|
|
case 0:
|
|
case 0:
|
|
@@ -311,7 +323,7 @@ function blast_protein_form_submit($form, &$form_state) {
|
|
}
|
|
}
|
|
break;
|
|
break;
|
|
case 2:
|
|
case 2:
|
|
- $matrix ="PAM250";
|
|
|
|
|
|
+ $matrix ="PAM250";
|
|
$gapKey = $form_state['values']['gapCost'];
|
|
$gapKey = $form_state['values']['gapCost'];
|
|
switch ($gapKey) {
|
|
switch ($gapKey) {
|
|
case 0:
|
|
case 0:
|
|
@@ -407,7 +419,7 @@ function blast_protein_form_submit($form, &$form_state) {
|
|
}
|
|
}
|
|
break;
|
|
break;
|
|
case 4:
|
|
case 4:
|
|
- $matrix ="BLOSUM62";
|
|
|
|
|
|
+ $matrix ="BLOSUM62";
|
|
$gapKey = $form_state['values']['gapCost'];
|
|
$gapKey = $form_state['values']['gapCost'];
|
|
switch ($gapKey) {
|
|
switch ($gapKey) {
|
|
case 0:
|
|
case 0:
|
|
@@ -457,7 +469,7 @@ function blast_protein_form_submit($form, &$form_state) {
|
|
}
|
|
}
|
|
break;
|
|
break;
|
|
case 5:
|
|
case 5:
|
|
- $matrix ="BLOSUM45";
|
|
|
|
|
|
+ $matrix ="BLOSUM45";
|
|
$gapKey = $form_state['values']['gapCost'];
|
|
$gapKey = $form_state['values']['gapCost'];
|
|
switch ($gapKey) {
|
|
switch ($gapKey) {
|
|
case 0:
|
|
case 0:
|
|
@@ -511,7 +523,7 @@ function blast_protein_form_submit($form, &$form_state) {
|
|
}
|
|
}
|
|
break;
|
|
break;
|
|
case 6:
|
|
case 6:
|
|
- $matrix ="BLOSUM50";
|
|
|
|
|
|
+ $matrix ="BLOSUM50";
|
|
$gapKey = $form_state['values']['gapCost'];
|
|
$gapKey = $form_state['values']['gapCost'];
|
|
switch ($gapKey) {
|
|
switch ($gapKey) {
|
|
case 0:
|
|
case 0:
|
|
@@ -577,7 +589,7 @@ function blast_protein_form_submit($form, &$form_state) {
|
|
}
|
|
}
|
|
break;
|
|
break;
|
|
case 7:
|
|
case 7:
|
|
- $matrix ="BLOSUM90";
|
|
|
|
|
|
+ $matrix ="BLOSUM90";
|
|
$gapKey = $form_state['values']['gapCost'];
|
|
$gapKey = $form_state['values']['gapCost'];
|
|
switch ($gapKey) {
|
|
switch ($gapKey) {
|
|
case 0:
|
|
case 0:
|
|
@@ -611,33 +623,44 @@ function blast_protein_form_submit($form, &$form_state) {
|
|
}
|
|
}
|
|
break;
|
|
break;
|
|
}
|
|
}
|
|
-
|
|
|
|
|
|
+
|
|
|
|
+ // If the query was submitted via the texrfield then create a file containing it
|
|
if ( isset($form_state['qFlag']) ) {
|
|
if ( isset($form_state['qFlag']) ) {
|
|
- if ( $form_state['qFlag'] == 'seqQuery' ) {
|
|
|
|
- $seq_content = $form_state['values']['FASTA'];
|
|
|
|
- $query = "/tmp/user__query_file.fasta";
|
|
|
|
- file_put_contents ( $query , $seq_content);
|
|
|
|
- } elseif ( $form_state['qFlag'] == 'upQuery' ) {
|
|
|
|
- $query = $form_state['upQuery_path'];
|
|
|
|
- }
|
|
|
|
|
|
+ if ( $form_state['qFlag'] == 'seqQuery' ) {
|
|
|
|
+ $seq_content = $form_state['values']['FASTA'];
|
|
|
|
+ $query = "/tmp/user__query_file.fasta";
|
|
|
|
+ file_put_contents ( $query , $seq_content);
|
|
|
|
+ }
|
|
|
|
+ elseif ( $form_state['qFlag'] == 'upQuery' ) {
|
|
|
|
+ $query = $form_state['upQuery_path'];
|
|
|
|
+ }
|
|
}
|
|
}
|
|
-
|
|
|
|
|
|
+
|
|
|
|
+ // If the BLAST database was uploaded then use it to run the BLAST
|
|
if ($form_state['dbFlag'] == 'upQuery') {
|
|
if ($form_state['dbFlag'] == 'upQuery') {
|
|
|
|
+
|
|
$subjectSeq = $form_state['upDB_path'];
|
|
$subjectSeq = $form_state['upDB_path'];
|
|
$subSeqOut = drupal_basename($form_state['upDB_path']) . rand(0, 10000);
|
|
$subSeqOut = drupal_basename($form_state['upDB_path']) . rand(0, 10000);
|
|
|
|
+
|
|
$blast_subj_cmd = "blastp -task blastp -query $query -subject $subjectSeq -out sites/default/files/$subSeqOut.blastp.html";
|
|
$blast_subj_cmd = "blastp -task blastp -query $query -subject $subjectSeq -out sites/default/files/$subSeqOut.blastp.html";
|
|
-// -evalue $eVal -word_size $wordSize -gapopen $gapOpen -gapextend $gapExtend -matrix $matrix -num_alignments 100 -html";
|
|
|
|
system($blast_subj_cmd);
|
|
system($blast_subj_cmd);
|
|
- } elseif ($form_state['dbFlag'] == 'blastdb') {
|
|
|
|
|
|
+
|
|
|
|
+ }
|
|
|
|
+ // Otherwise, we are using one of the website provided BLAST databases so form the
|
|
|
|
+ // BLAST command accordingly
|
|
|
|
+ elseif ($form_state['dbFlag'] == 'blastdb') {
|
|
|
|
+
|
|
$selected_db = $form_state['values']['SELECT_DB'];
|
|
$selected_db = $form_state['values']['SELECT_DB'];
|
|
$blastdb_node = node_load($selected_db);
|
|
$blastdb_node = node_load($selected_db);
|
|
- $blastdb_path = $blastdb_node->db_path;
|
|
|
|
|
|
+ $blastdb_path = $blastdb_node->db_path;
|
|
$blastdb_human_name = $form['DB']['SELECT_DB']['#options'][$selected_db];
|
|
$blastdb_human_name = $form['DB']['SELECT_DB']['#options'][$selected_db];
|
|
$subSeqOut = str_replace(' ','_', $blastdb_human_name) . rand(0, 10000);
|
|
$subSeqOut = str_replace(' ','_', $blastdb_human_name) . rand(0, 10000);
|
|
|
|
+
|
|
$blast_db_cmd = "blastp -task blastp -query $query -db $blastdb_path -out sites/default/files/$subSeqOut.blastp.html -evalue $eVal -word_size $wordSize -gapopen $gapOpen -gapextend $gapExtend -matrix $matrix -num_alignments 100 -html";
|
|
$blast_db_cmd = "blastp -task blastp -query $query -db $blastdb_path -out sites/default/files/$subSeqOut.blastp.html -evalue $eVal -word_size $wordSize -gapopen $gapOpen -gapextend $gapExtend -matrix $matrix -num_alignments 100 -html";
|
|
system($blast_db_cmd,$input);
|
|
system($blast_db_cmd,$input);
|
|
}
|
|
}
|
|
|
|
|
|
|
|
+ // Redirect to the BLAST Results page
|
|
$path = "$subSeqOut.blastp.html";
|
|
$path = "$subSeqOut.blastp.html";
|
|
drupal_goto("blast/report/$path");
|
|
drupal_goto("blast/report/$path");
|
|
}
|
|
}
|
|
@@ -646,46 +669,46 @@ function blast_protein_form_submit($form, &$form_state) {
|
|
/**
|
|
/**
|
|
* FASTA validating parser
|
|
* FASTA validating parser
|
|
*
|
|
*
|
|
- * @param $sequence
|
|
|
|
- * A string of characters to be validated. A sequence in FASTA format begins with a single-line description, followed by lines of sequence data.
|
|
|
|
- * The description line is distinguished from the sequence data by a greater-than (">") symbol in the first column.
|
|
|
|
- * The word following the ">" symbol is the identifier of the sequence, and the rest of the line is the description (both are optional).
|
|
|
|
- * There should be no space between the ">" and the first letter of the identifier. The sequence ends if another line starting with a ">" appears;
|
|
|
|
|
|
+ * @param $sequence
|
|
|
|
+ * A string of characters to be validated. A sequence in FASTA format begins with a single-line description, followed by lines of sequence data.
|
|
|
|
+ * The description line is distinguished from the sequence data by a greater-than (">") symbol in the first column.
|
|
|
|
+ * The word following the ">" symbol is the identifier of the sequence, and the rest of the line is the description (both are optional).
|
|
|
|
+ * There should be no space between the ">" and the first letter of the identifier. The sequence ends if another line starting with a ">" appears;
|
|
* this indicates the start of another sequence.
|
|
* this indicates the start of another sequence.
|
|
*
|
|
*
|
|
* @return
|
|
* @return
|
|
- * Return a boolean. 1 if the sequence does not pass the format valifation stage and 0 otherwise.
|
|
|
|
|
|
+ * Return a boolean. 1 if the sequence does not pass the format valifation stage and 0 otherwise.
|
|
*
|
|
*
|
|
*/
|
|
*/
|
|
function validateFasta($sequence) {
|
|
function validateFasta($sequence) {
|
|
$fastaIdRegEx = '/^>.*(\\n|\\r)/';
|
|
$fastaIdRegEx = '/^>.*(\\n|\\r)/';
|
|
$fastaSeqRegEx = '/[^acgturykmswbdhvnxACGTURYKMSWBDHVNX\*\-\n\r]/';
|
|
$fastaSeqRegEx = '/[^acgturykmswbdhvnxACGTURYKMSWBDHVNX\*\-\n\r]/';
|
|
if ( preg_match($fastaSeqRegEx,$sequence) && !(preg_match($fastaIdRegEx,$sequence)) ) {
|
|
if ( preg_match($fastaSeqRegEx,$sequence) && !(preg_match($fastaIdRegEx,$sequence)) ) {
|
|
- $flag = 1;
|
|
|
|
- } else {
|
|
|
|
- $flag = 0;
|
|
|
|
- }
|
|
|
|
-
|
|
|
|
|
|
+ $flag = 1;
|
|
|
|
+ } else {
|
|
|
|
+ $flag = 0;
|
|
|
|
+ }
|
|
|
|
+
|
|
return $flag;
|
|
return $flag;
|
|
}
|
|
}
|
|
|
|
|
|
|
|
|
|
/**
|
|
/**
|
|
- * Load the preexisting blast databases
|
|
|
|
|
|
+ * Generate an array of BLAST database options based on existing nodes of type BlastDB
|
|
*
|
|
*
|
|
* @return
|
|
* @return
|
|
- * Return human readble names of the preexisting blast databases
|
|
|
|
|
|
+ * Return human readble names of the pre-existing blast databases
|
|
*
|
|
*
|
|
*/
|
|
*/
|
|
function DB_options() {
|
|
function DB_options() {
|
|
$type = 'blastdb';
|
|
$type = 'blastdb';
|
|
$nodes = node_load_multiple(array(), array('type'=> $type));
|
|
$nodes = node_load_multiple(array(), array('type'=> $type));
|
|
-
|
|
|
|
|
|
+
|
|
$options = array();
|
|
$options = array();
|
|
|
|
|
|
foreach ($nodes as $node) {
|
|
foreach ($nodes as $node) {
|
|
if ( isset($node) && isset($node->db_dbtype) ) {
|
|
if ( isset($node) && isset($node->db_dbtype) ) {
|
|
- if ( ($node->db_dbtype=='p') ) {
|
|
|
|
|
|
+ if ( ($node->db_dbtype=='p') ) {
|
|
$options[$node->nid] = $node->db_name;
|
|
$options[$node->nid] = $node->db_name;
|
|
}
|
|
}
|
|
}
|
|
}
|
|
@@ -698,15 +721,14 @@ function DB_options() {
|
|
}
|
|
}
|
|
|
|
|
|
/**
|
|
/**
|
|
- * Fill the first dropdown list with appropriate options
|
|
|
|
|
|
+ * Fill the first dropdown list with appropriate options
|
|
*
|
|
*
|
|
* @return
|
|
* @return
|
|
- * An array consisting of matrices name for the first dropdown list
|
|
|
|
- *
|
|
|
|
|
|
+ * An array consisting of matrices name for the first dropdown list
|
|
*/
|
|
*/
|
|
function _ajax_example_get_first_dropdown_options() {
|
|
function _ajax_example_get_first_dropdown_options() {
|
|
return drupal_map_assoc(array(
|
|
return drupal_map_assoc(array(
|
|
- t('PAM30'),
|
|
|
|
|
|
+ t('PAM30'),
|
|
t('PAM70'),
|
|
t('PAM70'),
|
|
t('PAM250'),
|
|
t('PAM250'),
|
|
t('BLOSUM80'),
|
|
t('BLOSUM80'),
|
|
@@ -718,11 +740,10 @@ function _ajax_example_get_first_dropdown_options() {
|
|
}
|
|
}
|
|
|
|
|
|
/**
|
|
/**
|
|
- * Fill the second dropdown list with appropriate options
|
|
|
|
|
|
+ * Fill the second dropdown list with appropriate options
|
|
*
|
|
*
|
|
* @return
|
|
* @return
|
|
- * An array containing open and extension gap values for the chosen matrix (to fill the second dropdown list)
|
|
|
|
- *
|
|
|
|
|
|
+ * An array containing open and extension gap values for the chosen matrix (to fill the second dropdown list)
|
|
*/
|
|
*/
|
|
function _ajax_example_get_second_dropdown_options($key = '') {
|
|
function _ajax_example_get_second_dropdown_options($key = '') {
|
|
$options = array(
|
|
$options = array(
|
|
@@ -829,11 +850,7 @@ function _ajax_example_get_second_dropdown_options($key = '') {
|
|
}
|
|
}
|
|
|
|
|
|
/**
|
|
/**
|
|
- * Respond to Ajax dropdown call
|
|
|
|
- *
|
|
|
|
- * @return
|
|
|
|
- *
|
|
|
|
- *
|
|
|
|
|
|
+ * Respond to Ajax dropdown call
|
|
*/
|
|
*/
|
|
function ajax_example_dependent_dropdown_callback($form, $form_state) {
|
|
function ajax_example_dependent_dropdown_callback($form, $form_state) {
|
|
return $form['ALG']['SParam']['gapCost'];
|
|
return $form['ALG']['SParam']['gapCost'];
|