in reply to Find a sequence in a multifasta files and motifs
#! /usr/bin/perl use warnings; use strict; use feature qw{ say }; my $seq; while (<DATA>) { if (my $line_number = (/^>JJ22$/ ... /^>/)) { chomp; $seq .= $_ unless $line_number == 1 || $line_number =~ /E0/; } } my $count = () = $seq =~ /[MT].[KA]/g; say $count; __DATA__ >JJ57 MKIKLVTVGDAKEEYLIQGINEYLKRLNSYAKRETIEVPDEKAPEKLSDAEMLQVKEKEGEYVFVLAI NGKQLSSEEFSKEIFQTGISGKSNLTFTTCFSLGLSDSVLQRIMKGEPYHKL >JJ22 MDQNGASGSHPNRASTRKGAHARERGATVSAMSANRSNIIDEMAKICEADRQTFAIARRTRNESQ FFGFRTASNKAIEITEAMEKRGAMFLTQSKATDQLNGWQPSDEPDKTSAESEPWFRGKQLSSEEFS KEIFQTGISGKSNLTFTTCFSLGLSDSVLQRIMKGEPYHKL >JJ41 MWKTVAPIFAAIFAVGLCGTFRTNTRKGEPTTKCFVFVHDTKARIYQCTFKTWSCPWLNNIVSAQF QFVTGANYKIVVKLVGELFTETALFNWSSPTTIFTGLGTLITADKTLDCDSNML
The first check in the unless part excludes the header, the E0 is appended to the last line number, so it excludes the next header.
|
|---|
| Replies are listed 'Best First'. | |
|---|---|
|
Re^2: Find a sequence in a multifasta files and motifs
by rebkirl (Acolyte) on Jun 28, 2019 at 10:12 UTC |