<p>All that is needed is to align the sequences within a formal aligner like MAFFT or muscle. Then pass that through the PDGSpider converter. Arlequin should just take it from there.</p>
<p>Manual annotation of the alignment is normal e.g. <code>seqoTron</code> and is traditionally regarded as the best approach rather than do this algorithmically. Thus I'm suggesting (strongly suggesting) do this prior the PGDSpider conversion. If an alignment has not been performed Arlequin should throw an error.</p>
<p>If there are trailing gaps (indels or '-') at the beginning and end of the alignment Arlequin will/should ignore these. You can of course simply cut them out with the alignment editor (Mega X will also work for alignment editing). The important point is that Arlequin must correctly identify homologues its an allele frequency package and nothing more. If there are errors in the alignment Arlequin will believe there is population structure when it could be absent.</p>
<p>I can't honestly believe that Arlequin will count an indel '-' as an informative site. It should just ignore those sites, because it is a very respectable package that has been around for years and years.</p>
<p>Okay, I am now signing off for Christmas... with occasional exceptions.</p>

CC BY-SA 4.0
2022-12-20T16:29:03.920Z
2022-12-20T16:29:03.920Z
M__
<p>I have mtDNA sequences of the <em>cox1</em> gene from various populations of one species. I want to perform an AMOVA allele frequency analysis via Arlequin. This allele frequency package however requires my sequences to have the same length.</p>
<ul>
<li>Is it possible to use my sequences?</li>
<li>Do I have to align them first maybe?</li>
</ul>
<p>The length difference is most likely due to indels/deletions. Do I have to use specific parts of the sequence for each gene of the subunit?</p>

CC BY-SA 4.0
2022-12-19T10:14:22.227Z
Please [edit] your question and give us some more details. Are you talking about protein sequences? DNA? mRNA? Arlequin can perform various types of analysis, what analysis are you trying to do?
terdon
CC BY-SA 4.0
2022-12-19T10:46:08.213Z
Thanks, that helps. I don't know what you mean by subunit though, do you mean subunit of the gene's protein product, of cytochrome c oxidase? Or are you working on some sort of protein complex formed by the protein products of multiple genes? In any case, yes, you would very likely need to align the sequences first. Why not at least try that and see if it works?
terdon
CC BY-SA 4.0
2022-12-19T12:52:00.597Z
Oh BTW @Nickmofoe the correct name for the gene is cox1. Please avoid using the old nomenclature of COI, this can create confusion.
M__
CC BY-SA 4.0
2022-12-19T19:15:11.497Z
There are various workflows that do this that are easy to google, for example [here](https://stat.ethz.ch/pipermail/r-sig-genetics/2017-January/000135.html). Can you clarify how these don't work for you?
Maximilian Press
<phylogenetics><genetics><allele-frequency><populations>
How can I use sequences with different lengths in arlequin?
# How can I use sequences with different lengths in arlequin?

(Asked by: username_0 on 2022-12-19)

I have mtDNA sequences of the *cox1* gene from various populations of one species. I want to perform an AMOVA allele frequency analysis via Arlequin. This allele frequency package however requires my sequences to have the same length.

* Is it possible to use my sequences?
* Do I have to align them first maybe?

The length difference is most likely due to indels/deletions. Do I have to use specific parts of the sequence for each gene of the subunit?

**Tags:** phylogenetics, genetics, allele-frequency, populations
> **Comments:**
> **username_1**: Please [edit] your question and give us some more details. Are you talking about protein sequences? DNA? mRNA? Arlequin can perform various types of analysis, what analysis are you trying to do?  
> **username_1**: Thanks, that helps. I don't know what you mean by subunit though, do you mean subunit of the gene's protein product, of cytochrome c oxidase? Or are you working on some sort of protein complex formed by the protein products of multiple genes? In any case, yes, you would very likely need to align the sequences first. Why not at least try that and see if it works?  
> **username_2**: Oh BTW @Nickmofoe the correct name for the gene is cox1. Please avoid using the old nomenclature of COI, this can create confusion.  
> **username_3**: There are various workflows that do this that are easy to google, for example [here](https://stat.ethz.ch/pipermail/r-sig-genetics/2017-January/000135.html). Can you clarify how these don't work for you?  

---

## 1 answers

### Answer by username_2 on 2022-12-20. Score: 0 | Accepted answer

All that is needed is to align the sequences within a formal aligner like MAFFT or muscle. Then pass that through the PDGSpider converter. Arlequin should just take it from there.

Manual annotation of the alignment is normal e.g. `seqoTron` and is traditionally regarded as the best approach rather than do this algorithmically. Thus I'm suggesting (strongly suggesting) do this prior the PGDSpider conversion. If an alignment has not been performed Arlequin should throw an error.

If there are trailing gaps (indels or '-') at the beginning and end of the alignment Arlequin will/should ignore these. You can of course simply cut them out with the alignment editor (Mega X will also work for alignment editing). The important point is that Arlequin must correctly identify homologues its an allele frequency package and nothing more. If there are errors in the alignment Arlequin will believe there is population structure when it could be absent.

I can't honestly believe that Arlequin will count an indel '-' as an informative site. It should just ignore those sites, because it is a very respectable package that has been around for years and years.

Okay, I am now signing off for Christmas... with occasional exceptions.


---

