Skip to content

Latest commit

 

History

History
56 lines (46 loc) · 1.16 KB

reformat.md

File metadata and controls

56 lines (46 loc) · 1.16 KB

Goalign: toolkit and api for alignment manipulation

API

reformat

Reformat a fasta input alignment into different formats

package main

import (
	"bufio"
	"fmt"
	"io"

	"github.com/evolbioinfo/goalign/align"
	"github.com/evolbioinfo/goalign/io/fasta"
	"github.com/evolbioinfo/goalign/io/nexus"
	"github.com/evolbioinfo/goalign/io/paml"
	"github.com/evolbioinfo/goalign/io/phylip"
	"github.com/evolbioinfo/goalign/io/clustal"
	"github.com/evolbioinfo/goalign/io/utils"
)

func main() {
	var fi io.Closer
	var r *bufio.Reader
	var err error
	var al align.Alignment

	/* Get reader (plain text or gzip) */
	fi, r, err = utils.GetReader("align.fa")
	if err != nil {
		panic(err)
	}

	/* Parse Fasta */
	al, err = fasta.NewParser(r).Parse()
	if err != nil {
		panic(err)
	}
	fi.Close()

	/* Printing FASTA alignment */
	fmt.Println(fasta.WriteAlignment(al))
	/* Printing PHYLIP alignment */
	fmt.Println(phylip.WriteAlignment(al, false))
	/* Printing NEXUS alignment */
	fmt.Println(nexus.WriteAlignment(al))
	/* Printing PAML format alignment */
	fmt.Println(paml.WriteAlignment(al))
    /* Printing Clustal format alignment */
	fmt.Println(clustal.WriteAlignment(al))
}