Question: Bowtie alignment of piRNAs to find potential targets
1
gravatar for mujupas
4 months ago by
mujupas50
Japan
mujupas50 wrote:

Hi, I have a set of piRNAs in fasta format and I want to find their potential targets. According to literature nucleotides in position 1-10 in the piRNAs sequences should constitute the seed and do not tolerate any mismatch while nt 11-20 tolerate up to 3 mismatches. Bowtie can specify the length of the seed, but I'm having troubles finding the right combination of parameters to translate the piRNA target pairing logic into a single bowtie command line.

Is there any solution to get this done with a single alignment run or do I have to run a first alignment guided by the seed 1-10nt, and successively another alignment for the 11-20 region allowing for up to 3 mismatches?

Thank you in advance for your help!

bowtie pirna alignment • 234 views
ADD COMMENTlink modified 4 months ago • written 4 months ago by mujupas50
1
gravatar for Charles Plessy
4 months ago by
Charles Plessy2.6k
Japan
Charles Plessy2.6k wrote:

How about running the alignment tolerating 3 mismatches anywhere, and then filtering the results using the CIGAR string ? (You would need an aligner that produces CIGAR strings distinguishing matches (=) and mismatches (X); I do not remember if bowtie does).

ADD COMMENTlink written 4 months ago by Charles Plessy2.6k
1
gravatar for Charles Plessy
4 months ago by
Charles Plessy2.6k
Japan
Charles Plessy2.6k wrote:

Or maybe you can use LAST with a custom seeding scheme ?

ADD COMMENTlink written 4 months ago by Charles Plessy2.6k
0
gravatar for mujupas
4 months ago by
mujupas50
Japan
mujupas50 wrote:

Thank you, I'll try both options!

ADD COMMENTlink written 4 months ago by mujupas50
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1317 users visited in the last hour