#!/usr/bin/perl -w open (IN, "$ARGV[0]"); $/=">"; $i=0; while () { if(/^(\S+)?.*?\n(.*)/ms) { $i++; $seq[$i]=$2; $name[$i]=$1; } } close IN; for $j(1..$i){ $j++; tr/>//,$name[$j]; tr/>//,$seq[$j]; @sequence=split(/N+/,$seq[$j]); $n=length(@sequence); $seqs_to_print=$n-1; # print "$name[$j]\tLength of seqs:$n\n"; for $x(0..$seqs_to_print){ print ">$name[$j].$x\n"; print "$sequence[$x]\n"; } }