Question: how to get input file of closest-features a bed-format file
0
gravatar for MEITUO
9 weeks ago by
MEITUO10
MEITUO10 wrote:

HI everyone!

i am trying to calculate the gene distance of assembly contigs by closest-features. but this tool requires a bed format input file. how can i get this bed-format file?

sequencing gene genome • 130 views
ADD COMMENTlink modified 9 weeks ago by jordi.planells330 • written 9 weeks ago by MEITUO10

What kind of data do you have? Could you please share more details?

ADD REPLYlink written 9 weeks ago by jordi.planells330

metagenic data, the contigs after assembly

ADD REPLYlink written 8 weeks ago by MEITUO10

if you have the bam format you can use bedtools bamtobed (see here).

Then you can use bedtools closest (see here) to intersect the bed file obtained earlier with the features (also in bed). If you specify the -d flag, it will report the distance as well.

ADD REPLYlink written 8 weeks ago by jordi.planells330

thanks! and is the bam file that the qc reads mapping back to assembly contig? do this need filter low mapped reads or just can use it directly after bwa or bowite??

ADD REPLYlink written 8 weeks ago by MEITUO10

That I didn't get it sorry... Please explain better what you mean. If you have doubts about what a bam file is, you should spend some time familiarizing with the file types used in bioinformatics. See here for a short description of the most common file types in bioinformatics.

ADD REPLYlink written 8 weeks ago by jordi.planells330
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: 1510 users visited in the last hour
_