lift-over hg19 to GRch38
1
1
Entering edit mode
8.5 years ago
kanwarjag ★ 1.2k

I am trying to convert 450 K methylation coordinates (Single C location) of given hg19 to GRCH38 following suggestions in previous posts (GRCh37/38(NCBI) vs hg19/hg38(UCSC) and Converting Genome Coordinates From One Genome Version To Another (Ucsc Liftover, Ncbi Remap, Ensembl Api))

using default parameters which are-

Minimum ratio of bases that must remap:0.95:

BED 4 to BED 6 OptionsAllow multiple output regions:off Minimum hit size in query:0;

Minimum chain size in target:0 BED 12 OptionsMin ratio of alignment blocks or exons that must map:1.00If thickStart/thickEnd is not mapped, use the closest mapped base:off

However I get an error :

Deleted in new:
    Sequence intersects no chains
Partially deleted in new:
    Sequence insufficiently intersects one chain
Split in new:
    Sequence insufficiently intersects multiple chains
Duplicated in new:
    Sequence sufficiently intersects multiple chains
Boundary problem:
    Missing start or end base in an exon

Here is one line of input;

chr4    54974682    54974682    +

Any suggestion what is wrong?

Thanks

liftover error • 6.3k views
ADD COMMENT
0
Entering edit mode
8.5 years ago

the error is "deleted in new"

chr4:54974682-54974682 is part of a genomic contig on hg19 that is not anymore mapped on GRch38 because the new assembly is now better built without including this contig.

ADD COMMENT
0
Entering edit mode

I was having around 1Kb region that is strange that all is deleted or modified.

ADD REPLY

Login before adding your answer.

Traffic: 3094 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6