if (($seq =~/\w{2}A\w{6}T\w{2}[ACT]\w[GCA{4}T{4}]{4}[AT]/) && ($gcper >= 30 && $gcper <= 52)) { print "$seq\t$seqpos\t$gcper\n"; } #### GCAGGTGGATCTATTTCAT 3201-3220 42.11 TAAGAGGTGTTATTTGGAA 3268-3287 31.58 ATACGATGCTTCAAGAGAA 3346-3365 36.84 CAAGCTCATCATACTGGCT 1201-1220 47.37 GGTACTGACTTTGCTTGCT 2923-2942 47.37 CGTAGTGTTAAGTTATAGT 3003-3022 31.58 GTATGGGTAGGGTAAATCA 3248-3267 42.11 CCTGCTGTGATACGATGCT 3337-3356 52.63 CCTGCGCGCGCGCGATGCT 3300-3318 50.63