Off topic:Blast Short Sequences Containing Repeats
1
0
Entering edit mode
11.7 years ago
grafabian • 0

I try to detect repeats using blastn. I use the same sequence as query and subject. The sequence contains a 21bp tandem duplication. I expect three hits:

  1. The whole sequence aligns itself.
  2. First occurance of duplication in query aligns second occurance in subject.
  3. Second occurance of duplication in query aligns first occurance in subject.

This works very well for Sequence1 in the web version and in my standalone version. But for Sequence2 this is not working. I tried this for many sequences and it is working very well, but for that case it fails. Any suggestions?

Sequence1:

>gtacgtgcattttaaagattttccaatggaaaagaaatgctgcagaaacatttggcacattccattcttaccttggcacattccattcttaccaaactctaaattttctcttggaaa

Sequence2:

>gagtacgtgcattttaaagattttccaatggaaaagaaatgctgcagaaacatttggcacattccattcttaccttggcacattccattcttaccaaactctaaattttctcttggaaa

Blastn command: blastn -query seq.fa -subject seq.fa -word_size 21

blast repeats • 3.2k views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 2182 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6