R which

 setwd("E:/courses/molecular biology/homework1st")
 genes <- read.table('genes.txt',sep = '\t',header = T)

 chrID<-c("chr1","chr2","chr3","chr4","chr5","chr6","chr7","chr8","chr9","chr10","chr11","chr12",
          "chr13","chr14","chr15","chr16","chr17","chr18","chr19","chr20","chr21","chr22","chrX","chrY")
 pro_counts=vector()
 for(chr in chrID){
   chr<-genes[genes$chrom==chr,]
   protein_coding<-chr$cdsEnd-chr$cdsStart
   pro_counts<-c(pro_counts,length(which(protein_coding>)))
 }