Hi, I was trying to use bedtools to align a bed file with a Fasta file.
The first line of bed look like this:
DN38.contig00001 654 2138
When I was checking the output of bedtools, the result sequence starts at 655, instead of 654.
I am wondering that is the result generated by bedtools always between the 'start' and 'end', instead of starting at the 'start' and ending at 'end'?
Thanks!
1 answer
This is because bedtools uses base-0 for start and base-1 for the end... that's not confusing at all.
BED starts are zero-based and BED ends are one-based.
bedtools users are sometimes confused by the way the start and end of BED features are represented. Specifically, bedtools uses the UCSC Genome Browser’s internal database convention of making the start position 0-based and the end position 1-based: (http://genome.ucsc.edu/FAQ/FAQtracks#tracks1) In other words, bedtools interprets the “start” column as being 1 basepair higher than what is represented in the file.
https://bedtools.readthedocs.io/en/latest/content/overview.html
I didn't know this! Thanks for asking :)
Log in to answer this question.
-
How to make bam from start pos of read 1 and end pos of read 2 from paired end bam files
written by genomics_student •
-
calculate nucleotides frequencies from mapped reads in a BAM file
written by xiaoleiusc
-
Draw coverage plot from bed file
written by ben.kunfang •
-
Getting different results from bedtools and samtools
written by wstla27 •
-
How to get gene annotation for a list of genes
written by wstla27 •
-
bedtools slop set start to -1
written by Medhat
-
Features for detection of gene.
written by anu014
-
Increasing widths of the peaks in a bed file from ChIP-seq
written by morovatunc
-
Augustus result interpretation
written by bioinformaticssrm2011
-
What Does Thickstart (Col 7) Or Thickend (Col 8) Mean In A Bed File?
written by Jordan
there is no tool to "align a bed with a fasta". which sub tool are you using ? what was the cmd line ? did you read Cheat Sheet For One-Based Vs Zero-Based Coordinate Systems ? and How To Ask Good Questions On Technical And Scientific Forums ?