This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Generating a Bed file from a Fasta file

I created many random sequences by using NullSeq 2.0, and saved those sequences in a fasta format. Now I want to get the genomic coordinates and make a bedfile. Any ideas on how to do that?

bed nullseq
• 1,111 views • link
updated by Ram Ram • written by pirku • pirku

you could run blast or any other aligner to get the coordinate(s) of your fasta on a reference genome.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

I tried running BLAT, but the output looks like this. Seems like none of the sequences were mapped/aligned

match   mis-    rep.    N's Q gap   Q gap   T gap   T gap   strand  Q           Q       Q       Q   T           T       T       T   block   blockSizes  qStarts  tStarts
        match   match       count   bases   count   bases           name        size    start   end name        size    start   end count
---------------------------------------------------------------------------------------------------------------------------------------------------------------
• 0 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by pirku • pirku

there is not enough information. you have to show us the command you used, a snapshot of your fasta...

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

0 answers

No answers yet.

Log in to answer this question.

More posts like this