How to extract conserved regions from multi sequence alignment by linux terminal?
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
1 answer
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.
Log in to answer this question.
More posts like this
-
How do N's in DNA sequences affect multi-sequence alignments for PCR assay testing?
written by Ryan •
-
How to design one PCR primer for 5 different sequences by linux command line ?
written by Ahmad •
-
extract multiples regions from IGV
written by Morris_Chair
-
Extract N amino acids from fasta file
written by martha.chapa.mc18 •
-
Clipping primer after alignment
written by finswimmer
-
how to identify degenerate primers in reads
written by genya35 •
-
Designing primers to detect different fungi species from a complex soil sample
written by ser •
-
Aligning gene sequences
written by Anand Rao
-
Tool For Recognize Conserved Structural Motifs Among Dna Sequences
written by Anima Mundi
-
How To Merge Overlapping Two Or More Dna Sequence?
written by Geparada
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?
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
http://trimal.cgenomics.org/introduction
Thanks, but
I need command line not GUI
It does provide command-line interface: http://trimal.cgenomics.org/use_of_the_command_line_trimal_v1.2
Ok, thank you I will try