Skip to content

Result disagrees with ucsc's liftover #10

Open
@carbocation

Description

@carbocation

Example go program:

package main

import (
	"fmt"

	glo "github.com/blachlylab/GLO"
)

func main() {
	liftover := new(glo.LiftOver)
	liftover.Init()

	liftover.LoadChainFile("hg19", "hg38", "hg19ToHg38.over.chain")

	target := glo.NewChainInterval("chr1", 120887250, 120887251)
	response := liftover.Lift("hg19", "hg38", target)
	fmt.Println(response)
}

Result:

  • [chr1:143924375-143924377]

Meanwhile, providing the same range to UCSC's liftover (chr1 120887250 120887251) yields:

  • chr1 143924010 143924011

Metadata

Metadata

Assignees

Labels

No labels
No labels

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions