It would be helpful to have several test cases with actual input and actual matching output for each case.
#!/usr/bin/perl # http://perlmonks.org/?node_id=1188439 use strict; use warnings; $_ = join('', grep /^[\sACGT]+$/, <DATA>) =~ tr/ACGT//cdr; my $count = 0; while( /(?=(.{21}GG))/g and $count < 1000 ) { /(?=$1).+$1/ or print ">crispr_@{[++$count]}\n$1\n"; } __DATA__ > >2L type=chromosome_arm; loc=2L:1..23011544; ID=2L; dbxref=REFSEQ: +NT_033779,GB:AE014134; MD5=bfdfb99d39fa5174dae1e2ecd8a231cd; length=2 +3011544; release=r5.54; species=Dmel; CCTCTCATTTTCTCTCCCATATTATAGGGAGAAATATGATCGCGTATGCG AGAGTAGTGCCAACATATTGTGCTCTTTGATTTTTTGGCAACCCAAAATG TAAATTCATTGCAACGTTAAATACAGCACAATATATGATCGCGTATGCGA GAGTAGTGCCAACATATTGTGCTAATGAGTGCCTCTCGTTCTCTGTCTTA TATTACCGCAAACCCAAAAAGACAATACACGACAGAGAGAGAGAGCAGCG GAGATATTTAGATTGCCTATTAAATATGATCGCGTATGCGAGAGTAGTGC CAACATATTGTGCTCTCTATATAATGACTGCCTCTCATTCTGTCTTATTT TACCGCAAACCCAAATCGACAATGCACGACAGAGGAAGCAGAACAGATAT
In reply to Re: program to look for specific K-mer sequence
by tybalt89
in thread program to look for specific K-mer sequence
by pearllearner315
| For: | Use: | ||
| & | & | ||
| < | < | ||
| > | > | ||
| [ | [ | ||
| ] | ] |