[Biopython-dev] Bio.AlignIO, Bio.Nexus, MrBayes, polymorphic sites, maximum line length
    Nick Loman 
    n.j.loman at bham.ac.uk
       
    Thu Dec  2 17:11:13 UTC 2010
    
    
  
Peter wrote:
> Did you ask Bio.Nexus to write interleaved output?
> I've got MrBayes 3.1.2, and this seems to fix your example:
>
> diff --git a/Bio/AlignIO/NexusIO.py b/Bio/AlignIO/NexusIO.py
> index 72550b1..c3b1649 100644
> --- a/Bio/AlignIO/NexusIO.py
> +++ b/Bio/AlignIO/NexusIO.py
> @@ -107,7 +107,7 @@ class NexusWriter(AlignmentWriter):
>          n.alphabet = alignment._alphabet
>          for record in alignment:
>              n.add_sequence(record.id, record.seq.tostring())
> -        n.write_nexus_data(self.handle)
> +        n.write_nexus_data(self.handle, interleave=True)
>
>      def _classify_alphabet_for_nexus(self, alphabet):
>          """Returns 'protein', 'dna', 'rna' based on the alphabet (PRIVATE).
>
>
> Does that work for you?
>   
Hi Peter,
Yes, that does the trick! I wonder if perhaps that's what Cymon meant to 
say rather than interleave_by_partition (hence the boolean problem) ?
Excellent.
Cheers
Nick
    
    
More information about the Biopython-dev
mailing list