# HG changeset patch # User bernhardlutz # Date 1403116391 14400 # Node ID 893617de95e5dde90618f9cc0f5bc29e3a81ac81 # Parent d0df760f1d32dc2b5134a6a4450ebb30b3174785 Uploaded diff -r d0df760f1d32 -r 893617de95e5 annotateBed.xml --- a/annotateBed.xml Wed Jun 18 13:13:51 2014 -0400 +++ b/annotateBed.xml Wed Jun 18 14:33:11 2014 -0400 @@ -13,7 +13,7 @@ $bed.input #end for - #if names.names_select == 'yes': + #if $names.names_select == 'yes': -names #for $bed in $names.beds: $bed.inputName diff -r d0df760f1d32 -r 893617de95e5 closestBed.xml --- a/closestBed.xml Wed Jun 18 13:13:51 2014 -0400 +++ b/closestBed.xml Wed Jun 18 14:33:11 2014 -0400 @@ -7,7 +7,6 @@ closestBed - $split $strand $addition -t $ties @@ -20,13 +19,13 @@ - - - + + + - - + + diff -r d0df760f1d32 -r 893617de95e5 getfastaBed.xml --- a/getfastaBed.xml Wed Jun 18 13:13:51 2014 -0400 +++ b/getfastaBed.xml Wed Jun 18 14:33:11 2014 -0400 @@ -7,7 +7,6 @@ bedtools getfasta - $name $tab $strand @@ -17,22 +16,20 @@ -fo $output - - - + + - - - + + + - **What it does** -bedtools getfasta will extract the sequence defined by the coordinates in a BED interval and create a new FASTA entry in the output file for each extracted sequence. By default, the FASTA header for each extracted sequence will be formatted as follows: “:-”. +bedtools getfasta will extract the sequence defined by the coordinates in a BED interval and create a new FASTA entry in the output file for each extracted sequence. By default, the FASTA header for each extracted sequence will be formatted as follows: “<chrom>:<start>-<end>”. .. image:: $PATH_TO_IMAGES/getfasta-glyph.png diff -r d0df760f1d32 -r 893617de95e5 shuffleBed.xml --- a/shuffleBed.xml Wed Jun 18 13:13:51 2014 -0400 +++ b/shuffleBed.xml Wed Jun 18 14:33:11 2014 -0400 @@ -30,8 +30,7 @@ - - + @@ -40,6 +39,7 @@ + diff -r d0df760f1d32 -r 893617de95e5 slopBed.xml --- a/slopBed.xml Wed Jun 18 13:13:51 2014 -0400 +++ b/slopBed.xml Wed Jun 18 14:33:11 2014 -0400 @@ -38,7 +38,7 @@ **What it does** -bedtools slop will increase the size of each feature in a feature file by a user-defined number of bases. While something like this could be done with an awk '{OFS="\t" print $1,$2-,$3+}', bedtools slop will restrict the resizing to the size of the chromosome (i.e. no start < 0 and no end > chromosome size). +bedtools slop will increase the size of each feature in a feature file by a user-defined number of bases. While something like this could be done with an awk '{OFS="\t" print $1,$2-<slop>,$3+<slop>}', bedtools slop will restrict the resizing to the size of the chromosome (i.e. no start < 0 and no end > chromosome size). .. image:: $PATH_TO_IMAGES/slop-glyph.png .. class:: warningmark diff -r d0df760f1d32 -r 893617de95e5 tagBed.xml --- a/tagBed.xml Wed Jun 18 13:13:51 2014 -0400 +++ b/tagBed.xml Wed Jun 18 14:33:11 2014 -0400 @@ -35,7 +35,6 @@ - diff -r d0df760f1d32 -r 893617de95e5 test-data/0.bed --- a/test-data/0.bed Wed Jun 18 13:13:51 2014 -0400 +++ b/test-data/0.bed Wed Jun 18 14:33:11 2014 -0400 @@ -1,4 +1,4 @@ -chr1 100 200 -chr1 180 250 -chr1 250 500 -chr1 501 1000 +chr1 100 200 +chr1 180 250 +chr1 250 500 +chr1 501 1000 diff -r d0df760f1d32 -r 893617de95e5 test-data/0_result.bed --- a/test-data/0_result.bed Wed Jun 18 13:13:51 2014 -0400 +++ b/test-data/0_result.bed Wed Jun 18 14:33:11 2014 -0400 @@ -1,2 +1,2 @@ -chr1 100 500 -chr1 501 1000 +chr1 100 500 +chr1 501 1000 diff -r d0df760f1d32 -r 893617de95e5 test-data/1.bed --- a/test-data/1.bed Wed Jun 18 13:13:51 2014 -0400 +++ b/test-data/1.bed Wed Jun 18 14:33:11 2014 -0400 @@ -1,4 +1,4 @@ -chr1 100 200 a1 1 + -chr1 180 250 a2 2 + -chr1 250 500 a3 3 - -chr1 501 1000 a4 4 + +chr1 100 200 a1 1 + +chr1 180 250 a2 2 + +chr1 250 500 a3 3 - +chr1 501 1000 a4 4 + diff -r d0df760f1d32 -r 893617de95e5 test-data/1_result.bed --- a/test-data/1_result.bed Wed Jun 18 13:13:51 2014 -0400 +++ b/test-data/1_result.bed Wed Jun 18 14:33:11 2014 -0400 @@ -1,3 +1,3 @@ -chr1 100 250 + -chr1 501 1000 + -chr1 250 500 - +chr1 100 250 + +chr1 501 1000 + +chr1 250 500 - diff -r d0df760f1d32 -r 893617de95e5 test-data/2.bed --- a/test-data/2.bed Wed Jun 18 13:13:51 2014 -0400 +++ b/test-data/2.bed Wed Jun 18 14:33:11 2014 -0400 @@ -1,4 +1,4 @@ -chr1 100 200 -chr1 180 250 -chr1 250 500 -chr1 501 1000 +chr1 100 200 +chr1 180 250 +chr1 250 500 +chr1 501 1000 diff -r d0df760f1d32 -r 893617de95e5 test-data/2_result.bed --- a/test-data/2_result.bed Wed Jun 18 13:13:51 2014 -0400 +++ b/test-data/2_result.bed Wed Jun 18 14:33:11 2014 -0400 @@ -1,2 +1,2 @@ -chr1 100 500 3 -chr1 501 1000 1 +chr1 100 500 3 +chr1 501 1000 1 diff -r d0df760f1d32 -r 893617de95e5 test-data/3.bed --- a/test-data/3.bed Wed Jun 18 13:13:51 2014 -0400 +++ b/test-data/3.bed Wed Jun 18 14:33:11 2014 -0400 @@ -1,2 +1,2 @@ -chr1 100 200 -chr1 501 1000 +chr1 100 200 +chr1 501 1000 diff -r d0df760f1d32 -r 893617de95e5 test-data/3_result_1000.bed --- a/test-data/3_result_1000.bed Wed Jun 18 13:13:51 2014 -0400 +++ b/test-data/3_result_1000.bed Wed Jun 18 14:33:11 2014 -0400 @@ -1,1 +1,1 @@ -chr1 100 200 1000 +chr1 100 200 1000 diff -r d0df760f1d32 -r 893617de95e5 test-data/annotateBed1.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/annotateBed1.bed Wed Jun 18 14:33:11 2014 -0400 @@ -0,0 +1,4 @@ +chr1 100 200 nasty 1 - +chr2 500 1000 ugly 2 + +chr3 1000 5000 big 3 - + diff -r d0df760f1d32 -r 893617de95e5 test-data/annotateBed2.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/annotateBed2.bed Wed Jun 18 14:33:11 2014 -0400 @@ -0,0 +1,3 @@ +chr1 150 200 geneA 1 + +chr1 175 250 geneB 2 + +chr3 0 10000 geneC 3 - diff -r d0df760f1d32 -r 893617de95e5 test-data/annotateBed3.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/annotateBed3.bed Wed Jun 18 14:33:11 2014 -0400 @@ -0,0 +1,3 @@ +chr1 0 10000 cons1 1 + +chr2 700 10000 cons2 2 - +chr3 4000 10000 cons3 3 + diff -r d0df760f1d32 -r 893617de95e5 test-data/annotateBed4.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/annotateBed4.bed Wed Jun 18 14:33:11 2014 -0400 @@ -0,0 +1,4 @@ +chr1 0 120 known1 - +chr1 150 160 known2 - +chr2 0 10000 known3 + + diff -r d0df760f1d32 -r 893617de95e5 test-data/closestBedA.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/closestBedA.bed Wed Jun 18 14:33:11 2014 -0400 @@ -0,0 +1,1 @@ +chr1 100 200 diff -r d0df760f1d32 -r 893617de95e5 test-data/closestBedB.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/closestBedB.bed Wed Jun 18 14:33:11 2014 -0400 @@ -0,0 +1,2 @@ +chr1 500 1000 +chr1 1300 2000 diff -r d0df760f1d32 -r 893617de95e5 test-data/groupbyBed.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/groupbyBed.bed Wed Jun 18 14:33:11 2014 -0400 @@ -0,0 +1,14 @@ +chr21 9719758 9729320 variant1 chr21 9719768 9721892 ALR/Alpha 1004 + +chr21 9719758 9729320 variant1 chr21 9721905 9725582 ALR/Alpha 1010 + +chr21 9719758 9729320 variant1 chr21 9725582 9725977 L1PA3 3288 + +chr21 9719758 9729320 variant1 chr21 9726021 9729309 ALR/Alpha 1051 + +chr21 9729310 9757478 variant2 chr21 9729320 9729809 L1PA3 3897 - +chr21 9729310 9757478 variant2 chr21 9729809 9730866 L1P1 8367 + +chr21 9729310 9757478 variant2 chr21 9730866 9734026 ALR/Alpha 1036 - +chr21 9729310 9757478 variant2 chr21 9734037 9757471 ALR/Alpha 1182 - +chr21 9795588 9796685 variant3 chr21 9795589 9795713 (GAATG)n 308 + +chr21 9795588 9796685 variant3 chr21 9795736 9795894 (GAATG)n 683 + +chr21 9795588 9796685 variant3 chr21 9795911 9796007 (GAATG)n 345 + +chr21 9795588 9796685 variant3 chr21 9796028 9796187 (GAATG)n 756 + +chr21 9795588 9796685 variant3 chr21 9796202 9796615 (GAATG)n 891 + +chr21 9795588 9796685 variant3 chr21 9796637 9796824 (GAATG)n 621 + diff -r d0df760f1d32 -r 893617de95e5 test-data/groupbyinput.bed --- a/test-data/groupbyinput.bed Wed Jun 18 13:13:51 2014 -0400 +++ b/test-data/groupbyinput.bed Wed Jun 18 14:33:11 2014 -0400 @@ -1,14 +1,14 @@ -chr21 9719758 9729320 variant1 chr21 9719768 9721892 ALR/Alpha 1004 + -chr21 9719758 9729320 variant1 chr21 9721905 9725582 ALR/Alpha 1010 + -chr21 9719758 9729320 variant1 chr21 9725582 9725977 L1PA3 3288 + -chr21 9719758 9729320 variant1 chr21 9726021 9729309 ALR/Alpha 1051 + -chr21 9729310 9757478 variant2 chr21 9729320 9729809 L1PA3 3897 - -chr21 9729310 9757478 variant2 chr21 9729809 9730866 L1P1 8367 + -chr21 9729310 9757478 variant2 chr21 9730866 9734026 ALR/Alpha 1036 - -chr21 9729310 9757478 variant2 chr21 9734037 9757471 ALR/Alpha 1182 - -chr21 9795588 9796685 variant3 chr21 9795589 9795713 (GAATG)n 308 + -chr21 9795588 9796685 variant3 chr21 9795736 9795894 (GAATG)n 683 + -chr21 9795588 9796685 variant3 chr21 9795911 9796007 (GAATG)n 345 + -chr21 9795588 9796685 variant3 chr21 9796028 9796187 (GAATG)n 756 + -chr21 9795588 9796685 variant3 chr21 9796202 9796615 (GAATG)n 891 + -chr21 9795588 9796685 variant3 chr21 9796637 9796824 (GAATG)n 621 + +chr21 9719758 9729320 variant1 chr21 9719768 9721892 ALR/Alpha 1004 + +chr21 9719758 9729320 variant1 chr21 9721905 9725582 ALR/Alpha 1010 + +chr21 9719758 9729320 variant1 chr21 9725582 9725977 L1PA3 3288 + +chr21 9719758 9729320 variant1 chr21 9726021 9729309 ALR/Alpha 1051 + +chr21 9729310 9757478 variant2 chr21 9729320 9729809 L1PA3 3897 - +chr21 9729310 9757478 variant2 chr21 9729809 9730866 L1P1 8367 + +chr21 9729310 9757478 variant2 chr21 9730866 9734026 ALR/Alpha 1036 - +chr21 9729310 9757478 variant2 chr21 9734037 9757471 ALR/Alpha 1182 - +chr21 9795588 9796685 variant3 chr21 9795589 9795713 (GAATG)n 308 + +chr21 9795588 9796685 variant3 chr21 9795736 9795894 (GAATG)n 683 + +chr21 9795588 9796685 variant3 chr21 9795911 9796007 (GAATG)n 345 + +chr21 9795588 9796685 variant3 chr21 9796028 9796187 (GAATG)n 756 + +chr21 9795588 9796685 variant3 chr21 9796202 9796615 (GAATG)n 891 + +chr21 9795588 9796685 variant3 chr21 9796637 9796824 (GAATG)n 621 + diff -r d0df760f1d32 -r 893617de95e5 test-data/mapBedA.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/mapBedA.bed Wed Jun 18 14:33:11 2014 -0400 @@ -0,0 +1,4 @@ +chr1 10 20 a1 1 + +chr1 50 60 a2 2 - +chr1 80 90 a3 3 - + diff -r d0df760f1d32 -r 893617de95e5 test-data/mapBedB.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/mapBedB.bed Wed Jun 18 14:33:11 2014 -0400 @@ -0,0 +1,5 @@ +chr1 12 14 b1 2 + +chr1 13 15 b2 5 - +chr1 16 18 b3 5 + +chr1 82 85 b4 2 - +chr1 85 87 b5 3 + diff -r d0df760f1d32 -r 893617de95e5 test-data/shuffleBedA.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/shuffleBedA.bed Wed Jun 18 14:33:11 2014 -0400 @@ -0,0 +1,2 @@ +chr1 0 100 a1 1 + +chr1 0 1000 a2 2 - diff -r d0df760f1d32 -r 893617de95e5 test-data/shuffleBedGenome.genome --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/shuffleBedGenome.genome Wed Jun 18 14:33:11 2014 -0400 @@ -0,0 +1,5 @@ +chr1 10000 +chr2 8000 +chr3 5000 +chr4 2000 + diff -r d0df760f1d32 -r 893617de95e5 test-data/t --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/t Wed Jun 18 14:33:11 2014 -0400 @@ -0,0 +1,4 @@ +chr1 100 200 nasty 1 - +chr2 500 1000 ugly 2 + +chr3 1000 5000 big 3 - +