Provided by: genometools_1.5.1-2ubuntu1_amd64 bug


       gt-extractfeat - Extract features given in GFF3 file from sequence file.


       gt extractfeat [option ...] [GFF3_file]


           set type of features to extract (default: undefined)

           join feature sequences in the same subgraph into a single one (default: no)

           translate the features (of a DNA sequence) into protein (default: no)

           add sequence ID of extracted features to FASTA descriptions (default: no)

           add target ID(s) of extracted features to FASTA descriptions (default: no)

           set the sequence file from which to extract the features (default: undefined)

           set the sequence files from which to extract the features use -- to terminate the list
           of sequence files

           match the sequence descriptions from the input files for the desired sequence IDs (in
           GFF3) (default: no)

           use sequence descriptions to map the sequence IDs (in GFF3) to actual sequence
           entries. If a description contains a sequence range (e.g., III:1000001..2000000), the
           first part is used as sequence ID (III) and the first range position as offset
           (1000001) (default: no)

           set file containing sequence-region to sequence file mapping (default: undefined)

           be verbose (default: no)

           set output width for FASTA sequence printing (0 disables formatting) (default: 0)

           redirect output to specified file (default: undefined)

           write gzip compressed output file (default: no)

           write bzip2 compressed output file (default: no)

           force writing to output file (default: no)

           display help and exit

           display version information and exit


       File format for option -regionmapping:

       The file supplied to option -regionmapping defines a “mapping”. A mapping maps the
       sequence-regions given in the GFF3_file to a sequence file containing the corresponding
       sequence. Mappings can be defined in one of the following two forms:

       mapping = { chr1 = "hs_ref_chr1.fa.gz", chr2 = "hs_ref_chr2.fa.gz" }


       function mapping(sequence_region) return "hs_ref_"..sequence_region..".fa.gz" end

       The first form defines a Lua ( table named “mapping” which maps each
       sequence region to the corresponding sequence file. The second one defines a Lua function
       “mapping”, which has to return the sequence file name when it is called with the
       sequence_region as argument.


       Report bugs to <>.

                                            10/07/2012                          GT-EXTRACTFEAT(1)