1 new changeset in galaxy-central:
http://bitbucket.org/galaxy/galaxy-central/changeset/10c1c31c03bc/
changeset: 10c1c31c03bc
branches:
user: greg
date: 2011-05-31 18:03:52
summary: Fix for importing library datasets into a new named history. Library datasets will no longer be imported into both the current history and the new named history ( only the latter ) in this scenario.
affected #: 1 file (53 bytes)
Repository URL: https://bitbucket.org/galaxy/galaxy-central/
--
This is a commit notification from bitbucket.org. You are receiving
this because you have the service enabled, addressing the recipient of
this email.
1 new changeset in galaxy-central:
http://bitbucket.org/galaxy/galaxy-central/changeset/3e7f7c40b4fe/
changeset: 3e7f7c40b4fe
branches:
user: richard_burhans
date: 2011-05-27 15:22:27
summary: Updated help text for genome diversity tools
affected #: 3 files (501 bytes)
--- a/tools/genome_diversity/extract_flanking_dna.xml Fri May 27 11:19:16 2011 +0200
+++ b/tools/genome_diversity/extract_flanking_dna.xml Fri May 27 09:22:27 2011 -0400
@@ -54,7 +54,9 @@
<help>
**What it does**
- If extracts flanking dna for the SNPs in the dataset.
+ It reports a DNA segment containing each SNP, with up to 200 nucleotides on
+ either side of the SNP position, which is indicated by "n". Fewer nucleotides
+ are reported if the SNP is near an end of the assembled genome fragment.
-----
--- a/tools/genome_diversity/select_restriction_enzymes.xml Fri May 27 11:19:16 2011 +0200
+++ b/tools/genome_diversity/select_restriction_enzymes.xml Fri May 27 09:22:27 2011 -0400
@@ -57,7 +57,9 @@
<help>
**What it does**
- It selects SNPs in the dataset by restriction enzyme.
+ It selects the SNPs that are differentially cut by at least one of the
+ specified restriction enzymes. The enzymes are required to cut the amplified
+ segment (for the specified PCR primers) only at the SNP.
-----
--- a/tools/genome_diversity/select_snps.xml Fri May 27 11:19:16 2011 +0200
+++ b/tools/genome_diversity/select_snps.xml Fri May 27 09:22:27 2011 -0400
@@ -51,7 +51,9 @@
<help>
**What it does**
- If selects a specified number of SNPs from the dataset.
+ It attempts to select a specified number of SNPs from the dataset, making them
+ approximately uniformly spaced relative to the reference genome. The number
+ actually selected may be slightly more than the specified number.
-----
Repository URL: https://bitbucket.org/galaxy/galaxy-central/
--
This is a commit notification from bitbucket.org. You are receiving
this because you have the service enabled, addressing the recipient of
this email.
1 new changeset in galaxy-central:
http://bitbucket.org/galaxy/galaxy-central/changeset/cc686fadddbc/
changeset: cc686fadddbc
branches:
user: kanwei
date: 2011-05-27 10:47:17
summary: [Peter Cock] Filtering tool: Only cast used columns. Fix valid lines and skipped info lines. [Kanwei] Add new testcases to test this behavior. Fixes #537. Fixes #535
affected #: 4 files (364 bytes)
--- a/tools/stats/filtering.py Thu May 26 13:29:00 2011 -0400
+++ b/tools/stats/filtering.py Fri May 27 10:47:17 2011 +0200
@@ -32,7 +32,7 @@
cond_text = sys.argv[3]
try:
in_columns = int( sys.argv[4] )
- assert sys.argv[5] #check to see that the column types varaible isn't null
+ assert sys.argv[5] #check to see that the column types variable isn't null
in_column_types = sys.argv[5].split( ',' )
except:
stop_err( "Data does not appear to be tabular. This tool can only be used with tab-delimited data." )
@@ -60,22 +60,25 @@
except:
if operand in secured:
stop_err( "Illegal value '%s' in condition '%s'" % ( operand, cond_text ) )
-
-# Find the largest column used in the filter.
-largest_col_index = -1
-for match in re.finditer( 'c(\d)+', cond_text ):
- col_index = int( match.group()[1:] )
- if col_index > largest_col_index:
- largest_col_index = col_index
+
+# Work out which columns are used in the filter (save using 1 based counting)
+used_cols = sorted(set(int(match.group()[1:]) \
+ for match in re.finditer('c(\d)+', cond_text)))
+largest_col_index = max(used_cols)
# Prepare the column variable names and wrappers for column data types. Only
-# prepare columns up to largest column in condition.
+# cast columns used in the filter.
cols, type_casts = [], []
for col in range( 1, largest_col_index + 1 ):
col_name = "c%d" % col
cols.append( col_name )
col_type = in_column_types[ col - 1 ]
- type_cast = "%s(%s)" % ( col_type, col_name )
+ if col in used_cols:
+ type_cast = "%s(%s)" % ( col_type, col_name )
+ else:
+ #If we don't use this column, don't cast it.
+ #Otherwise we get errors on things like optional integer columns.
+ type_cast = col_name
type_casts.append( type_cast )
col_str = ', '.join( cols ) # 'c1, c2, c3, c4'
@@ -83,6 +86,7 @@
assign = "%s, = line.split( '\\t' )[:%i]" % ( col_str, largest_col_index )
wrap = "%s = %s" % ( col_str, type_cast_str )
skipped_lines = 0
+invalid_lines = 0
first_invalid_line = 0
invalid_line = None
lines_kept = 0
@@ -96,9 +100,6 @@
line = line.rstrip( '\\r\\n' )
if not line or line.startswith( '#' ):
skipped_lines += 1
- if not invalid_line:
- first_invalid_line = i + 1
- invalid_line = line
continue
try:
%s
@@ -107,7 +108,7 @@
lines_kept += 1
print >> out, line
except:
- skipped_lines += 1
+ invalid_lines += 1
if not invalid_line:
first_invalid_line = i + 1
invalid_line = line
@@ -132,5 +133,7 @@
print 'kept %4.2f%% of %d lines.' % ( 100.0*lines_kept/valid_lines, total_lines )
else:
print 'Possible invalid filter condition "%s" or non-existent column referenced. See tool tips, syntax and examples.' % cond_text
- if skipped_lines > 0:
- print 'Skipped %d invalid lines starting at line #%d: "%s"' % ( skipped_lines, first_invalid_line, invalid_line )
+ if invalid_lines:
+ print 'Skipped %d invalid line(s) starting at line #%d: "%s"' % ( invalid_lines, first_invalid_line, invalid_line )
+ if skipped_lines:
+ print 'Skipped %i comment (starting with #) or blank line(s)' % skipped_lines
--- a/tools/stats/filtering.xml Thu May 26 13:29:00 2011 -0400
+++ b/tools/stats/filtering.xml Fri May 27 10:47:17 2011 +0200
@@ -1,4 +1,4 @@
-<tool id="Filter1" name="Filter" version="1.0.1">
+<tool id="Filter1" name="Filter" version="1.1.0"><description>data on any column using simple expressions</description><command interpreter="python">
filtering.py $input $out_file1 "$cond" ${input.metadata.columns} "${input.metadata.column_types}"
@@ -29,7 +29,11 @@
<param name="cond" value="c3=='chr1' and c5>5"/><output name="out_file1" file="filter1_test3.sam"/></test>
-
+ <test>
+ <param name="input" value="filter1_inbad.bed"/>
+ <param name="cond" value="c1=='chr22'"/>
+ <output name="out_file1" file="filter1_test4.bed"/>
+ </test></tests><help>
Repository URL: https://bitbucket.org/galaxy/galaxy-central/
--
This is a commit notification from bitbucket.org. You are receiving
this because you have the service enabled, addressing the recipient of
this email.
1 new changeset in galaxy-central:
http://bitbucket.org/galaxy/galaxy-central/changeset/bae2efef9ac7/
changeset: r5602:bae2efef9ac7
user: richard_burhans
date: 2011-05-26 19:29:00
summary: updated help section for the genome diversity extract primers tool
affected #: 1 file (153 bytes)
--- a/tools/genome_diversity/extract_primers.xml Thu May 26 08:27:52 2011 -0400
+++ b/tools/genome_diversity/extract_primers.xml Thu May 26 13:29:00 2011 -0400
@@ -49,7 +49,16 @@
<help>
**What it does**
- It extracts primers for the SNPs in the dataset.
+ This tool extracts primers for SNPs in the argument dataset using the
+ Primer3 program. The first line of output for a given SNP reports the
+ name of the assembled contig, the SNP's position in the contig, the
+ two variant nucleotides, and Primer3's "pair penalty". The next line,
+ if not blank, names restriction enzymes (from the user-adjustable list)
+ that differentially cut at that site, but do not cut at any other position
+ between and including the primer positions. The next lines show the SNP's
+ flanking regions, with the SNP position indicated by "n", including the
+ primer positions and an additional 3 nucleotides.
+
-----
@@ -57,30 +66,26 @@
- input file::
- chr2_75111355_75112576 314 A C L F chr2 75111676 C F 15 4 53 2 9 48 Y 96 0.369 0.355 0.396 0
- chr8_93901796_93905612 2471 A C A A chr8 93904264 A A 8 0 51 10 2 14 Y 961 0.016 0.534 0.114 2
- chr10_7434473_7435447 524 T C S S chr10 7435005 T S 11 5 90 14 0 69 Y 626 0.066 0.406 0.727 0
- chr14_80021455_80022064 138 G A H H chr14 80021593 G H 14 0 69 9 6 124 Y 377 0.118 0.997 0.195 1
- chr15_64470252_64471048 89 G A Y Y chr15 64470341 G Y 5 6 109 14 0 69 Y 312 0.247 0.998 0.393 0
- chr18_48070585_48071386 514 C T E K chr18 48071100 T K 7 7 46 14 0 69 Y 2 0.200 0.032 0.163 0
- chr18_50154905_50155664 304 A G Y C chr18 50155208 A Y 4 2 17 5 1 22 Y 8 0.022 0.996 0.128 0
- chr18_57379354_57380496 315 C T V V chr18 57379669 G V 11 0 60 9 6 62 Y 726 0.118 0.048 0.014 1
- chr19_14240610_14242055 232 C T A V chr19 14240840 C A 18 8 56 15 5 42 Y 73 0.003 0.153 0.835 0
- chr19_39866997_39874915 3117 C T P P chr19 39870110 C P 3 7 65 14 2 32 Y 6 0.321 0.911 0.462 4
- etc.
+ chr5_30800874_30802049 734 G A chr5 30801606 A 24 0 99 4 11 97 Y 496 0.502 0.033 0.215 6
+ chr8_55117827_55119487 994 A G chr8 55118815 G 25 0 102 4 11 96 Y 22 0.502 0.025 2.365 1
+ chr9_100484836_100485311 355 C T chr9 100485200 T 27 0 108 6 17 100 Y 190 0.512 0.880 2.733 4
+ chr12_3635530_3637738 2101 T C chr12 3637630 T 25 0 102 4 13 93 Y 169 0.554 0.024 0.366 4
- output file::
- > chr2_75111355_75112576 314 A C 0.613831
+ chr5_30800874_30802049 734 G A 0.352964
+ BglII,MboI,Sau3AI,Tru9I,XhoII
+ 1 CTGAAGGTGAGCAGGATTCAGGAGACAGAAAACAAAGCCCAGGCCTGCCCAAGGTGGAAA
+ >>>>>>>>>>>>>>>>>>>>
- 1 TTTGCCCTGAGGCAGACTTTTTAAAGTACTGTGTAATGTATGAAGTCCTTCTGCTCAAGC
- >>>>>>>>>>>>>>>>>>>>
-
- 61 AAATCATTGGCATGAAAACAGTTGCAAACTTATTGTGAGAGAAGAGTCCAAGAGTTTTAA
-
-
- 121 CAGTCTGTAAGTATATAGCCTGTGAGTTTGATTTCCTTCTTGTTTTTnTTCCAGAAACAT
-
+ 61 AGTCTAACAACTCGCCCTCTGCTTAnATCTGAGACTCACAGGGATAATAACACACTTGGT
+
+
+ 21 CAAGGAATAAACTAGATATTATTCACTCCTCTAGAAGGCTGCCAGGAAAATTGCCTGACT
+ <<<<<<<
+
+ 181 TGAACCTTGGCTCTGA
+ <<<<<<<<<<<<<
etc.
</help></tool>
Repository URL: https://bitbucket.org/galaxy/galaxy-central/
--
This is a commit notification from bitbucket.org. You are receiving
this because you have the service enabled, addressing the recipient of
this email.