Question: CRAM reference registry and the GRch38 reference genome
0
gravatar for Mehulsharma.253
17 months ago by
Mehulsharma.25310 wrote:

Hi. I'm downloading some CRAM files from 1000genomes for use in variant calling. I will convert them to BAM (since most tools can't call from CRAM files). I'm a bit confused about how to go about this process since I assumed a straightforward Samtools based conversion based on a given reference FASTA.

What's the deal with the reference registry and why do I require it ?

I went through the README document on the FTP site but I'm still quite confused

I've already downloaded the hg38 FASTA from Human Genome Resources

ADD COMMENTlink modified 17 months ago by h.mon29k • written 17 months ago by Mehulsharma.25310
2
gravatar for h.mon
17 months ago by
h.mon29k
Brazil
h.mon29k wrote:

I've already downloaded the hg38 FASTA from Human Genome Resources

You have to download the files referenced by the checksums found on the cram headers. To decompress the cram files, you need exactly the same reference as used for compression. To ensure the correct reference is used, the 1000genomes cram files contain the identity of the reference contigs used for compression, this identity is given by MD5 or SHA1 checksums.

ADD COMMENTlink written 17 months ago by h.mon29k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2137 users visited in the last hour