Is it okay to map reads from WES designed for GRCh37 to GRCh38?
Suppose I have WES reads (in fastq format), and the capture regions are based on GRCh37 coordinate. Can I map the WES reads using bwa to GRCh38 coordinate? If not, what problems would arise?
• 1,079 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
samuelandjw •
•
written
by
samuelandjw •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Can I Use SureSelect v7 on Exome Data Captured with SureSelect v6?
written by George •
-
Should we trust genotypes called in simple tandem repeat regions?
written by samuelandjw •
-
VCF from GRCH37 to GRCH38: how to manage reference differences
written by jpuntomarcos •
-
GATK site filtering: when should I use VQSR or hard filtering?
written by samuelandjw •
-
Novel regions in GRCh38 compared to GRCh37
written by sf21 •
-
BWA mapping parameters for DNA capture sequencing protocol
written by graeme.thorn
-
Weird BWA MEM Behavior: Zero-Length Read
written by aostern •
-
MYH2 gene mutation location
written by Monos •
-
Genomic coordinates of Exome Capture (SureSelect_Human_All_Exon_V4) in GRCh38
written by gil.hornung
-
SNP coordinate systems in UCSC LiftOver (hg19), NCBI Remap (GRCh37(hg19)) and dbSNP (GRCh38(hg38))
written by curious •
there won't be any problem; The reads don't "know" where they've been designed. The only thing to do is to liftover your bed capture. May be you'll lost a few regions in the process.