6
6
copy (full_name) to this field, :*=*(full_name)
7
7
append (full_name) to this field, :*=*1 *(full_name)
8
8
calculate sequence checksum, seqcheck#|sequence|"ARB_";checksum(exclude=.-n;toupper)
10
count ambiguities (nry. ...) (filter: ECOLI), ambig#|sequence|filter(SAI=ECOLI,exclude=.-)|mid(10,-10)|len("-ACGTU")
9
count ambiguities (nry. ...) (filter: ECOLI), ambig#|sequence|filter(SAI=ECOLI,exclude=.-)|len("-ACGTU")
11
10
truncate seq. ranges outside ecoli range, |filter(SAI=ECOLI,exclude=.)
12
11
increase G+C content in helical regions by 5%, ali_*/data#|change(SAI=HELIX,exclude=.-,change=5,to=GC)
13
12
store taxonomy, ARB_treetax#|taxonomy(1)