Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

0.3.0 does not support hg38 aligned BAM file? #10

Open
nswh opened this issue May 27, 2021 · 1 comment
Open

0.3.0 does not support hg38 aligned BAM file? #10

nswh opened this issue May 27, 2021 · 1 comment

Comments

@nswh
Copy link

nswh commented May 27, 2021

Tested on 0.3.0 0.2.1 0.1.3
Only 0.1.3 passed
However, 0.1.3 was developed based on hg19 annotation?

$ docker run --rm -it -v $(pwd):/home drmjc/mity:0.3.0 call --prefix NA12878-2403 NA12878-2403.bam
mity version 0.3.0
Calling mitochondrial variants
Running FreeBayes in sensitive mode
FreeBayes failed: b"unable to find FASTA index entry for 'chrM'\n"

$ docker run --rm -it -v $(pwd):/home drmjc/mity:0.2.1 call --prefix NA12878-2403 NA12878-2403.bam
mity version 0.2.1
Calling mitochondrial variants
Traceback (most recent call last):
  File "/usr/local/bin/mity", line 29, in <module>
    args.func(args)
  File "/usr/local/lib/python3.7/site-packages/mitylib/commands.py", line 45, in _cmd_call
    genome = select_reference_genome(args.reference, None)
  File "/usr/local/lib/python3.7/site-packages/mitylib/util.py", line 305, in select_reference_genome
    assert len(res) == 1
AssertionError

$ docker run --rm -it -v $(pwd):/home drmjc/mity:0.1.3 call --prefix NA12878-2403 NA12878-2403.bam
mity version 0.1.3
Calling mitochondrial variants
Running FreeBayes in sensitive mode

BAM file was aligned using Broad GATK reference genome bundle https://console.cloud.google.com/storage/browser/genomics-public-data/resources/broad/hg38/v0/

@drmjc
Copy link
Member

drmjc commented May 31, 2021

Hi,
Did you try specifying the reference genome?

docker run --rm -it -v $(pwd):/home drmjc/mity:0.3.0 call --prefix NA12878-2403 --reference hg38 NA12878-2403.bam

PS --normalise is recommended

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants