Samtools depth with 0-based BED file
Hi everyone,
I have a question about using samtools with a BED file as input for coordinates. I used samtools depth with a zero-based BED file and I get one-based coordinates as output. For example, for the range chr1 9 10 in my BED file, I get the position chr1 10 as output from samtools. I assume that I would have to subtract 1 from the SAM positions to reformat to zero-based coordinates.
Am I doing this correctly?
Best
Lynn
1 answer
Yes, samtools output is 1-based. See e.g. https://github.com/samtools/samtools/issues/1063
Log in to answer this question.
More posts like this
-
0-based coordinates from samtools depth
written by marco.barr •
-
How to get MD tag sequence of a read from a sam file given a bed file?
written by n.a.tekkey •
-
per sample mean coverage and standard deviation
written by ijlal.hyder2012 •
-
What are "unused reference sequences" in a bam file?
written by anniejeanforster •
-
How does bedtools bamtobed calculate end position
written by timing •
-
Obtaining depth for all positions for different exonic ranges from mpileup file
written by kspata
-
Creating alternate Bed file by incorporating INDELs locations
written by ersgupta •
-
Problems In Running Samtools Mpileup/Depth For Specific Regions
written by C Shao •
-
Samtools Depth Can Only Have Stdout?
written by Chris •
-
Bed Coordinates
written by Florianino •
