Question: bed file conversion and comparing with the genome to find its region
1
gravatar for kapil.joshi036
3.4 years ago by
Student ,School of life sciences, Manipal University, Manipal, India
kapil.joshi03680 wrote:

I have a bed file with the co-ordinates includes chr start and end columns.

#chr    start   end
chrM    1        1
chrM    2        2
chrM    3        3
chrM    4        4
chrM    5        5
chrM    7        7
chrM    8        8
chrM    9        9
chrM    10      10
chrM    11      11
chrM    12      12
chrM    13      13
chrM    14      14

here there are 2 region are present i want to get output like

#chr    start   end
 chrM   1        5
 chrM    7       14

which tool should i use to perform this after this i need to find its location into the mitochorial genome

bed genome • 1.1k views
ADD COMMENTlink modified 3.4 years ago by michael.ante3.4k • written 3.4 years ago by kapil.joshi03680
2
gravatar for Pierre Lindenbaum
3.4 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum123k wrote:

I'm afraid your BED file only describes deletions because a BED in half-open : see chromEnd in https://genome.ucsc.edu/FAQ/FAQformat.html#format1

the tool you need is "bedtools merge" http://bedtools.readthedocs.io/en/latest/content/tools/merge.html

which tool should i use to perform this after this i need to find its location into the mitochorial genome

uh ???

ADD COMMENTlink modified 3.4 years ago • written 3.4 years ago by Pierre Lindenbaum123k

sir , actually i jst have the co-ordiantes of mitocondrial sequence so i wannt to find its location into the mitochondrial genome

ADD REPLYlink written 3.4 years ago by kapil.joshi03680

If you have co-ordinates then aren't they the location in the genome?

ADD REPLYlink written 3.4 years ago by genomax73k

i found the solution for it download the mtDNA gene bed file and annotate with bedtools

ADD REPLYlink written 3.4 years ago by kapil.joshi03680
1
gravatar for michael.ante
3.4 years ago by
michael.ante3.4k
Austria/Vienna
michael.ante3.4k wrote:

Hi,

To merge a bed-file, just use bedtools merge.

bedtools merge -i my_file.bed > merged_file.bed

[EDIT] regarding your second question, you may want to have a look at bedtools annotate or bedtools closest.

Cheers,

Michael

ADD COMMENTlink modified 3.4 years ago • written 3.4 years ago by michael.ante3.4k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1643 users visited in the last hour