This is a test version of Biostars. For the public version, visit https://www.biostars.org.
CrossMap lift over bed but not the same vcf file

Trying to see what is happening from the previous post CrossMap liftover failed

I downloaded a vcf file from ftp://ftp.ensemblgenomes.org/pub/plants/release-31/vcf/zea_mays/zea_mays.vcf.gz and try to lift it over using CrossMap also it did not work
But when I convert the vcf file to bed file first using vcf2bed from bedops package and the I lift it over using CrossMap it worked !!

any suggestion here? or explanation for this behavior

snp vcf bed sequence
• 2,533 views • link
written by Medhat Medhat

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • help with CrossMap
    written by raalsuwaidi • raalsuwaidi
    Hello all, I would really appreciate your help as I am new to working with different file builds and having a setback lifting a vcf …
  • Missing SNP in gnomAD hg38 liftover vcfs
    written by Nicolas Rosewick Nicolas Rosewick
    Hi, On gnomAD website there are both hg19 and hg38 vcf. hg38 vcfs are liftover from hg19. https://gnomad.broadinstitute.org/downloads I was analyzing some data using both …
  • STRINGTIE WARNING: no reference transcripts were found for the genomic sequences where reads were m…
    written by hirahameed222 • hirahameed222
    Hi everybody! I am using RNA-Seq protocol for identifying differentially expressed gene in RNASEQ data of 6 wheat varieties. I am following this protocol https://www.nature.com/articles/nprot.2016.095. …
  • LiftOver a VCF file
    written by Susmita Mandal Susmita Mandal
    Hey all, I have been trying to liftover a particular VCF file from GRCm38 to NCBIm37. I have used UCSC LiftOver tool, Ensembl API, CrossMap …
  • "makeTxDbFromGFF" error ?!
    written by dina.hesham139 dina.hesham139
    Hey, I'm using deseq2 after kallisto to analyze rice data. I'm using an ensembl gtf and I want to create a txdb. I used this …
  • bedops vcf2bed core dump
    written by Ram Ram
    Hello again, I'm extracting BED from a 77G VCF file using `bedops` `vcf2bed`, and it does a bunch of `core dump`s. This is because I'm …
  • bedops vcf2bed eats up all allocated memory
    written by Ram Ram
    Hello, I'm trying to run `bedops` `vcf2bed` on a huge (27G) VCF file. I'd like to extract a BED of the deletions in the VCF. …
  • unable to download file from plant ensembl
    written by Medhat Medhat
    I am working on Maize AGV4 release 33, trying to download vep cash from ```wget ftp://ftp.ensemblgenomes.org/pub/plants/release-33/vep/zea*``` It gives error > No such directory ‘pub/plants/release-33/vep I …
  • CrossMap liftover failed
    written by Medhat Medhat
    I am trying to liftover HapMap.vcf dowlonaded from [here](http://cbsusrv04.tc.cornell.edu/users/panzea/download.aspx?filegroupid=15) for maize from AGPv3 to AGPv4 using this command; ``` CrossMap.py vcf AGPv3_to_AGPv4.chain.gz all_chr_hmp321_withDPGL_imputed.vcf.gz genome_AGPv4.fa ./all_chr_hmp321_withDPGL_imputed_AGPv4.vcf …
  • Lift-over on a VCF
    written by win win
    <p>Hi all,</p> <p>I have a VCF file that was generated using HG18 and I want to convert it to VCF but that has HG19 data.</p> …