Questions tagged [multiple-sequence-alignment]
The multiple-sequence-alignment tag has no usage guidance.
29
questions
0
votes
1
answer
67
views
How to highlight the specific peptide sequences after performing multiple sequence alignment for the fasta file?
I have the peptide sequences and fasta files separately. I first aligned the fasta files using msa package. After that I'm trying to highlight the peptide sequences ...
3
votes
3
answers
187
views
Read Clustal file in Python
This question was also asked on StackOverflow
I have a multiple sequence alignment (MSA) file derived from mafft in clustal format which I want to import into Python and save into a PDF file. I need ...
1
vote
1
answer
81
views
How to identify sequence origin in DNA shuffle reads?
The data
I have a large number of reads from sequences that were generated by randomly shuffling regions of two parent sequences together. See the following image:
The regions that are shuffled are ...
- 13
1
vote
0
answers
63
views
difficult in runing the R code for ggtree
I tried to run the below code but ended having issue. I have checked thoroughly and I didn't find any mismatch in fasta sequence
...
- 31
2
votes
1
answer
133
views
drawing multiple sequence alignment as tree in R
I am very new here. Please help me with my codes as am trying to draw multiple sequence alignment as tree suing ggtree package in R and i have following error.
...
- 31
3
votes
2
answers
149
views
Select representetive sequence within MSA (protein)
I am trying to pick the best representative sequence from an amino acid MSA. I do not want to use the consensus, as it generates an artificial one. Here is my approach:
I calculate a Percentage ...
- 31
6
votes
4
answers
611
views
multi-sequence alignment of samples with multiple contigs each
I want to perform a multi-sequence alignment on 12 samples that clustered based on cgMLST. Ultimate goal is to find out whether they differ by the presence of a few genes.
I have performed multi-...
- 283
1
vote
2
answers
71
views
Generate consensus protein sequence from relatively gappy alignment?
I would like to know what would be the best way to generate a full-length consensus sequence?
I want to obtain a single representative sequence from a relatively gappy multiple sequence alignment of ...
- 2,294
2
votes
1
answer
264
views
Sequence alignment with MUSCLE
For protein multiple sequence alignment (MSA), I am using MUSCLE (v5.1). The default setting without any specifications gives a 'bad' output even though my proteins are very similar (UniProt sequence ...
- 95
-1
votes
2
answers
104
views
ClustalOmega alignment
I am performing multiple sequence alignment using Clustal Omega for 600 RNA covid-sequences, i.e. cDNA on Genbank, with input characters around 30574 characters for each sequence.
I am running it on ...
2
votes
2
answers
100
views
How to make MSA for contact prediction?
I'd like to try a contact prediction method (like TripletRes, CCMpred etc) using a custom set of sequences. There are about 1.5M sequences, all very similar to the canonical sequence; their ...
- 121
1
vote
4
answers
871
views
Remove gaps from alignment?
I have an MSA (protein sequence) and I have used various programmes (Clustal, Aliview, MEGA11 etc) to align. However, in all programmes I get many gaps which is not ideal as I am trying to construct a ...
- 95
0
votes
1
answer
34
views
find conservative regions relative to a specific accession
I have 7 RefSeq accessions which I performed multiple sequence aligment on it using CLUSTALO....
- 111
0
votes
2
answers
77
views
Use biopython to align SeqRecords stored in dict
I'd like to perform multiple alignments, where a gene from each sample was read in from fasta files. The fasta file represented one sample and had multiple genes. I have read in each sample fasta file ...
- 131
1
vote
2
answers
123
views
Multiple Sequence Alignment
I applied three multiple sequence alignment methods of muscle, mafft and clistal omega on 30 sequences and 400 sequences . In the case of 30 seq, the output of alignment from the 3 methods was in the ...
1
vote
1
answer
116
views
Generating sequences from HMM
I encountered this question in a textbook (source: https://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.721.2540&rep=rep1&type=pdf) and was trying to workout the solution for it. Was ...
0
votes
0
answers
33
views
Fastest way to count where two bam files are homozygote reference (inverted variant calling)
Maybe I'm having a brain-freeze so excuse me if this is a waste of everyone's time...
I am working with an organism with ~700Mb genome.
I have bam files of two individuals that are whole genome ...
- 121
2
votes
0
answers
31
views
Alignment with inserts and keeping the indexing of ref seq intact
Parts of sequences are given below-
Reference sequence (pre-alignment):
ATTAAAGGTTTATACCTTCCCAGGTAACAAACCAACCAACTTTCGATCTCTTGTAGATCTGTTCTCTAAACGAACTTTAAAATCTGTGT
...
- 145
4
votes
0
answers
146
views
How to get phylogenetic tree from multiple genes?
I constructed a phylogenetic tree using a gene (example - secA). I had to gather the same gene sequence for all the required species from public database-NCBI and then constructed the tree after ...
- 73
0
votes
0
answers
147
views
Splitting .bam files into separate samples for tophat2
I have RNA seq data of 12 samples (paired end, so 24 files) and I ran tophat2 for alignment on all of them which gave me one bam file for all the samples. Now I need separate bam files for all the ...
0
votes
2
answers
373
views
Remove from Multi-FASTA by Sequence ID
I want to remove sequence VRE32514 – it doesn’t belong and thus is the reason it lacks additional metadata. However I tried implementing this code from a similar question:
...
- 63
1
vote
1
answer
35
views
BEAUti not accepting taxa tip dates
I'm using BEAUti as part of my pipeline and the sequence names include | to separate the metadata fields - with culture date being at the very end.
When I attempt ...
- 63
1
vote
2
answers
153
views
Getting exact mismatch data from multiple sequence alignment
TL;DR Can anyone point me in the direction of a tool that can be used to input different viral substrain reference sequences (fasta) as well as sequence data (fastq/bam/vcf) and find the most likely ...
- 21
2
votes
1
answer
455
views
How do I plot clusters based on DNA sequence alignment?
I have a group of about 1050 sequences that have been aligned using blastn. I am now trying to create a cluster plot of the sequences in order to identify similar sequences, something like this:
...
- 23
1
vote
1
answer
88
views
File format of substitution matrix in clustalw
I need to set the substitution matrix used by command line CLUSTALW when comparing DNA sequences to:
0 -1 -1 -1
-1 0 -1 -1
-1 -1 0 -1
-1 -1 -1 0
from my ...
- 159
1
vote
1
answer
123
views
Median string problem & multiple sequence alignment
I read about the median string problem as an introduction to the multiple sequence alignment, however none of the MSA algorithms used seems to be using the idea of finding the median string.
To my ...
- 159
2
votes
0
answers
60
views
Why do I obtain different output results with blast vs awk commands
I have an awk command that identifies 30 pb from two multifasta files.
When I used two input files:
E.g. 100 sequences each, I get the same result with the ...
- 41
2
votes
2
answers
70
views
Multiple Alignment cost application
Is there any biological application/case where someone would be interested in the estimated total cost of the alignment between a set of sequences (genes or amino acids) without the aligned sequence ...
- 159
1
vote
0
answers
29
views
Call NCBI Multiple Sequence Alignment Viewer with target URL?
I am storing a sequence aligment file on a webserver and now would like to visualize the results. I used the NCBI Multiple Sequence Alignment Viewer . The issue is that there are several manual steps ...
- 1,285