2 answers
I ended up using samtools and bcftools like so:
samtools mpileup -uf human_g1k_v37.fasta sample.bam --region MT | bcftools call -m --ploidy 1 > sample_mt_only.vcf
This puts the entire MT in vcf. If I would supply "-v" parameter to bcftools, it would only store the variants.
samtools --version
samtools 1.7
bcftools --version
bcftools 1.7
We published a caller for mitochondrial ancient DNA but it is suitable for mitochondrial modern DNA as well. For mitochondrial modern DNA the problem is even simpler. The name of the program is endoCaller part of the (schmutzi package)
You can call it as such:
endoCaller -seq consensus.fa -log consensus.log -name [name sample here] chrMT.fa inputMT.bam
You might want to use the option -qual, this will filter bases according to a PHRED scale. e.g. -qual 50 means the confidence will be 1/100,000.
I have editted my post with "mitochondrial modern" and "mitochondrial ancient" to make it clear that I mean mitochondrial. Apologies for the confusion!
Log in to answer this question.
-
Best way to deal with chimeric reads?
written by deniselavezzari •
-
Fasta sequence from Bam file
written by rse
-
what does the result show in Hisat ?
written by XBria
-
Search for reads aligned to multiple chromosomes
written by oleglovky •
-
Entire Sequence from Bam File using Python Pysam Module
written by BioICoder •
-
Split bam by multiple chromosome to a single bam file
written by l0o0
-
Easy way to calculate NQS (Neighborhood Quality Standard) Score?
written by jik664 •
-
Tool for identifying mitochondrial SNP from exome seq data
written by peris
-
using GATK UnifiedGenotyper on a single chromosome
written by stefano.iantorno •
-
Wgs Using Solid And Illumina Platform
written by Ashutosh Pandey
Sure, tell the variant caller to only consider that chromosome. For more information, please read the manual of the variant caller. Most variant callers have suh options.