Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: CRAM reference registry and the GRch38 reference genome
0
gravatar for Mehulsharma.253
17 months ago by
Mehulsharma.253 • 10
Mehulsharma.253 • 10 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
grch38 variantcalling samtools bam cram • 858 views
ADD COMMENT • link •
modified 17 months ago by h.mon ♦ 29k • written 17 months ago by Mehulsharma.253 • 10
2
gravatar for h.mon
17 months ago by
h.mon ♦ 29k
Brazil
h.mon ♦ 29k 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 COMMENT • link written 17 months ago by h.mon ♦ 29k
Please log in to add an answer.

Similar posts • Search »

Content
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