Be careful with coordinates, if your contig have 100 nucleotides and you specify the coord 50,101. the script only write to 100. Make sure the coordinates begin from lower value to upper e.g (1 10, 400 425, etc). IUPAC ambiguity codes of the two possible nucleotides are converted as ...