Mercurial > repos > fubar > jbrowse2
view maf2bed.pl @ 47:3e53204c2419 draft
planemo upload for repository https://github.com/usegalaxy-eu/temporary-tools/tree/master/jbrowse2 commit 3a43e9e0ffce0966101203102e769d1ced28618a-dirty
author | fubar |
---|---|
date | Wed, 06 Mar 2024 00:23:02 +0000 |
parents | 2de9f585505b |
children | e7a6f7a7148d |
line wrap: on
line source
#!/usr/bin/env perl # https://unpkg.com/browse/jbrowse-plugin-mafviewer@1.0.6/dist/ # license is Apache2_license.txt included here use warnings; use strict; $, = ' '; $\ = "\n"; $, = "\t"; my $id = 0; my $buffer = ''; my $start = 0; my $end = 0; my $score = 0; my $chrom = ''; while (<STDIN>) { chomp; next if /^$/; my @line = split('\s+'); if (/^s\s+$ARGV[0]/) { $chrom = $line[1]; $chrom =~ s/$ARGV[0]\.//; $start = $line[2]; $end = $line[2] + $line[3]; s/^s //; s/ +/:/g; my $temp = $_; $buffer = $buffer eq '' ? $temp : "$buffer,$temp"; } elsif (/^a/) { $score = +(s/^a score=//); if($id > 0) { print $chrom, $start, $end, "$ARGV[0]_$id", $score, $buffer; } $id += 1; $buffer = ''; } elsif (/^s/) { s/^s //; s/ +/:/g; my $temp = $_; $buffer = $buffer eq '' ? $temp : "$buffer,$temp"; } } print $chrom, $start, $end, "$ARGV[0]_$id", $score, $buffer;