Find patterns in DNA sequence
I want te be able to count the number of occurrences in a given sequence (for example ACTTTAG) in the GRCh38 reference genome. Is there an existing tool for doing this? Thanks!
• 3,232 views
•
link
4 answers
I wrote a simple script for finding patterns (regular expressions in fact) in fasta files, it's fastaRegexFinder.py and I also happen to mention it in this post Quadruplex sequence batch prediction
If you just want to count the number of occurrences you can do
fastaRegexFinder.py -f genome.fa -r 'ACTTTAG' | wc -l
• 0 views
•
link
You can also use bowtie1.
It is specially nice in finding (mapping) short sequences..
• 0 views
•
link
Log in to answer this question.