Can LastZ aligner work with IUPAC ambiguity in Fasta files?
0
0
Entering edit mode
5.6 years ago

Hi,

I am trying to compare different tools to align contigs to our reference assembly. The reference contains IUPAC ambiguity codes. LastZ doesn't seem to handle them with the default settings:

lastz  ../../lsalmonis/LSalAtl2s.fasta[multiple] GCA_001005205.1_lsal_atl_canada_female_v1_genomic.fna  --notransition --step=20 --nogapped --format=sam 
FAILURE: bad fasta character in LSalAtl2s.fasta, >LSalAtl2s1: W

LastZ has the 'quantum DNA' feature, that lets one align strings over arbitrary alphabets, but that requires to define a new substitution matrix. I am uncertain if this is the right way to do this.

lastz aligment iupac ambiguity • 1.7k views
ADD COMMENT
2
Entering edit mode

I generally use the option --ambiguous=iupac in these cases

ADD REPLY
0
Entering edit mode

@Michael: Email Bob Harris (email can be found on this page), the author of lastz. He is good about responding to lastz queries.

ADD REPLY

Login before adding your answer.

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