Depth of Coverage in Mitochondrial Genome with Samtools
Hi all
I want to ask you how to calculate mitochondrial genome's depth of coverage with samtools on linux ubuntu? I have done this work about some regions of chromosomes on BAM file with samtools, but those commands don't answer about mitochondrial genome.
Thanks
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to calculate depth of coverage and applying read depth filter
written by Sd •
-
Can I use vg tools to calculate the depth for a bed range file like with samtools depth?
written by carolynzy •
-
Calculating coverage/depth for nanopore dRNAseq reads with reference transcriptome/genome
written by kenneditodd •
-
Samtools coverage vs samtools depth for % coverage
written by melissachua90 •
-
no output after running samtools depth
written by Ana
-
Basic SAMtools tutorial
written by statfa
-
coverage per base using samtools depth and mpileup
written by tonja.r
-
samtools idxstats versus samtools view command
written by frida.danielsson •
-
Depth of coverage for specific target from bam file
written by jfjiang •
-
Where can I find the coordinates of mitochondrial DNA on the human genome on UCSC Genome Browser?
written by quachtina96 •
what does that mean ?
also check you're using the correct chromosome name (chrM instead of chrMT...)
I meant that those commands were not executed.
Yes, thank you. I tried this "chrM", but it gave me an error. In command to calculate the depth of coverage over the entire length of the mitochondrial genome, what regions should I type in the command after "chrM"?
?
what is the command ?
The commands that I used to calculate the depth of coverage in different regions of different chromosomes did not work for the mitochondrial genome.
Is there a special command to calculate the depth of coverage of certain regions in mitochondrial genome?
that's what I asked : what is the command ?
The command was: samtools depth -a -H -r chr number:from-to -o FILE restored BAM file
Now my problem about mitochondrial genome is solved with this command: samtools depth -a -H -r MT:from-to -o restored BAM file
Thanks a lot for your help
You may check first if the BAM file contains contigo corresponding to chrM/MT with or without the chr prefix. Like: samtools view -H bam_file | grep SQ
I did it. thank you for your help Now my problem about mitochondrial genome is solved with this command: samtools depth -a -H -r MT:from-to -o restored BAM file