Updated with code to parse sequence name from fasta header, and to continue work on a file rather than restart using an output file parameter
TreeSequence global alignment