code: def simulate(args):
      
      if args.fasta:
          name = None
          seq = ""
          reads = dict()
          with open(args.fasta) as in_handle:
              for line in in_handle:
                  if line.startswith(">"):
                      if name:
                          reads.update(_generate_reads(seq, name))
                      seq = ""
                      name = line[1:-1]
                  else:
                      seq += line.strip()
  
          reads.update(_generate_reads(seq, name))
      _write_reads(reads, args.out)
docs: Main function that manage simulatin of small RNAs