How to extract conserved regions from multi sequence alignment by linux terminal?
1
0
Entering edit mode
2.8 years ago
Ahmad ▴ 10

Hello everyone, I hope you are well I want to make common PCR primer for group of sequences , So I want to extract conserved regions from multi sequence alignment by linux terminal. Is there a tool or loop that can help me with that?

Thanks for your attention

conserved linux Multisequence Alignment primer regions • 1.7k views
ADD COMMENT
0
Entering edit mode

I realise your question is about finding conserved region but your main task is primer design. So why not use a tool dedicated to designing primers? eg primer3/primerblast?

ADD REPLY
0
Entering edit mode

I think that these tools do not work on the conserved regions or similar regions in the sequences, but here they work individually on each sequence, and do not produce a common primers

ADD REPLY
0
Entering edit mode
ADD REPLY
0
Entering edit mode

Thanks, but

I need command line not GUI

ADD REPLY
0
Entering edit mode

It does provide command-line interface: http://trimal.cgenomics.org/use_of_the_command_line_trimal_v1.2

ADD REPLY
0
Entering edit mode

Ok, thank you I will try

ADD REPLY
0
Entering edit mode
2.8 years ago
Joe 21k

You could relatively simply adapt code like here: Trim sequences based on alignment in python

Instead of searching for the first non-gap character, you could simply substitute in a test that checks if every character in a given column is the same.

ADD COMMENT

Login before adding your answer.

Traffic: 1373 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