Error running example with hg19 genome about deepvariant HOT 6 CLOSED

Are you saying you've swapped hg19.fa.gz for our hs37d5.fa.gz fasta? That won't work because the chromosome names are different - 20 in hs37d5 and it sounds like chr20 in hg19.fa.gz. DeepVariant requires the BAM file and the reference genome to have consistent contigs.

from deepvariant.

I do not think this is the case. The example BAM file, NA12878_S1.chr20.10_10p1mb.bam, contains reads mapped to chromosome 20 which show chr20 in RNAME for all the reads. This means the fasta file should contain chromosomes named as 'chr'.

Morever, the example fasta, ucsc.hg19.chr20.unittest.fasta, comes from hg19 and it works. Any idea why using the entire genome would make the run fail?

from deepvariant.

The only real requirement is that the fasta and the BAM need to have a large number of consistent contigs (meaning same length in bp and with the same name) and DeepVariant will process the dataset. If you do samtools view -H on the BAM and cat the hg19.fa.gz.fai file, do these have the same set of contigs with the same lengths?

from deepvariant.

Ok, might be doing something wrong. How do you create your .gz.fai file? Is it coming from indexing a fasta file compressed with bgzip or how?

from deepvariant.

I already have a copy of hs37d5.fa on my machine, so to create a gzipped, indexed version I do:

cp hs37d5.fa /tmp
bgzip /tmp/hs37d5.fa
samtools faidx /tmp/hs37d5.fa.gz
head /tmp/hs37d5.fa.gz.fai
1       249250621       52      60      61
2       243199373       253404903       60      61
3       198022430       500657651       60      61
4       191154276       701980507       60      61
5       180915260       896320740       60      61
6       171115067       1080251307      60      61
7       159138663       1254218344      60      61
8       146364022       1416009371      60      61
9       141213431       1564812846      60      61
10      135534747       1708379889      60      61

from deepvariant.

Ok, I think that solves the issue. When indexing over the gz only indexed the first chromosome for some reason in my machine

from deepvariant.