利用perl从大量数据中提取需要的有用信息
open (YEAST,"f:\\perl\\yeast_set2.S288C.stampy.pe.recal.variants.filter.vcf")||die("can not open the file!");while(<YEAST>){if($_=~/^##/){}elsif($_=~/^#/){@category= split/\s/,$_;open(COUNTBASE,">>f:\\perl\\a.csv")||die("can not open the file!");print COUNTBASE " @category[0] @category[1] @category[3] @category[4] @category[5] \n"}else{@information=split/\s/,$_;$number=8;until($number>22){$number++;if(@information[$number]=~/^0\/0/){open(COUNTBASE,">>f:\\perl\\a.csv")||die("can not open the file!");printf COUNTBASE "@category[$number] @information[0] @information[1] @information[3] @information[4] @information[5] S\n" ; }elsif(@information[9]=~/^0\/1/){open(COUNTBASE,">>f:\\perl\\a.csv")||die("can not open the file!");printf COUNTBASE "@category[$number] @information[0] @information[1] @information[3] @information[4] @information[5] H\n"; }elsif(@information[9]=~/^1\/1/){open(COUNTBASE,">>f:\\perl\\a.csv")||die("can not open the file!");printf COUNTBASE "@category[$number] @information[0] @information[1] @information[3] @information[4] @information[5] Y\n"; }}}}