This is a test version of Biostars. For the public version, visit https://www.biostars.org.
FIMO output to bed

Trying to figure out a good way to covert fimo output files such as the gff file to a standard bed4-6 file for use elsewhere.

Was doing the following but it seems to generate an unusually structured bed file:

sortBed -i fimo.gff | gff2bed > fimo.bed
fimo meme bedops bedtools
• 4,515 views • link
updated by BioinfGuru BioinfGuru • written by rbronste rbronste

My apologies neglected to search back, however your response adds to my understanding of how to do this so thanks.

• 0 views • link
written by rbronste rbronste

1 answer

A typical FIMO output looks like following (NOTE: This is with version 4.11.4, might slightly differ with other version)

#pattern name   sequence name   start   stop    strand  score   p-value q-value matched sequence
Homeodomain.UP00163_1   CLYBL|chr13|100529567|100530067 231     247     -       11.982  4.03e-05                TTCTTTAATTAATACAA
Homeodomain.UP00163_1   CLYBL|chr13|100529567|100530067 232     248     +       14.036  9.32e-06                TGTATTAATTAAAGAAT
Homeodomain.UP00163_1   CLYBL|chr13|100450472|100450972 247     263     +       10.5405 9.06e-05                TAACCTAATTAGATTCT

You can try the following to get a BED file

cat fimo_result.txt | grep -v pattern | cut -f2 | tr '|' '\t' | cut -f2-4 | sort -k1,1V -k2,2n > fimo_to_bed.bed

output:

chr1    202997705       202998205
chr2    207998147       207998647
chr3    140986720       140987220
chr13   100450472       100450972
• 0 views • link
written by venu venu

Unless I;m mistaken, this is describing the txt file...not the gff file.

The gff file needs parsing. to get out the p and q values

Why don't you just use the fimo.txt file?

Just cut the first 2 columns from it. The txt file contains all the information in the gff file but also filters by the p/q value you selected when running fimo

• 0 views • link
written by BioinfGuru BioinfGuru

Your command just pulls out the factor name instead of the interval, seems like the current format for FIMO output is a little different - can you suggest on how to adjust? Thanks!:

 # motif_id motif_alt_id    sequence_name   start   stop    strand  score   p-value q-value matched_sequence
MA0032.2    FOXC1   chr5:29248245-29248604  202 212 +   15.44   4.51e-07    0.18    tatgtaaacat
MA0032.2    FOXC1   chr18:47026418-47027081 24  34  +   15.3    9.99e-07    0.18    TATGTAAATAT
MA0032.2    FOXC1   chr11:17198451-17199123 170 180 +   15.3    9.99e-07    0.18    tatgtaaatat
MA0032.2    FOXC1   chr15:17197994-17198595 230 240 -   15.3    9.99e-07    0.18    TATGTAAATAT
MA0032.2    FOXC1   chr2:24838945-24839358  294 304 -   15.3    9.99e-07    0.18    TATGTAAATAT
• 0 views • link
written by rbronste rbronste

Try this

cat fimo_result.txt | grep -v motif_alt_id | cut -f3 | tr ':-' '\t' | sort -k1,1V -k2,2n > fimo_to_bed.bed
• 0 views • link
written by venu venu

Log in to answer this question.

More posts like this
  • awk for interval midpoints and length
    written by rbronste rbronste
    Kind of a basic awk question but just trying to take the following intervals: chr1 10001 10150 chr1 10008 10090 chr1 10025 10207 chr1 10058 …
  • Limiting interval sizes with awk
    written by rbronste rbronste
    Kind of a basic awk question. I have a bed file where I have been doing the following to limit interval sizes to those under …
  • BEDOPS sort-bed BED_LINE_LEN question
    written by rbronste rbronste
    So I am just trying to sort a bed file in the following way: sort-bed GSE63137_ATAC-seq_PV_neurons_HOMER_peaks.bed > GSE63137_ATAC-seq_PV_neurons_HOMER_peaks_sorted.bed It is only a 3 column bed …
  • FIMO output that preserves padj values
    written by rbronste rbronste
    Hi I am generating a fasta file from a DESeq2 output in the following way: bedtools getfasta -fi /son/genomes/mm10.fa -bed treatment_pw_padj.001_sorted.bed -fo treatment_pw_padj.001_sorted.fa The above …
  • bwtool stalled liftover
    written by rbronste rbronste
    I am trying to liftover some mm9 bigwigs to mm10 with bwtool in the following way: bwtool lift GSE81122_121+122_Amy_CK_120min_1M_H3K27ac.ucsc.bigWig mm9ToMm10.over.chain.gz test1.bigwig In doing so the …
  • Fimo output for classes of TFs from specific database
    written by rbronste rbronste
    Currently if getting a fimo database for specific TFs can do the following: fimo -max-stored-scores 10000000 --thresh 1e-6 -oc /path/to/output/directory -motif /path/to/motif/file \ /path/to/HOCOMOCOv11_HUMAN_mono_meme_format.meme /path/to/whole_genome.fa …
  • Extracting the exact motif location from FIMO
    written by rbronste rbronste
    So I am using FIMO on some small datasets to look for motifs of interest and it works quite well to scan a PWM against …
  • Official gene symbol list to BED file
    written by rbronste rbronste
    Simply trying to generate an annotated bed file for a specific assembly only from a list of official gene symbols, in my case derived from …
  • bedops gff2bed error
    written by varsha619 varsha619
    Hello, I am trying to convert a .gff file to bed format using bedops gff2bed tool. Command used - sortBed -i in.gff | gff2bed > …
  • FIMO GFF output to standard BED
    written by rbronste rbronste
    Hi, Wondering how to convert FIMO output GFF to standard BED. Using one converter it still looks partially like the FIMO GFF output. Giving an …