Mercurial > repos > dcouvin > pmlst2
diff pmlst.pl @ 0:140d4f9e1f20 draft default tip
Uploaded
author | dcouvin |
---|---|
date | Mon, 06 Sep 2021 16:00:46 +0000 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/pmlst.pl Mon Sep 06 16:00:46 2021 +0000 @@ -0,0 +1,29 @@ +#!/usr/bin/perl -w +use strict; +use warnings; + +my $STpmlst = ""; + +my $path = $ARGV[0]; +my $seq = $ARGV[1]; +my $scheme = $ARGV[2]; + + +open PMLST, "python3 $path/pmlst/pmlst.py -i $seq -s $scheme -p $path/pmlst/pmlst_db/ |"; +while (<PMLST>) { + chomp ($_); + if ($_ =~ "sequence_type") { + #print "ST:$_ \n"; + my @table = split (/:/,$_); + $STpmlst = $table [1]; + #chomp ($STpmlst); $STpmlst = chop($STpmlst); + $STpmlst =~ s/\s+//; + $STpmlst = substr($STpmlst,1); + my @table2 = split (/'/,$STpmlst); + $STpmlst = $table2[0]; + #print "$STpmlst\n" ; + } +} + +#print "Sequence\tSequence Type (pMLST)\n" ; +print "pMLST for $seq:\t$STpmlst\n" ;