source: branches/stable/lib/sellists/mod_fields_gc.sellst

Last change on this file was 5390, checked in by westram, 16 years ago
  • TAB-Ex
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.6 KB
Line 
1Count nucleotides,                              nuc#|sequence|count("agctuAGCTU")
2Count nucleotides (incl. IUPAC),                nuc#|sequence|count("agctumrwsykvhdbnAGCTUMRWSYKVHDBN")
3Count nucleotides gene,                         nuc_term#|sequence|filter(SAI=termini,include=x)|count("ACGTU")
4Count amino acids,                              aa#|sequence|count("ACDEFGHIKLMNPQRSTVWYacdefghiklmnpqrstvwy")
5Count poison characters (RNA/DNA),              poison#|sequence|len;count("acgtumrwsykvhdbnACGTUMRWSYKVHDBN-.")|minus
6Count poison characters (PRO),                  poison#|sequence|len;count("abcdefghiklmnpqrstvwxyzABCDEFGHIKLMNPQRSTVWXYZ-.*")|minus
7Replace poison characters with gaps (RNA/DNA),  |translate("acgtumrwsykvhdbnACGTUMRWSYKVHDBN-.","acgtumrwsykvhdbnACGTUMRWSYKVHDBN-.","-")
8Replace poison characters with gaps (PRO),      |translate("abcdefghiklmnpqrstvwxyzABCDEFGHIKLMNPQRSTVWXYZ-.*","abcdefghiklmnpqrstvwxyzABCDEFGHIKLMNPQRSTVWXYZ-.*","-")
9start gene,                                     start#|sequence|filter(SAI=5_prime,exclude=x)|count("ACGTUMRWSYKVHDBNacgtumrwsykvhdbn")
10end gene,                                       end#|sequence|filter(SAI=3_prime,exclude=x)|count("ACGTUMRWSYKVHDBNacgtumrwsykvhdbn")
11calculate G+C content,                          GC#|sequence|count("gcGC");count("agctuAGCTU")|"G+C:";per_cent
12calculate G+C content helical parts only,       GCH#|sequence|filter(SAI=HELIX;exclude=.-)|count("gcGC");count("agctuAGCTU")|"GCH=";per_cent;"(";cut(2);")"
13calculate G+C content non helical parts only,   GCN#|sequence|filter(SAI=HELIX;include=.-)|count("gcGC");count("agctuAGCTU")|"GCN=";per_cent;"(";cut(2);")"
14
Note: See TracBrowser for help on using the repository browser.