Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • running Grch37 to grch38 liftover tool locally

    Hi all,

    I am trying to run the UCSC ligtover tool on my local machine (linux version), but am running into some difficulties.

    my input file is a standard bed:

    Code:
    chr1    17369   17391   miR-6859-3p
    and to run the tool I am doing this (as specified by their usage):

    Code:
    ./liftOver -errorHelp my_hg19.bed hg38ToHg19.over.chain my_new.bed unMapped
    I keep getting this error:

    Code:
    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

    Yet if I run this directly on the USCS website, it runs just fine.

    I must be doing something wrong, but can't quite figure it out.....

  • #2
    were do they specify this usage?

    you simply need to remove "-errorHelp"
    This just prints the help for error messages (which is your ouptput), everything after it is ignored

    Comment


    • #3
      What direction are you trying to go with liftover? Title of your post says Grch37 to Grch38 but you seem to be using the chain file for going the other direction.

      Comment


      • #4
        Yes, I am trying to go from Grch37 to Grch38. I can't seem to find the chain file to go the other (and didn't realize the I had the wrong one).

        The specified usage is:
        Code:
        liftOver - Move annotations from one assembly to another
        usage:
           liftOver oldFile map.chain newFile unMapped
        oldFile and newFile are in bed format by default, but can be in GFF and
        maybe eventually others with the appropriate flags below.
        The map.chain file has the old genome as the target and the new genome
        as the query.
        
        ***********************************************************************
        WARNING: liftOver was only designed to work between different
                 assemblies of the same organism. It may not do what you want
                 if you are lifting between different organisms. If there has
                 been a rearrangement in one of the species, the size of the
                 region being mapped may change dramatically after mapping.
        ***********************************************************************
        
        options:
           -minMatch=0.N Minimum ratio of bases that must remap. Default 0.95
           -gff  File is in gff/gtf format.  Note that the gff lines are converted
                 separately.  It would be good to have a separate check after this
                 that the lines that make up a gene model still make a plausible gene
                 after liftOver
           -genePred - File is in genePred format
           -sample - File is in sample format
           -bedPlus=N - File is bed N+ format
           -positions - File is in browser "position" format
           -hasBin - File has bin value (used only with -bedPlus)
           -tab - Separate by tabs rather than space (used only with -bedPlus)
           -pslT - File is in psl format, map target side only
           -ends=N - Lift the first and last N bases of each record and combine the
                     result. This is useful for lifting large regions like BAC end pairs.
           -minBlocks=0.N Minimum ratio of alignment blocks or exons that must map
                          (default 1.00)
           -fudgeThick    (bed 12 or 12+ only) If thickStart/thickEnd is not mapped,
                          use the closest mapped base.  Recommended if using 
                          -minBlocks.
           -multiple               Allow multiple output regions
           -minChainT, -minChainQ  Minimum chain size in target/query, when mapping
                                   to multiple output regions (default 0, 0)
           -minSizeT               deprecated synonym for -minChainT (ENCODE compat.)
           -minSizeQ               Min matching region size in query with -multiple.
           -chainTable             Used with -multiple, format is db.tablename,
                                       to extend chains from net (preserves dups)
           -errorHelp              Explain error messages
        Even removing the 'errorHelp' just gives back this usage help.

        Comment


        • #5
          Originally posted by lre1234 View Post
          Yes, I am trying to go from Grch37 to Grch38. I can't seem to find the chain file to go the other (and didn't realize the I had the wrong one).
          Here is the chain file for hg19toGrch38: http://hgdownload.soe.ucsc.edu/golde....over.chain.gz

          Comment


          • #6
            Thank you. That works!!

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Strategies for Sequencing Challenging Samples
              by seqadmin


              Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
              03-22-2024, 06:39 AM
            • seqadmin
              Techniques and Challenges in Conservation Genomics
              by seqadmin



              The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

              Avian Conservation
              Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
              03-08-2024, 10:41 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, Yesterday, 06:37 PM
            0 responses
            10 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, Yesterday, 06:07 PM
            0 responses
            9 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-22-2024, 10:03 AM
            0 responses
            49 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-21-2024, 07:32 AM
            0 responses
            67 views
            0 likes
            Last Post seqadmin  
            Working...
            X