# HG changeset patch # User nml # Date 1507736785 14400 # Node ID 1b3c339fd390730e899dd3f20ba709d093c79ed6 planemo upload for repository https://github.com/phac-nml/galaxy_tools commit a0204b99a722240fe9b03b78a0786b30aa8ecc96 diff -r 000000000000 -r 1b3c339fd390 fasta_to_bed.pl --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/fasta_to_bed.pl Wed Oct 11 11:46:25 2017 -0400 @@ -0,0 +1,25 @@ +#!/usr/bin/env perl +use strict; +use warnings; +use Bio::SeqIO; +use Bio::Seq; +use Getopt::Long; + +my ($fasta_file, $out); + +GetOptions( + "i|input=s" => \$fasta_file, #contigs fasta file + "o|output=s" => \$out +); + +my $seqio_object= Bio::SeqIO->new (-format =>'fasta', -file=>$fasta_file); +open(my $out_fh, ">", $out) || die "Could write to file '$out'\n"; + +while (my $seq_object = $seqio_object->next_seq()) { + my $seq_id = $seq_object->display_id(); + my $length = $seq_object->length(); + print $out_fh $seq_id . "\t" . "1" . "\t" . $length . "\t" . $seq_id . "\n"; + +} + + diff -r 000000000000 -r 1b3c339fd390 fasta_to_bed.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/fasta_to_bed.xml Wed Oct 11 11:46:25 2017 -0400 @@ -0,0 +1,28 @@ + + Create Bed file from Multiple fasta file + + perl-bioperl + perl-getopt-long + + + + + + + + + + + + + + + + Create a bed file format using the start and end position of one or more fata record in a fasta file. + + + diff -r 000000000000 -r 1b3c339fd390 test-data/contigs.fasta --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/contigs.fasta Wed Oct 11 11:46:25 2017 -0400 @@ -0,0 +1,40 @@ +>lcl|AF319579.1_cds_AAN14420.1_1 [gene=dnaA] [protein=DnaA] [protein_id=AAN14420.1] [location=1090..2430] [gbkey=CDS] +TTGGAAAATATCTCTGATTTATGGAATAGTGCCTTAAAAGAATTAGAAAAAAAGGTAAGCAAGCCAAGTT +ATGAAACCTGGTTAAAATCAACAACTGCTCATAACTTAAAAAAAGATGTATTAACTATTACAGCTCCAAA +TGAATTCGCTCGTGACTGGCTAGAATCTCATTATTCGGAACTAATCTCCGAAACACTTTATGATTTAACA +GGGGCAAAACTAGCTATTCGTTTTATCATTCCCCAAAGTCAGGCAGAAGAGGATATAGATCTCCCTTCTG +TTAAGCAAAAACACGCACATGATGAGTCTAATCACTTACCACAGAGTATGTTAAACCCAAAATATACATT +TGATACATTTGTTATTGGCTCTGGTAATCGTTTTGCACATGCTGCTTCTTTAGCTGTAGCTGAAGCACCA +GCTAAAGCCTATAATCCGCTCTTTATTTACGGGGGCGTTGGACTTGGAAAAACACATTTAATGCACGCAA +TTGGGCATTATGTAATTGAACACAACCCAAATGCAAAAGTTGTATATTTATCATCTGAAAAATTTACAAA +TGAATTTATTAATTCAATTCGTGATAATAAAGCTGTCGATTTTCGTAATAAATATCGTAATGTAGATGTC +TTATTGATAGATGATATTCAATTTCTTGCTGGAAAAGAACAAACACAGGAAGAGTTTTTCCATACATTTA +ATGCGTTACATGAAGAAAGTAAACAAATTGTGATTTCAAGCGATAGACCACCAAAAGAAATCCCAACTTT +AGAAGATCGCCTTCGTTCTCGCTTTGAATGGGGGTTAATTACGGACATTACGCCACCAGATTTAGAAACC +AGAATCGCAATTTTACGTAAAAAAGCAAAAGCAGAAGGTCTTGATATACCAAATGAAGTTATGCTTTATA +TCGCAAATCAAATCGATTCAAATATTCGTGAACTAGAGGGTGCACTTATTCGAGTTGTAGCTTATTCATC +TTTAATTAATAAAGATATGAATGCTGATTTAGCGGCTGAAGCGCTTAAAAATATTATTCCAAACTCTATA +CCTAGAATTATTTCTATTTCTGATATTCAAAAGGCTGTTGGGGGCGTCTATCAAGTAAAATTAGAAGATT +TCAAAGCAAAAAAACGGACAAAATCCGTTGCATTCCCTCGCCAAATTGCGATGTATTTATCACGTGAACT +TACGGACTCATCTTTGCCAAAAATAGGTGAAGAATTTGGTGGCCGCGATCATACAACAGTAATCCATGCT +CACGAAAAAATTTCGAAGCTATTAAAAACGGATACTCAATTACAAAAACACGTTGAAGAAGTTAAGGATA +TTTTAAAATAG +>lcl|AF319579.1_cds_AAN14421.1_2 [gene=dnaN] [protein=DNA polymerase III beta chain] [protein_id=AAN14421.1] [location=2613..3758] [gbkey=CDS] +ATGCGTTTTACAATACAAAAAGACTATCTTGTAAGAAGTGTACAAGATGTAATGAAGGCAGTTTCTTCTC +GTACAACAATTCCAATTCTTACAGGGATTAAAGTTGTAGCAACTGAAGAAGGAGTTACTTTAACAGGTAG +TGATGCAGACATCTCTATCGAATCATTTATCCCAGTCGAAGATGCTGGAAAAGAAATTGTGGAAATTGAA +CAAACAGGAAGTATTATTTTACAAGCAAAATATTTTAGTGAGATCGTAAAAAAATTACCTAAAGAAACTG +TCGAAATTTCTGTAGAAAATCATTTTATGACAAAAATAAAATCTGGAAAATCAGAGTTCAACTTAAATGG +TTTAGATGCTGCTGAATATCCGTTATTACCACAAATCGAAGAACATCATGTATTTAAGATTCCAACAGAT +TTACTAAAGCACATGATTCGTCAAACAGTATTTGCAGTTTCAAGTTCTGAAACAAGACCGATCTTGACAG +GTGTAAACTGGAAGGTATATAACAGCGAGCTAACTTGTATTGCAACAGATAGTCACAGACTAGCACTTCG +TAAAGCAAAAATCGAAGGGTATAATATTGCCGATGAATTTCAAGCGAATGTCGTTATTCCTGGTAAGAGC +TTAAGTGAATTAAGTAAAATTCTAGATGAATCTGAAGAAATGGTAGATATCGTTATTACGGAGTATCAAG +TATTATTCCGTACAAAACATTTATTATTCTTCTCAAGATTATTAGAAGGAAATTATCCAGATACAACGCG +TTTAATTCCAGCTGAAAGTAAAACTGATATATTTGTAAATACAAAAGAATTTTTACAAGCGATTGATCGT +GCGTCGCTGTTAGCAAGAGATGGTCGTAACAACGTTGTTAAATTATCGACATTAGAACAACAGATGTTAG +AGATTTCTTCGAATGCACCAGAAATCGGAAAAGTAGTAGAAGAGGTTCAATGTGAAAATGTAGATGGAGA +AGAATTAAAAATATCTTTTAGTGCAAAATATATGATGGACGCATTAAAGGCTTTAGACAGTACAGAAATT +AAAGTTAGTTTTACTGGAGCGATGAGACCGTTTTTAATTCGCACAGTAAATGATAATTCCATTATCCAAT +TAATTTTACCAGTCCGTACTTACTAA + diff -r 000000000000 -r 1b3c339fd390 test-data/results.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/results.bed Wed Oct 11 11:46:25 2017 -0400 @@ -0,0 +1,2 @@ +lcl|AF319579.1_cds_AAN14420.1_1 1 1341 lcl|AF319579.1_cds_AAN14420.1_1 +lcl|AF319579.1_cds_AAN14421.1_2 1 1146 lcl|AF319579.1_cds_AAN14421.1_2