This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Find cooccurences of factors across multiple files

Hey,

I have several hundred files that look like this:

a.file
#genomic positions
1, 3, 4, 9, 10
#factors
a, b, d, g

Each file holds multiple genomic positions (numbers) and factors (characters). I want to investigate, whether there are genomic positions that frequently co-occur with factors across all files, but I have no idea on how to approach this. Can someone point me into the right direction? Is there a tool or a library that might help? Thank you!

chip-seq sequencing gene
• 868 views • link
updated by jordi.planells jordi.planells • written by tirichl • tirichl

1 answer

bedtools intersect accepts multiple file to be intersected. Have you tried with it? You can report the number of occurrences with -c flag.

bedtools intersect -c -a your_file -b factor1 factor2 factorN

Then you could print the lines with more than X occurrences with awk.

awk 'BEGIN{FS="\t";OFS="\t"}{if($4 > X) print $0}'

Hope it helps!

• 0 views • link
written by jordi.planells jordi.planells

Log in to answer this question.

More posts like this