Re: Row order from multiple sequence alignment
Spencer Bliven <[email protected]> Tue, 1 Mar 2016 11:16:12 +0100
| Newsgroups | gmane.comp.java.bio.general |
|---|---|
| Message-ID | <CA+P6ar=5r9ohS=RcFkHPKJgnQwJy2-JbPrfHKEP9w44qKyoC1A@mail.gmail.com> |
The original sequences are also in a permuted order. I took your suggestion to store the order in the accessions. That worked fine. Thanks, Spencer On Tue, Mar 1, 2016 at 7:43 AM, Andreas Prlic <[email protected]> wrote: > Hi Spencer, > > You might have to compare the accession codes of the sequences... > > // 1 is the first sequence in the alignment, not 0... > AlignedSequence<ProteinSequence, AminoAcidCompound> as = > profile.getAlignedSequence(1); > System.out.println(as.getAccession()); > > Alternatively, what is the order of sequences that comes out of profile.getOriginalSequences() > ? > > > A > > > > On Mon, Feb 29, 2016 at 9:15 AM, Spencer Bliven <[email protected]> > wrote: > >> I'm creating a multiple sequence alignment using >> Alignments.getMultipleSequenceAlignment, as described in the cookbook >> <http://biojava.org/wiki/BioJava:CookBook3:MSA>. The problem is that the >> returned profile has rows in a different order than the input array of >> sequenced. How can I map from an index in the inputs to the rows of the >> profile? >> >> Thanks, >> Spencer >> >> _______________________________________________ >> Biojava-l mailing list - [email protected] >> http://mailman.open-bio.org/mailman/listinfo/biojava-l >> > > > > -- > ----------------------------------------------------------------------- > Dr. Andreas Prlic > RCSB PDB Protein Data Bank > Technical & Scientific Team Lead > University of California, San Diego > > Editor Software Section > PLOS Computational Biology > > BioJava Project Lead > ----------------------------------------------------------------------- > _______________________________________________ Biojava-l mailing list - [email protected] http://mailman.open-bio.org/mailman/listinfo/biojava-l