application: est2genome [
  documentation: "Align EST sequences to genomic DNA sequence"
  groups: "Alignment:Global"
  gui: "yes"
  batch: "yes"
  cpu: "medium"
    relations: "EDAM_topic:0182 Sequence alignment"
    relations: "EDAM_operation:0292 Sequence alignment construction"
]

section: input [
  information: "Input section"
  type: "page"
]

  seqall: estsequence [
    parameter: "Y"
    type: "dna"
    information: "Spliced EST nucleotide sequence(s)"
    knowntype: "est sequence"
    relations: "EDAM_data:2887 Sequence record (nucleic acid)"
  ]

  sequence: genomesequence [
    parameter: "Y"
    type: "dna"
    information: "Unspliced genomic nucleotide sequence"
    knowntype: "genomic sequence"
    relations: "EDAM_data:2887 Sequence record (nucleic acid)"
  ]

endsection: input

section: additional [
  information: "Additional section"
  type: "page"
]

  integer: match [
    additional: "Y"
    default: "1"
    information: "Score for matching two bases"
    relations: "EDAM_data:1401 Match reward score"
  ]

  integer: mismatch [
    additional: "Y"
    default: "1"
    information: "Cost for mismatching two bases"
    relations: "EDAM_data:1402 Mismatch penalty score"
  ]

  integer: gappenalty [
    additional: "Y"
    default: "2"
    information: "Gap penalty"
    help: "Cost for deleting a single base in either sequence,
           excluding introns"
    relations: "EDAM_data:2137 Gap penalty"
  ]

  integer: intronpenalty [
    additional: "Y"
    default: "40"
    information: "Intron penalty"
    help: "Cost for an intron, independent of length."
    relations: "EDAM_data:1772 Score or penalty"
  ]

  integer: splicepenalty [
    additional: "Y"
    default: "20"
    information: "Splice site penalty"
    help: "Cost for an intron, independent of length and
           starting/ending on donor-acceptor sites"
    relations: "EDAM_data:1772 Score or penalty"
  ]

  integer: minscore [
    additional: "Y"
    default: "30"
    information: "Minimum accepted score"
    help: "Exclude alignments with scores below this threshold score."
    relations: "EDAM_data:2146 Threshold"
  ]

endsection: additional

section: advanced [
  information: "Advanced section"
  type: "page"
]

  boolean: reverse [
    information: "Reverse orientation"
    help: "Reverse the orientation of the EST sequence"
    relations: "EDAM_data:2527 Parameter"
  ]

  boolean: usesplice [
    default: "Y"
    information: "Use donor and acceptor splice sites"
    help: "Use donor and acceptor splice sites. If you want to ignore
           donor-acceptor sites then set this to be false."
    relations: "EDAM_data:2527 Parameter"
  ]

  list: mode [
    default: "both"
    information: "Comparison mode"
    help: "This determines the comparison mode. The default value is
           'both', in which case both strands of the est are compared
           assuming a forward gene direction (ie GT/AG splice sites), and the
           best comparison redone assuming a reversed (CT/AC) gene splicing
           direction. The other allowed modes are 'forward', when just the
           forward strand is searched, and 'reverse', ditto for the reverse
           strand."
    values: "both:Both strands; forward:Forward strand only;
             reverse:Reverse strand only"
    relations: "EDAM_data:2527 Parameter"
  ]

  boolean: best [
    default: "Y"
    information: "Print out only best alignment"
    help: "You can print out all comparisons instead of just the best
           one by setting this to be false."
    relations: "EDAM_data:2527 Parameter"
  ]

  float: space [
    default: "10.0"
    information: "Space threshold (in megabytes)"
    help: "For linear-space recursion. If product of sequence lengths
           divided by 4 exceeds this then a divide-and-conquer strategy is
           used to control the memory requirements. In this way very long
           sequences can be aligned. \
           If you have a machine with plenty of memory you can raise this
           parameter (but do not exceed the machine's physical RAM)"
    relations: "EDAM_data:2527 Parameter"
  ]

  integer: shuffle [
    information: "Shuffle"
    relations: "EDAM_data:2527 Parameter"
  ]

  integer: seed [
    default: "20825"
    information: "Random number seed"
    relations: "EDAM_data:2527 Parameter"
  ]

endsection: advanced

section: output [
  information: "Output section"
  type: "page"
]

  outfile: outfile [
    parameter: "y"
    knowntype: "est2genome output"
    relations: "EDAM_data:1313 Nucleic acid features (coding sequence)"
  ]

  boolean: align [
    information: "Show the alignment"
    help: "Show the alignment. The alignment includes the first and
           last 5 bases of each intron, together with the intron width. The
           direction of splicing is indicated by angle brackets (forward or
           reverse) or ???? (unknown)."
    relations: "EDAM_data:2527 Parameter"
  ]

  integer: width [
    default: "50"
    information: "Alignment width"
    relations: "EDAM_data:2527 Parameter"
  ]

endsection: output
