|
|
cutseq |
Please help by correcting and extending the Wiki pages.
This simple editing program allows you to cut out a region from your sequence by specifying the begin and end positions of the region to remove. It removes the sequence from the specified start to the end positions (inclusive) and writes the remaining sequence to the output file.
To remove bases 10 to 12 from a database entry and write to the new sequence file 'gatta2.seq':
% cutseq tembl:x13776 gatta2.seq -from=10 -to=12 Remove a section from a sequence |
Go to the input files for this example
Go to the output files for this example
Example 2
To remove the first 20 bases from 'tembl:x13776' and write it to 'jsh.seq':
% cutseq tembl:x13776 -from=1 -to=20 -out=jsh.seq Remove a section from a sequence |
Go to the output files for this example
Example 3
If the default start and end positions are accepted, then all of the sequence is removed!
% cutseq tembl:x13776 starta.seq -sbeg=-1000 -send=1290 Remove a section from a sequence Start of region to delete [1168]: End of region to delete [1290]: |
Go to the output files for this example
Remove a section from a sequence
Version: EMBOSS:6.5.6.0
Standard (Mandatory) qualifiers:
[-sequence] sequence (Gapped) sequence filename and optional
format, or reference (input USA)
-from integer [Start of sequence (0)] This is the start
position (inclusive) of the section of the
sequence that you wish to remove. (Any
integer value)
-to integer [End of sequence (0)] This is the end
position (inclusive) of the section of the
sequence that you wish to remove. (Any
integer value)
[-outseq] seqout [
|
| Qualifier | Type | Description | Allowed values | Default |
|---|---|---|---|---|
| Standard (Mandatory) qualifiers | ||||
| [-sequence] (Parameter 1) |
sequence | (Gapped) sequence filename and optional format, or reference (input USA) | Readable sequence | Required |
| -from | integer | This is the start position (inclusive) of the section of the sequence that you wish to remove. | Any integer value | Start of sequence (0) |
| -to | integer | This is the end position (inclusive) of the section of the sequence that you wish to remove. | Any integer value | End of sequence (0) |
| [-outseq] (Parameter 2) |
seqout | Sequence filename and optional format (output USA) | Writeable sequence | <*>.format |
| Additional (Optional) qualifiers | ||||
| (none) | ||||
| Advanced (Unprompted) qualifiers | ||||
| (none) | ||||
| Associated qualifiers | ||||
| "-sequence" associated sequence qualifiers | ||||
| -sbegin1 -sbegin_sequence |
integer | Start of the sequence to be used | Any integer value | 0 |
| -send1 -send_sequence |
integer | End of the sequence to be used | Any integer value | 0 |
| -sreverse1 -sreverse_sequence |
boolean | Reverse (if DNA) | Boolean value Yes/No | N |
| -sask1 -sask_sequence |
boolean | Ask for begin/end/reverse | Boolean value Yes/No | N |
| -snucleotide1 -snucleotide_sequence |
boolean | Sequence is nucleotide | Boolean value Yes/No | N |
| -sprotein1 -sprotein_sequence |
boolean | Sequence is protein | Boolean value Yes/No | N |
| -slower1 -slower_sequence |
boolean | Make lower case | Boolean value Yes/No | N |
| -supper1 -supper_sequence |
boolean | Make upper case | Boolean value Yes/No | N |
| -scircular1 -scircular_sequence |
boolean | Sequence is circular | Boolean value Yes/No | N |
| -sformat1 -sformat_sequence |
string | Input sequence format | Any string | |
| -iquery1 -iquery_sequence |
string | Input query fields or ID list | Any string | |
| -ioffset1 -ioffset_sequence |
integer | Input start position offset | Any integer value | 0 |
| -sdbname1 -sdbname_sequence |
string | Database name | Any string | |
| -sid1 -sid_sequence |
string | Entryname | Any string | |
| -ufo1 -ufo_sequence |
string | UFO features | Any string | |
| -fformat1 -fformat_sequence |
string | Features format | Any string | |
| -fopenfile1 -fopenfile_sequence |
string | Features file name | Any string | |
| "-outseq" associated seqout qualifiers | ||||
| -osformat2 -osformat_outseq |
string | Output seq format | Any string | |
| -osextension2 -osextension_outseq |
string | File name extension | Any string | |
| -osname2 -osname_outseq |
string | Base file name | Any string | |
| -osdirectory2 -osdirectory_outseq |
string | Output directory | Any string | |
| -osdbname2 -osdbname_outseq |
string | Database name to add | Any string | |
| -ossingle2 -ossingle_outseq |
boolean | Separate file for each entry | Boolean value Yes/No | N |
| -oufo2 -oufo_outseq |
string | UFO features | Any string | |
| -offormat2 -offormat_outseq |
string | Features format | Any string | |
| -ofname2 -ofname_outseq |
string | Features file name | Any string | |
| -ofdirectory2 -ofdirectory_outseq |
string | Output directory | Any string | |
| General qualifiers | ||||
| -auto | boolean | Turn off prompts | Boolean value Yes/No | N |
| -stdout | boolean | Write first file to standard output | Boolean value Yes/No | N |
| -filter | boolean | Read first file from standard input, write first file to standard output | Boolean value Yes/No | N |
| -options | boolean | Prompt for standard and additional values | Boolean value Yes/No | N |
| -debug | boolean | Write debug output to program.dbg | Boolean value Yes/No | N |
| -verbose | boolean | Report some/full command line options | Boolean value Yes/No | Y |
| -help | boolean | Report command line options and exit. More information on associated and general qualifiers can be found with -help -verbose | Boolean value Yes/No | N |
| -warning | boolean | Report warnings | Boolean value Yes/No | Y |
| -error | boolean | Report errors | Boolean value Yes/No | Y |
| -fatal | boolean | Report fatal errors | Boolean value Yes/No | Y |
| -die | boolean | Report dying program messages | Boolean value Yes/No | Y |
| -version | boolean | Report version number and exit | Boolean value Yes/No | N |
The input is a standard EMBOSS sequence query (also known as a 'USA').
Major sequence database sources defined as standard in EMBOSS installations include srs:embl, srs:uniprot and ensembl
Data can also be read from sequence output in any supported format written by an EMBOSS or third-party application.
The input format can be specified by using the command-line qualifier -sformat xxx, where 'xxx' is replaced by the name of the required format. The available format names are: gff (gff3), gff2, embl (em), genbank (gb, refseq), ddbj, refseqp, pir (nbrf), swissprot (swiss, sw), dasgff and debug.
See: http://emboss.sf.net/docs/themes/SequenceFormats.html for further information on sequence formats.
ID X13776; SV 1; linear; genomic DNA; STD; PRO; 2167 BP.
XX
AC X13776; M43175;
XX
DT 19-APR-1989 (Rel. 19, Created)
DT 14-NOV-2006 (Rel. 89, Last updated, Version 24)
XX
DE Pseudomonas aeruginosa amiC and amiR gene for aliphatic amidase regulation
XX
KW aliphatic amidase regulator; amiC gene; amiR gene.
XX
OS Pseudomonas aeruginosa
OC Bacteria; Proteobacteria; Gammaproteobacteria; Pseudomonadales;
OC Pseudomonadaceae; Pseudomonas.
XX
RN [1]
RP 1167-2167
RA Rice P.M.;
RT ;
RL Submitted (16-DEC-1988) to the INSDC.
RL Rice P.M., EMBL, Postfach 10-2209, Meyerhofstrasse 1, 6900 Heidelberg, FRG.
XX
RN [2]
RP 1167-2167
RX DOI; 10.1016/0014-5793(89)80249-2.
RX PUBMED; 2495988.
RA Lowe N., Rice P.M., Drew R.E.;
RT "Nucleotide sequence of the aliphatic amidase regulator gene (amiR) of
RT Pseudomonas aeruginosa";
RL FEBS Lett. 246(1-2):39-43(1989).
XX
RN [3]
RP 1-1292
RX PUBMED; 1907262.
RA Wilson S., Drew R.;
RT "Cloning and DNA sequence of amiC, a new gene regulating expression of the
RT Pseudomonas aeruginosa aliphatic amidase, and purification of the amiC
RT product";
RL J. Bacteriol. 173(16):4914-4921(1991).
XX
RN [4]
RP 1-2167
RA Rice P.M.;
RT ;
RL Submitted (04-SEP-1991) to the INSDC.
RL Rice P.M., EMBL, Postfach 10-2209, Meyerhofstrasse 1, 6900 Heidelberg, FRG.
XX
DR GOA; Q51417.
DR InterPro; IPR003211; AmiSUreI_transpt.
DR UniProtKB/Swiss-Prot; Q51417; AMIS_PSEAE.
[Part of this file has been deleted for brevity]
FT /note="ClaI fragment deleted in pSW36, constitutive
FT phenotype"
FT misc_feature 1
FT /note="last base of an XhoI site"
FT misc_feature 648..653
FT /note="end of 658bp XhoI fragment, deletion in pSW3 causes
FT constitutive expression of amiE"
FT misc_difference 1281
FT /replace="g"
FT /note="conflict"
FT /citation=[3]
XX
SQ Sequence 2167 BP; 363 A; 712 C; 730 G; 362 T; 0 other;
ggtaccgctg gccgagcatc tgctcgatca ccaccagccg ggcgacggga actgcacgat 60
ctacctggcg agcctggagc acgagcgggt tcgcttcgta cggcgctgag cgacagtcac 120
aggagaggaa acggatggga tcgcaccagg agcggccgct gatcggcctg ctgttctccg 180
aaaccggcgt caccgccgat atcgagcgct cgcacgcgta tggcgcattg ctcgcggtcg 240
agcaactgaa ccgcgagggc ggcgtcggcg gtcgcccgat cgaaacgctg tcccaggacc 300
ccggcggcga cccggaccgc tatcggctgt gcgccgagga cttcattcgc aaccgggggg 360
tacggttcct cgtgggctgc tacatgtcgc acacgcgcaa ggcggtgatg ccggtggtcg 420
agcgcgccga cgcgctgctc tgctacccga ccccctacga gggcttcgag tattcgccga 480
acatcgtcta cggcggtccg gcgccgaacc agaacagtgc gccgctggcg gcgtacctga 540
ttcgccacta cggcgagcgg gtggtgttca tcggctcgga ctacatctat ccgcgggaaa 600
gcaaccatgt gatgcgccac ctgtatcgcc agcacggcgg cacggtgctc gaggaaatct 660
acattccgct gtatccctcc gacgacgact tgcagcgcgc cgtcgagcgc atctaccagg 720
cgcgcgccga cgtggtcttc tccaccgtgg tgggcaccgg caccgccgag ctgtatcgcg 780
ccatcgcccg tcgctacggc gacggcaggc ggccgccgat cgccagcctg accaccagcg 840
aggcggaggt ggcgaagatg gagagtgacg tggcagaggg gcaggtggtg gtcgcgcctt 900
acttctccag catcgatacg cccgccagcc gggccttcgt ccaggcctgc catggtttct 960
tcccggagaa cgcgaccatc accgcctggg ccgaggcggc ctactggcag accttgttgc 1020
tcggccgcgc cgcgcaggcc gcaggcaact ggcgggtgga agacgtgcag cggcacctgt 1080
acgacatcga catcgacgcg ccacaggggc cggtccgggt ggagcgccag aacaaccaca 1140
gccgcctgtc ttcgcgcatc gcggaaatcg atgcgcgcgg cgtgttccag gtccgctggc 1200
agtcgcccga accgattcgc cccgaccctt atgtcgtcgt gcataacctc gacgactggt 1260
ccgccagcat gggcggggga ccgctcccat gagcgccaac tcgctgctcg gcagcctgcg 1320
cgagttgcag gtgctggtcc tcaacccgcc gggggaggtc agcgacgccc tggtcttgca 1380
gctgatccgc atcggttgtt cggtgcgcca gtgctggccg ccgccggaag ccttcgacgt 1440
gccggtggac gtggtcttca ccagcatttt ccagaatggc caccacgacg agatcgctgc 1500
gctgctcgcc gccgggactc cgcgcactac cctggtggcg ctggtggagt acgaaagccc 1560
cgcggtgctc tcgcagatca tcgagctgga gtgccacggc gtgatcaccc agccgctcga 1620
tgcccaccgg gtgctgcctg tgctggtatc ggcgcggcgc atcagcgagg aaatggcgaa 1680
gctgaagcag aagaccgagc agctccagga ccgcatcgcc ggccaggccc ggatcaacca 1740
ggccaaggtg ttgctgatgc agcgccatgg ctgggacgag cgcgaggcgc accagcacct 1800
gtcgcgggaa gcgatgaagc ggcgcgagcc gatcctgaag atcgctcagg agttgctggg 1860
aaacgagccg tccgcctgag cgatccgggc cgaccagaac aataacaaga ggggtatcgt 1920
catcatgctg ggactggttc tgctgtacgt tggcgcggtg ctgtttctca atgccgtctg 1980
gttgctgggc aagatcagcg gtcgggaggt ggcggtgatc aacttcctgg tcggcgtgct 2040
gagcgcctgc gtcgcgttct acctgatctt ttccgcagca gccgggcagg gctcgctgaa 2100
ggccggagcg ctgaccctgc tattcgcttt tacctatctg tgggtggccg ccaaccagtt 2160
cctcgag 2167
//
|
The output is a standard EMBOSS sequence file.
The results can be output in one of several styles by using the command-line qualifier -osformat xxx, where 'xxx' is replaced by the name of the required format. The available format names are: embl, genbank, gff, pir, swiss, dasgff, debug, listfile, dbmotif, diffseq, excel, feattable, motif, nametable, regions, seqtable, simple, srs, table, tagseq.
See: http://emboss.sf.net/docs/themes/SequenceFormats.html for further information on sequence formats.
>X13776 X13776.1 Pseudomonas aeruginosa amiC and amiR gene for aliphatic amidase regulation ggtaccgctcgagcatctgctcgatcaccaccagccgggcgacgggaactgcacgatcta cctggcgagcctggagcacgagcgggttcgcttcgtacggcgctgagcgacagtcacagg agaggaaacggatgggatcgcaccaggagcggccgctgatcggcctgctgttctccgaaa ccggcgtcaccgccgatatcgagcgctcgcacgcgtatggcgcattgctcgcggtcgagc aactgaaccgcgagggcggcgtcggcggtcgcccgatcgaaacgctgtcccaggaccccg gcggcgacccggaccgctatcggctgtgcgccgaggacttcattcgcaaccggggggtac ggttcctcgtgggctgctacatgtcgcacacgcgcaaggcggtgatgccggtggtcgagc gcgccgacgcgctgctctgctacccgaccccctacgagggcttcgagtattcgccgaaca tcgtctacggcggtccggcgccgaaccagaacagtgcgccgctggcggcgtacctgattc gccactacggcgagcgggtggtgttcatcggctcggactacatctatccgcgggaaagca accatgtgatgcgccacctgtatcgccagcacggcggcacggtgctcgaggaaatctaca ttccgctgtatccctccgacgacgacttgcagcgcgccgtcgagcgcatctaccaggcgc gcgccgacgtggtcttctccaccgtggtgggcaccggcaccgccgagctgtatcgcgcca tcgcccgtcgctacggcgacggcaggcggccgccgatcgccagcctgaccaccagcgagg cggaggtggcgaagatggagagtgacgtggcagaggggcaggtggtggtcgcgccttact tctccagcatcgatacgcccgccagccgggccttcgtccaggcctgccatggtttcttcc cggagaacgcgaccatcaccgcctgggccgaggcggcctactggcagaccttgttgctcg gccgcgccgcgcaggccgcaggcaactggcgggtggaagacgtgcagcggcacctgtacg acatcgacatcgacgcgccacaggggccggtccgggtggagcgccagaacaaccacagcc gcctgtcttcgcgcatcgcggaaatcgatgcgcgcggcgtgttccaggtccgctggcagt cgcccgaaccgattcgccccgacccttatgtcgtcgtgcataacctcgacgactggtccg ccagcatgggcgggggaccgctcccatgagcgccaactcgctgctcggcagcctgcgcga gttgcaggtgctggtcctcaacccgccgggggaggtcagcgacgccctggtcttgcagct gatccgcatcggttgttcggtgcgccagtgctggccgccgccggaagccttcgacgtgcc ggtggacgtggtcttcaccagcattttccagaatggccaccacgacgagatcgctgcgct gctcgccgccgggactccgcgcactaccctggtggcgctggtggagtacgaaagccccgc ggtgctctcgcagatcatcgagctggagtgccacggcgtgatcacccagccgctcgatgc ccaccgggtgctgcctgtgctggtatcggcgcggcgcatcagcgaggaaatggcgaagct gaagcagaagaccgagcagctccaggaccgcatcgccggccaggcccggatcaaccaggc caaggtgttgctgatgcagcgccatggctgggacgagcgcgaggcgcaccagcacctgtc gcgggaagcgatgaagcggcgcgagccgatcctgaagatcgctcaggagttgctgggaaa cgagccgtccgcctgagcgatccgggccgaccagaacaataacaagaggggtatcgtcat catgctgggactggttctgctgtacgttggcgcggtgctgtttctcaatgccgtctggtt gctgggcaagatcagcggtcgggaggtggcggtgatcaacttcctggtcggcgtgctgag cgcctgcgtcgcgttctacctgatcttttccgcagcagccgggcagggctcgctgaaggc cggagcgctgaccctgctattcgcttttacctatctgtgggtggccgccaaccagttcct cgag |
>X13776 X13776.1 Pseudomonas aeruginosa amiC and amiR gene for aliphatic amidase regulation tgctcgatcaccaccagccgggcgacgggaactgcacgatctacctggcgagcctggagc acgagcgggttcgcttcgtacggcgctgagcgacagtcacaggagaggaaacggatggga tcgcaccaggagcggccgctgatcggcctgctgttctccgaaaccggcgtcaccgccgat atcgagcgctcgcacgcgtatggcgcattgctcgcggtcgagcaactgaaccgcgagggc ggcgtcggcggtcgcccgatcgaaacgctgtcccaggaccccggcggcgacccggaccgc tatcggctgtgcgccgaggacttcattcgcaaccggggggtacggttcctcgtgggctgc tacatgtcgcacacgcgcaaggcggtgatgccggtggtcgagcgcgccgacgcgctgctc tgctacccgaccccctacgagggcttcgagtattcgccgaacatcgtctacggcggtccg gcgccgaaccagaacagtgcgccgctggcggcgtacctgattcgccactacggcgagcgg gtggtgttcatcggctcggactacatctatccgcgggaaagcaaccatgtgatgcgccac ctgtatcgccagcacggcggcacggtgctcgaggaaatctacattccgctgtatccctcc gacgacgacttgcagcgcgccgtcgagcgcatctaccaggcgcgcgccgacgtggtcttc tccaccgtggtgggcaccggcaccgccgagctgtatcgcgccatcgcccgtcgctacggc gacggcaggcggccgccgatcgccagcctgaccaccagcgaggcggaggtggcgaagatg gagagtgacgtggcagaggggcaggtggtggtcgcgccttacttctccagcatcgatacg cccgccagccgggccttcgtccaggcctgccatggtttcttcccggagaacgcgaccatc accgcctgggccgaggcggcctactggcagaccttgttgctcggccgcgccgcgcaggcc gcaggcaactggcgggtggaagacgtgcagcggcacctgtacgacatcgacatcgacgcg ccacaggggccggtccgggtggagcgccagaacaaccacagccgcctgtcttcgcgcatc gcggaaatcgatgcgcgcggcgtgttccaggtccgctggcagtcgcccgaaccgattcgc cccgacccttatgtcgtcgtgcataacctcgacgactggtccgccagcatgggcggggga ccgctcccatgagcgccaactcgctgctcggcagcctgcgcgagttgcaggtgctggtcc tcaacccgccgggggaggtcagcgacgccctggtcttgcagctgatccgcatcggttgtt cggtgcgccagtgctggccgccgccggaagccttcgacgtgccggtggacgtggtcttca ccagcattttccagaatggccaccacgacgagatcgctgcgctgctcgccgccgggactc cgcgcactaccctggtggcgctggtggagtacgaaagccccgcggtgctctcgcagatca tcgagctggagtgccacggcgtgatcacccagccgctcgatgcccaccgggtgctgcctg tgctggtatcggcgcggcgcatcagcgaggaaatggcgaagctgaagcagaagaccgagc agctccaggaccgcatcgccggccaggcccggatcaaccaggccaaggtgttgctgatgc agcgccatggctgggacgagcgcgaggcgcaccagcacctgtcgcgggaagcgatgaagc ggcgcgagccgatcctgaagatcgctcaggagttgctgggaaacgagccgtccgcctgag cgatccgggccgaccagaacaataacaagaggggtatcgtcatcatgctgggactggttc tgctgtacgttggcgcggtgctgtttctcaatgccgtctggttgctgggcaagatcagcg gtcgggaggtggcggtgatcaacttcctggtcggcgtgctgagcgcctgcgtcgcgttct acctgatcttttccgcagcagccgggcagggctcgctgaaggccggagcgctgaccctgc tattcgcttttacctatctgtgggtggccgccaaccagttcctcgag |
>X13776 X13776.1 Pseudomonas aeruginosa amiC and amiR gene for aliphatic amidase regulation |
Qualifiers that are in-built for sequence datatypes allow the input and output files to be specified in more detail, for example, the format for the output files.
| Program name | Description |
|---|---|
| aligncopy | Read and write alignments |
| aligncopypair | Read and write pairs from alignments |
| biosed | Replace or delete sequence sections |
| codcopy | Copy and reformat a codon usage table |
| degapseq | Remove non-alphabetic (e.g. gap) characters from sequences |
| descseq | Alter the name or description of a sequence |
| entret | Retrieve sequence entries from flatfile databases and files |
| extractalign | Extract regions from a sequence alignment |
| extractfeat | Extract features from sequence(s) |
| extractseq | Extract regions from a sequence |
| featcopy | Read and write a feature table |
| featmerge | Merge two overlapping feature tables |
| featreport | Read and write a feature table |
| feattext | Return a feature table original text |
| listor | Write a list file of the logical OR of two sets of sequences |
| makenucseq | Create random nucleotide sequences |
| makeprotseq | Create random protein sequences |
| maskambignuc | Mask all ambiguity characters in nucleotide sequences with N |
| maskambigprot | Mask all ambiguity characters in protein sequences with X |
| maskfeat | Write a sequence with masked features |
| maskseq | Write a sequence with masked regions |
| newseq | Create a sequence file from a typed-in sequence |
| nohtml | Remove mark-up (e.g. HTML tags) from an ASCII text file |
| noreturn | Remove carriage return from ASCII files |
| nospace | Remove whitespace from an ASCII text file |
| notab | Replace tabs with spaces in an ASCII text file |
| notseq | Write to file a subset of an input stream of sequences |
| nthseq | Write to file a single sequence from an input stream of sequences |
| nthseqset | Read and write (return) one set of sequences from many |
| pasteseq | Insert one sequence into another |
| revseq | Reverse and complement a nucleotide sequence |
| seqcount | Read and count sequences |
| seqret | Read and write (return) sequences |
| seqretsetall | Read and write (return) many sets of sequences |
| seqretsplit | Read sequences and write them to individual files |
| sizeseq | Sort sequences by size |
| skipredundant | Remove redundant sequences from an input set |
| skipseq | Read and write (return) sequences, skipping first few |
| splitsource | Split sequence(s) into original source sequences |
| splitter | Split sequence(s) into smaller sequences |
| trimest | Remove poly-A tails from nucleotide sequences |
| trimseq | Remove unwanted characters from start and end of sequence(s) |
| trimspace | Remove extra whitespace from an ASCII text file |
| union | Concatenate multiple sequences into a single sequence |
| vectorstrip | Remove vectors from the ends of nucleotide sequence(s) |
| yank | Add a sequence reference (a full USA) to a list file |
Please report all bugs to the EMBOSS bug team (emboss-bug © emboss.open-bio.org) not to the original author.