#!/usr/bin/perl -w open (AN_ANNOT, "<$ARGV[0]")or die "can not open vcf file\n"; open (OUTDGV,">$ARGV[0].mcCorrolcnv.txt") || die "can not open $ARGV[0].dgvcnv.txt ($!)\n"; @Annotation_file = ; close AN_ANNOT; foreach $annotation_line(@Annotation_file) { if ( $annotation_line =~ /^chr/ ) { chomp $annotation_line; my( $chr, $start, $end, $description ) = split( /\s+/ , $annotation_line ); if ( $description =~ /dgv_Cnv=True,overlap=\d+\.\d+\%\;/ ) { print OUTDGV "$chr\t$start\t$end\t$&\n"; } } }