Hello, I wanted to inspect the (non-N) regions that were added in, in GRCh38/hg38 that were not there in GRCh37/hg19. Is there a hg38 bed file with those "new" regions somewhere or is there an easy way to get that without lifting over again? Thanks for your help!
1 answer
I used the excellent remap service from NCBI: http://www.ncbi.nlm.nih.gov/genome/tools/remap
to go from my hg19 assembly to my hg38 assembly.
I made a dummy bed file using the chromosome lengths, which you can get from here:
mysql --user=genome --host=genome-mysql.cse.ucsc.edu -A -e "select chrom, size from hg19.chromInfo" > hg19.genome
and selected:
Minimum ratio of bases that must be remapped: 0
Maximum ratio for difference between source length and target length: 100000
and that gave me all possible mappings.
Log in to answer this question.
-
1000 Genome: Cause of duplicate variants with different genotypes
written by JourneyToAbyss
-
Tutorial: Automated dbSNP lookup by rsID position, plus genome build liftover
written by Kevin Blighe
-
convert genome coordinates from b36 to hg19
written by @anna •
-
What aren't hg19-hg38 liftOver chains symmetrical?
written by d-cameron
-
Remap GRCh38.p13 to GRCh38
written by schmau •
-
Count all (non N) nucleotides within regions defined in BED file
written by William
-
Novel regions in GRCh38 compared to GRCh37
written by sf21 •
-
Using hg38 assembly with hg19 exome capture data
written by jdm204 •
-
DNase data for GRCh38
written by Adam Przedniczek •
-
EnsGene.txt (UCSC/hg19) and Homo_sapiens.GRCh38.76.gtf (GRCH38) positions do not match
written by pwg46
See these posts:
A: How to determine reference genome that was used for alignment given BAM with obs
A: Convert BED(UCSC) to VCF