code: def parse_insertion(insertion, gff):
      
      offset = insertion['offset']
      for ins in parse_fasta(insertion['insertion sequence'].split('|')):
          strand = insertion['strand']
          ID = ins[0].split('>')[1].split()[0]
          Start, End = [int(i) for i in ins[0].split('gene-pos=', 1)[1].split()[0].split('-')]
          Start, End = abs(Start + offset), abs(End + offset)
          if strand == '-':
              Start, End = End, Start
          gff['#seqname'].append(insertion['ID'])
          gff['source'].append(insertion['source'])
          gff['feature'].append('IVS')
          gff['start'].append(Start)
          gff['end'].append(End)
          gff['score'].append('.')
          gff['strand'].append(strand) # same as rRNA
          gff['frame'].append('.')
          gff['attribute'].append('ID=%s' % (ID))
      return gff
docs: parse insertion to gff format