func TestReadGTFWithPadding(t *testing.T) {
	gtfRecords := ReadGTF(context.Background(), testutil.GetFilePath(
		"//go/src/github.com/grailbio/bio/fusion/parsegencode/testdata/annotation.gtf"), false, 20, false,
		0)
	assert.EQ(t, len(gtfRecords), 4)
	assert.EQ(t, findGene(gtfRecords, "ENSG1.1").transcripts["ENST1.1"].exons[0].start, 80)
	assert.EQ(t, findGene(gtfRecords, "ENSG1.1").transcripts["ENST1.1"].exons[0].stop, 170)
	assert.EQ(t, findGene(gtfRecords, "ENSG2.1").transcripts["ENST3.1"].exons[0].start, 130)
	assert.EQ(t, findGene(gtfRecords, "ENSG2.1").transcripts["ENST3.1"].exons[0].stop, 320)
}
TestReadGTFWithoutPadding will test gtf reading functionality with 20bp padded exons
package parsegencode
import (
	"context"
	"testing"
	"github.com/grailbio/testutil"
	"github.com/grailbio/testutil/assert"
	"github.com/grailbio/testutil/expect"
)
func findGene(genes []*GencodeGene, name string) *GencodeGene {
	for _, gene := range genes {
		if gene.geneID == name {
			return gene
		}
	}
	panic(name)
}
func TestReadGTFWithoutPadding(t *testing.T) {
	gtfRecords := ReadGTF(context.Background(), testutil.GetFilePath(
		"//go/src/github.com/grailbio/bio/fusion/parsegencode/testdata/annotation.gtf"), false, 0, false,
		0)
	assert.EQ(t, len(gtfRecords), 4)
	expect.EQ(t, findGene(gtfRecords, "ENSG1.1").transcripts["ENST1.1"].havanaTranscript, "OTTHUMT1.1")
	expect.EQ(t, findGene(gtfRecords, "ENSG4.1").transcripts["ENST6.1"].exons[1].stop, 300)
	expect.EQ(t, findGene(gtfRecords, "ENSG3.1").strand, "-")
}