This is a test version of Biostars. For the public version, visit https://www.biostars.org.
What the last output of liftover means?

from pyliftover import LiftOver

lo = LiftOver('hg17', 'hg18')

lo.convert_coordinate('chr1', 1000000)

[('chr1', 949796, '+', 21057807908)]

What '+' and 21057807908 mean?

liftover python
• 1,669 views • link
updated by Wayne Wayne • written by ManuelDB ManuelDB

2 answers

This is the comment from author's script posted on github.

  Note that coordinates are 0-based, and even at negative strand are relative to the beginning of the genome.
    I.e. position 0 strand + is the first position of the genome. Position 0 strand - is also the first position of the genome 
    (and the last position of reverse-complemented genome).
• 0 views • link
written by supertech supertech

https://github.com/konstantint/pyliftover/blob/master/pyliftover/liftover.py#L66

Returns a list of possible conversions for a given chromosome position. The list may be empty (no conversion), have a single element (unique conversion), or several elements (position mapped to several chains). The list contains tuples (target_chromosome, target_position, target_strand, conversion_chain_score), where conversion_chain_score is the "alignment score" field specified at the chain used to perform conversion. If there are several possible conversions, they are sorted by decreasing conversion_chain_score.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this