gene_num=0;
uncertain_num=0;
for y=1:2
    if y==1
        f='./before2013/';
    else
        f='./after2013/';
    end
    list=dir(f);
    for i=3:1:size(list,1)
        load(['./gbffmat/',list(i).name(1:end-4),'.gbff.mat']);
        cds=g.CDS;
        gene_num=gene_num+size(cds,2);
        for j=1:1:size(cds,2)
            product=cds(j).product;
            product=reshape(product',1,size(product,1)*size(product,2));
            if isempty(product) || contains(product,'hypothetical') || contains(product,'putative')
                uncertain_num=uncertain_num+1;
            end
        end
    end
end

uncertain_num/gene_num