We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
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
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
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
htsjdk.samtools.SAMException: Input file must be bam file, not sam file
Is there a way to solve this?
The text was updated successfully, but these errors were encountered:
No branches or pull requests
Hello,
I am trying to run mutserve using .cram files as from documentation seems it supports it:
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?
The text was updated successfully, but these errors were encountered: