mutserve icon indicating copy to clipboard operation
mutserve copied to clipboard

Mutserve call does not work with CRAM file

Open dana-obk opened this issue 2 years ago • 0 comments

Hello,

I am trying to run mutserve using .cram files as from documentation seems it supports it:

Usage: mutserve call [--baq] [--deletions] [--help] [--insertions] [--no-ansi]
                     [--no-freq] [--version] [--write-fasta] [--write-raw]
                     [--alignQ=<alignQ>] [--baseQ=<baseQ>]
                     [--contig-name=<contig>] [--level=<level>] [--mapQ=<mapQ>]
                     [--mode=<mode>] --output=<output> --reference=<reference>
                     [--threads=<threads>] [<input>...]
Call homoplasmic and heteroplasmic positions.
      **[<input>...]          BAM/CRAM files**
      --alignQ=<alignQ>     Minimum Align Quality
                              Default: 30
      --baq                 Enable BAQ
                              Default: false
      --baseQ=<baseQ>       Minimum Base Quality
                              Default: 20
      --contig-name=<contig>
                            Specifify mtDNA contig name
                              Default: null
      --deletions           Call deletions (beta)
                              Default: false
      --help
      --insertions          Call insertions (beta)
                              Default: false
      --level=<level>       Minimum Heteroplasmy Level
      --mapQ=<mapQ>         Minimum Map Quality
                              Default: 20
      --mode=<mode>         Specifify mutserve mode
                              Default: mtdna
      --no-ansi             Disable ANSI support
                              Default: false
      --no-freq             Use Frequency File
                              Default: false
      --output=<output>     "Output (txt or vcf)
      --reference=<reference>
                            Reference
      --threads=<threads>   Number of threads
      --version
      --write-fasta         Write fasta file
                              Default: false
      --write-raw           Write raw file
                              Default: false

This is my code: mutserve call --reference rCRS.fa --output out_from_cram.vcf --threads 2 my_input.cram

But the program does not recognise the .cram file, and sees it as a sam file: htsjdk.samtools.SAMException: Input file must be bam file, not sam file

Is there a way to solve this?

dana-obk avatar May 06 '22 11:05 dana-obk