［求助］编程求助──有一定挑战性的问题

CODE:

#! /usr/bin/perl

use warnings;
use strict;

while(<DATA>)
{
        if(/^>/)
        {
                print;
                next;
        }
        chomp;
        while(/A(?=AAA)/g)
        {
                if($+[0] < 10)
                {
                        print substr($_, 0, $+[0] + 13) . "\n";
                }
                else
                {
                        print substr($_, $+[0] - 11,  24) . "\n";
                }
        }
}


__END__
>Y64G10A.7|Y64G10A.7|14302238|14319530|IV
ATATTTCAATTTCAAAATTAAATTTATTTTTATTTTCAACAATTCAATCTCTTGCCCAAATCTTTCTTCCCCTCATTTCTCATGTAAGAATCTCATTTTTTGACAGTTAGCCTAGTTATGTTACAGTTGATTAATTTAAATACTTTTACTGTTTCAGTACTAAAATTTATAAAAATCCGAGTGTCTTTTTCTTAATGAAATCAGTTTCAATTATTCCAAGCCATAATAACTTTCCATGTTTTCACACAGGTCTATGTGTTTCATTAAATAAATTT
>Y53F4B.31|gst-28|15166113|15167310|II
ATCGAAAATTTTCTTGAAATTTTGTTCAATGAAGTAATATAATTACAAATTCGAATAAAT
>ZK792.8|atg-4.2|11668005|11671460|IV
TTATTGACAATTTTTGAAATTCTGTACACTGAACATTTATTGGTGCTTGATATGTATTGATGACGCAATGTGAACATTATTAAATAAAA
>Y69H2.8|nhr-241|18699416|18700957|V
AGTCATGTTGCCAGGCACAGGCAAATTGTTT
>Y71F9AL.13|rpl-1|2875396|2876313|I
GGTAAAATGTCGAAGGTTTCCCGCGAGTCTCTTAATGAGGCTATCGCTGAGGTCCTTAAGGGATCCTCGGAAAAGCCACGCAAGTTCCGTGAGACTATCGAGCTCCAGATCGGTCTCAAGAACTACGACCCACAGAAGGACAAGCGTTTCAGCGGATCCATCAGGTAGCTACGACCTGTGAAAACAGGGAGTTGGGACTTTCAAAGCCTAACGGACCCCCACCCAGGGTCTTAAATTTTTGGGGAGGCAATTATTGCCGCTTTGAAATAATGATATGAGCTGGTATAATACTGAAGCACATTCCACGCCCAAACATGAAGGTGTGCGTCTTCGGAGACCAGCACCATCTTGACGAGGCTGCCGCCGGAGACATTCCATCGATGAGCGCCGATGACTTGAAGAAGCTTAACAAGCAGAAGAAGCTCATCAAGAAGCTCGCCAAGAGCTACGATGCTTTCATCGCTTCCGAATCCTTGATCAAGCAGATCCCACGTATCCTCGGTCCAGGACTCAACAAGGCCGGAAAATTCCCTTCCGTCGTTACCCACGGAGAATCTCTCCAAAGCAAGAGTGACGAGATCCGCGCCACCGTCAAGTTCCAGATGAAGAAG
>Y47G6A.4|Y47G6A.4|3547093|3552007|I
ACGGAAACCCCTCACACATTCATTTCATGGTATTAAATATAAGCATTCTAAAATCCTAATTGTACTAATAGTATTATACATGATATACTCTGTCTGTCTGCTTGGTACAAATTTACTTCATTCTCAATTTTCAAAACAATAAATCATATATTCTTT
>Y17G9B.4|Y17G9B.4|4752110|4753154|IV
AATAAATGTGTTTTCTTATTTAAAATATTATTCTGTTTCGAG
>Y65B4BR.5|Y65B4BR.5|535792|536676|I
CTTGTTTCCTGATGACCTTGCAGATACTCTTGT
>Y47D3B.11|Y47D3B.11|11456758|11462763|III
AAAATTAGCCAGAATACTCTATCAAAAAATATCATATCACACCTTTTCTTCTCCCAGTTTTTCTCTTTTTGTCGTGCCTTGCTCGTCATCTCTCCAACCCCGATCTCTGTATTTTTCTCCACGTATGACTGCGCCTTGTAGGCACGCAGGTAGGCATTTTTGTGCCTACGTGGATTAATTGCCTAAATTGTCTTAAAATGCTTAGTGTTTTTCAGGTGTCAAACTTTCCTCCCCATGTGTACCCAACGGGTCAATCACATTCTCTTCTTGCATATGAATACATTTTCATTTTTCTGTTGTTTTTTAGTTTTTTTTTAGACTACCATTTTTTTTTTAACTTTCTAGAAGTTTCTAGTAATAA

QUOTE:

原帖由 ly5066113 于 2008-11-12 22:52 发表
1. perl -ne 'print if /^>/ || /.{50}/' urfile

QUOTE:

原帖由 ly5066113 于 2008-11-12 23:10 发表
2.
#! /usr/bin/perl

use warnings;
use strict;

while()
{
        if(/^>/)
        {
                print;
                next;
        }
        chomp;
        while(/A(?=AAA)/g)
        {
                if($+[0] < 10)
                {
                        print substr($_, 0, $+[0] + 13) ...