0362494184981193a406895a488aee5e92ece803 markd Tue Feb 7 14:24:25 2023 -0800 support for converting HPRC GFF3 files from Ensembl. Not perfect, however, all genes get converted, so this can be fixed in post-preocessing diff --git src/hg/utils/gff3ToGenePred/tests/makefile src/hg/utils/gff3ToGenePred/tests/makefile index 30b9b15..c31607c 100644 --- src/hg/utils/gff3ToGenePred/tests/makefile +++ src/hg/utils/gff3ToGenePred/tests/makefile @@ -1,30 +1,30 @@ kentSrc = ../../../.. include ../../../../inc/common.mk gff3ToGenePred = ${DESTBINDIR}/gff3ToGenePred # sh commands to deal with gff3ToGenePred shouldhave failed and didn't cmdShouldFail = echo "Error: command should have failed" >&2; false all:: test:: geneMRnaTest noGeneMRnaTest geneMRnaHonorTest discontinuousTest multCdsOutOfExonTest \ noIdTest errCases1Test bogusQuotesTest noExonsTest geneTranscriptTest transcriptCdsParentTest \ minimalGenesTest geneDefaultStatusUnknownTest useNameTest nameAttrIdTest nameAttrNameTest \ frameShiftTest mm10GencodeTest ncbiSegmentsTest ncbiProblemsTest makeBadTest \ - transcriptOnlyTest + transcriptOnlyTest hprcTest geneMRnaTest: mkout ${gff3ToGenePred} input/geneMRna.gff3 output/$@.gp diff expected/$@.gp output/$@.gp noGeneMRnaTest: mkout ${gff3ToGenePred} input/noGeneMRna.gff3 output/$@.gp diff expected/$@.gp output/$@.gp geneMRnaHonorTest: mkout ${gff3ToGenePred} -honorStartStopCodons input/geneMRna.gff3 output/$@.gp diff expected/$@.gp output/$@.gp geneDefaultStatusUnknownTest: mkout ${gff3ToGenePred} -defaultCdsStatusToUnknown input/geneMRna.gff3 output/$@.gp @@ -114,21 +114,25 @@ diff expected/$@.out output/$@.out # -warnAndContinue on maker file with missing gene/mrna records. This test this option on # GFF3 parse tests, not just genePred conversion makeBadTest: mkout ${gff3ToGenePred} -warnAndContinue -attrsOut=output/$@.attrs input/makerBad.gff3 output/$@.gp > output/$@.out 2>&1 diff expected/$@.gp output/$@.gp diff expected/$@.attrs output/$@.attrs diff expected/$@.out output/$@.out # pseudogene annotations as transcripts-only. Yet another non-standard variation transcriptOnlyTest: mkout ${gff3ToGenePred} -allowMinimalGenes input/transcriptOnly.gff3 output/$@.gp diff expected/$@.gp output/$@.gp +# check features used by Ensembl HPRC annotations +hprcTest: mkout + ${gff3ToGenePred} -geneNameAttr=Name -rnaNameAttr=transcript_id input/hprc.gff3 output/$@.gp + diff expected/$@.gp output/$@.gp mkout: @mkdir -p output clean:: rm -rf output