This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Find a position in a fragment

Hi! I'm new in the bioinformatics world and I would like to analyse a data but I need help! Specifically, I want to know how if I can find the fragments or reads that contains a specific position in a file like this format:

chr start end
1       1       5
1       3       6
1       5       10
1       10      100

F.e: If I want to know how many fragments contains the position number 2... How can I do this in bash? THANK YOU VERY VERY MUCH IN ADVANCED!

gene sequencing read

there is no need to shout ! ;)

I'm thinking this can be solved by using one of the bedtools (bedIntersect?)

Thank you very much (without shouting) I will take a look to bedIntersect :)

0 answers

No answers yet.

Log in to answer this question.