@@ -0,0 +1,1213 @@
+#!/usr/bin/env python3%0A%22%22%22Investigate GAF reading error on saccharomyces%22%22%22%0A%0Afrom os import system%0Afrom os.path import join%0Afrom os.path import exists%0Afrom os.path import basename%0Aimport wget%0A%0Afrom tests.utils import REPO%0Afrom goatools.anno.gaf_reader import GafReader%0A%0A%0Adef test_i195():%0A    %22%22%22Investigate GAF reading error on saccharomyces%22%22%22%0A    fin_gaf1 = join(REPO, 'sgd.gaf')%0A    dnld_gaf1 = 'http://current.geneontology.org/annotations/sgd.gaf.gz'%0A    _dnld_gaf(fin_gaf1, dnld_gaf1)%0A%0A    fin_gaf2 = join(REPO, 'gene_association.sgd.gaf')%0A    dnld_gaf2 = 'http://downloads.yeastgenome.org/curation/literature/gene_association.sgd.gaf.gz'%0A    _dnld_gaf(fin_gaf2, dnld_gaf2)%0A%0A    # Read GAF%0A    print('READING: %7BGAF%7D'.format(GAF=basename(fin_gaf1)))%0A    ##DVK objanno_sgd1 = GafReader(fin_gaf1)%0A    objanno_sgd2 = GafReader(fin_gaf2)%0A    ##DVK assert objanno_sgd1 == objanno_sgd2%0A%0A%0Adef _dnld_gaf(fin_gaf, dnld_gaf):%0A    %22%22%22Download saccharomyces GAF%22%22%22%0A    if not exists(fin_gaf):%0A        system('rm -f %7BGAF%7D*'.format(GAF=fin_gaf))%0A        wget.download(dnld_gaf)%0A        fin_gaz = '%7BGAF%7D.gz'.format(GAF=fin_gaf)%0A        system('gunzip %7BZIP%7D'.format(ZIP=fin_gaz))%0A%0A%0Aif __name__ == '__main__':%0A    test_i195()%0A

Add test for https://github.com/tanghaibao/goatools/issues/195