/var/www/DASR_2011/filerepo/4db5770b732ff/4db5770d85ea1_overlapbedfile.txt [overlap_combined_columns] => /var/www/DASR_2011/filerepo/4db5770b732ff/4db5770d38c46_combinedcolumns.txt [overlaps_chart1] => https://chart.googleapis.com/chart?cht=bvs&chxt=x,y&chm=N,000000,0,-1,11&chxl=0:|chr1|chr2|chr3|chr4|chr5|chr6|chr7|chr8|chr9|chr10|chr11|chr12|chr13|chr14|chr15|chr16|chr17|chr18|chr19|chr20|chr21|chr22|chrX|chrY&chd=t:935,996,667,945,972,811,898,820,755,777,978,640,649,563,752,545,588,534,455,0,0,0,171,0&chxr=1,0,996&chds=0,996&chbh=25,15,4&chs=1000x200 [overlaps_chart2] => https://chart.googleapis.com/chart?cht=lc&chxt=x,y&chco=4D89F9,C6D9FD,&chd=t:0.2108,0.1879,0.2023,0.1963,0.1973,0.1898,0.1968,0.1972,0.1867,0.1982,0.1863,0.2097,0.215,0.211,0.2315,0.2095,0.1874,0.2215,0.17,0,0,0,0.1612,0|0.8408,0.8571,0.8056,0.884,0.8408,0.8844,0.9007,0.82,0.8779,0.8993,0.8964,0.8625,0.8484,0.8595,0.8963,0.8385,0.8963,0.8558,0.9137,0,0,0,0.7155,0&chm=N,FF0000,0,-1,11,1.0|N,0000FF,1,-1,11,-1.0&chdl=4db5770b739c9|4db5770b74199&chco=FF0000,0000FF&chxl=0:|chr1|chr2|chr3|chr4|chr5|chr6|chr7|chr8|chr9|chr10|chr11|chr12|chr13|chr14|chr15|chr16|chr17|chr18|chr19|chr20|chr21|chr22|chrX|chrY&chxr=1,0,0.2315&chds=0,0.2315&chbh=25,15,4&chs=1000x200 [temp] => /var/www/DASR_2011/filerepo/asdfasdfa ) */ //print_r($result); $metadata = array(); $db = getDB(); $bed1 = $db->getRow("select name, header, genome_id from bedfile where bedfile_id = ?", array($_GET['bed_id1'])); $bed2 = $db->getRow("select name, header, genome_id from bedfile where bedfile_id = ?", array($_GET['bed_id2'])); $bed1_name = $bed1['name']; //$db->getOne("select name from bedfile where bedfile_id = ?", $_GET['bed_id1']); $bed2_name = $bed2['name']; //$db->getOne("select name from bedfile where bedfile_id = ?", $_GET['bed_id2']); if($bed1['genome_id'] == $bed2['genome_id']) { $new_genome_id = $bed1['genome_id']; $new_genome = $db->getOne("select name from genome where genome_id = ?", $new_genome_id); } else { $new_genome_id = null; $new_genome = ''; } $dbl_col_filepath = FILE_REPOSITORY_PATH . '/' . uniqid() . date('YmdHis'); rename($result['overlap_combined_columns'], $dbl_col_filepath); $header = "track name=\"compare $bed1[name] overlap with $bed2[name]\""; //$bed1['header'] . " " . $bed2['header']; //trim(`head -1 $dbl_col_filepath`); $date = date('Y-m-d H:i:s'); $d = array('name' => "Comparison $bed1_name, $bed2_name, side-by-side comparison", 'descr' => "$header Overlap: $_GET[bp_overlap]"); $db->query("insert into bedfile (experiment_id, name, genome_id, filetype, descr, header, upload_date, overlap_counts_chart_url, overlaps_per_read_chart_url) values (?,?,?,?,?,?,?,?,?)", array($_GET['experiment_id'], $d['name'], $new_genome_id, "rawbed", $d['descr'], $header, $date, $result['overlaps_chart1'], $result['overlaps_chart2'])); $dbl_col_id = $db->getOne("select last_insert_id()"); $db->query("insert into rawbed (bedfile_id, filepath) values (?,?)", array($dbl_col_id, $dbl_col_filepath)); $metadata[] = array('name' => $d['name'], 'header' => $header, 'upload_date' => date('n/j/y', strtotime($date)), 'id' => $dbl_col_id, 'genome' => $new_genome, 'method' => '', 'descr' => $d['descr']); $overlap_filepath = FILE_REPOSITORY_PATH . '/' . uniqid() . date('YmdHis'); rename($result['overlap_bedfile'], $overlap_filepath); //$header = $bed1['header'] . " " . $bed2['header']; //trim(`head -1 $overlap_filepath`); $d = array('name' => "Comparison $bed1_name, $bed2_name, region of overlap", 'descr' => "$header Overlap: $_GET[bp_overlap]"); $db->query("insert into bedfile (experiment_id, name, genome_id, filetype, descr, header, upload_date, overlap_counts_chart_url, overlaps_per_read_chart_url) values (?,?,?,?,?,?,?,?,?)", array($_GET['experiment_id'], $d['name'], $new_genome_id, "rawbed", $d['descr'], $header, $date, $result['overlaps_chart1'], $result['overlaps_chart2'])); $overlap_id = $db->getOne("select last_insert_id()"); $db->query("insert into rawbed (bedfile_id, filepath) values (?,?)", array($overlap_id, $overlap_filepath)); $metadata[] = array('name' => $d['name'], 'header' => $header, 'upload_date' => date('n/j/y', strtotime($date)), 'id' => $overlap_id, 'genome' => $new_genome, 'method' => '', 'descr' => $d['descr']); $db->commit(); rrmdir($result['temp']); print json_encode($metadata); ?>