Mercurial > repos > dcouvin > pmlst
diff pmlst.pl @ 0:6dcb66689a7d draft default tip
Uploaded
author | dcouvin |
---|---|
date | Wed, 25 Aug 2021 21:12:34 +0000 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/pmlst.pl Wed Aug 25 21:12:34 2021 +0000 @@ -0,0 +1,28 @@ +#!/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 Type (pMLST)\t$STpmlst\n" ;