Legofit
infers population history from nucleotide site patterns.
|
Convert axt format to raf format.
The axt format is used to describe the segments of one genome that align to another. For example, there are axt files that describe the alignment of the chimpanzee reference sequence to human one. Given such a file, axt2raf.py will generate a .raf file that maps the bases of the chimpanzee sequence onto human coordinates. Because reference genomes are haploid, the output .raf file will have only 0s and 1s as reference allele frequencies.
usage: axt2raf.py [options] [inputfile] where options may include: -minlen <x> : set min alignment length to <x> base pairs -minqual <x> : set min alignment quality to <x> Program reads from "inputfile" if provided, or otherwise from standard input. Input should be in axt format. Writes to standard output.