GRCH38 and hg19

I did alignments with hg19 but I tried to use GRCH 38 for annotation but I get lot of undocumented genes like

AC010884.1
AC011366.3
AC011747.3
AC013275.2

Can we really use GRCh38 for hg19 or not at all. Should I be using GRCH37?

Thanks

hg19 human-genome

• 10,360 views • link

updated 3.6 years ago by Ram

• written 11.1 years ago by kanwarjag

1 answer

The reference genome used for alignment should be the same reference genome used for further processing, including annotation.

Consider either using the current aligned file with GRCh37 (hg19) for annotation or perform the alignment again using GRCh38 (hg38) and then you can use GRCh38 (hg38) during annotation.

• 0 views • link

updated 3.6 years ago by Ram

• written 11.1 years ago by Irene@Sequencing.com

More posts like this

How to use SNPlocs.Hsapiens.dbSNP150.GRCh38
written 2.0 years ago by learner__lin2003 •

I try to covert chr:pos to snp rsid, and I find that my gwas data was built on GRCh 38 and dbSNP 150. So I …
About hg19 reference genome
written 4.1 years ago by luckysardar171 •

Hello biostars I am trying to get some information from paper its genecount with genename. according to my research of interest I need chromosome information …
Converting Genome Reference Consortium to UCSC, i.e. GRCh* -> hg*
written 4.7 years ago by DavidStreid •

I am wondering how to convert a `GRCh*` file to its corresponding `hg*`. I would like to do this as part of the process to …
biomaRt: Timeout on getBM().
written 5.1 years ago by lincaijin1994 •

Hi. I want to get the exon annotation of a list of chromosome interval (hg19). Here is the code: getBM_value <- list( chromosome_name = bed_df$chromosome_name, …
Why does the iGenomes have decoy for only GRCh38.
written 5.6 years ago by Apprentice

Hi. I would like to use GRCh38 as a reference sequence for NGS data analysis. I tried to get the file of GRCh38 at iGenomes. …
Debug LiftOver Failures
written 5.9 years ago by afollette •

I often have to LiftOver genome coordinates for SNV from GRCh37 to GRCh38. Recently, I found a SNV in TP53 that did not LiftOver from …
Reference genome and associated gene model (gtf, gff3) file for splice junctions in RNA STAR in Gal…
written 5.9 years ago by nash52 •

Hi, I am trying to process RNA-seq data from fastq files to get raw counts using Galaxy. I want to use STAR for alignment purposes. …
VEP cache trancript versions
written 6.1 years ago by Chadi Saad

I have a same variant obtained once on GRCH37 and once on GRCH38 I want to annotate it with VEP to obtain the HGVS annotation. …
What is the difference between GRCh37 and hs37? And hg19?
written 8.0 years ago by juanfdelahoz •

Hi! I've been struggling with the naming conventions of human reference genomes... I know hg19 and GRCh37 are the same, but different names for each …
ncRNA gencode files
written 10.3 years ago by kanwarjag

In gencode how can I download ncRNA for GRCH37 ( (=comprehensive) gtf file I see GRCH38 are there but none of the link provided me …