]> Trim poly-A tails off EST sequences (EMBOSS) edit http://bioweb.pasteur.fr/docs/EMBOSS/trimest.html trimest &emboss_init; input Input section sequence sequence -- dna [sequences] (-sequence) perl " -sequence=$value -sformat=fasta" 1 dna 8 seqsfile perl 1 advanced Advanced section minlength Minimum length of a poly-A tail (-minlength) 4 perl (defined $value && $value != $vdef)? " -minlength=$value" : "" 2 This is the minimum length that a poly-A (or poly-T) tail must have before it is removed. If there are mismatches in the tail than there must be at least this length of poly-A tail before the mismatch for the mismatch to be considered part of the tail. 1 mismatches Number of contiguous mismatches allowed in a tail (-mismatches) 1 perl (defined $value && $value != $vdef)? " -mismatches=$value" : "" 3 If there are this number or fewer contiguous non-A bases in a poly-A tail then, if there are '-minlength' 'A' bases before them, they will be considered part of the tail and removed . <BR> For example the terminal 4 A's of GCAGAAAA would be removed with the default values of -minlength=4 and -mismatches=1 (There are not at least 4 A's before the last 'G' and so only the A's after it are considered to be part of the tail). The terminal 9 bases of GCAAAAGAAAA would be removed; There are at least -minlength A's preceeding the last 'G', so it is part of the tail. reverse Write the reverse complement when poly-T is removed (-reverse) 1 perl ($value)? "" : " -noreverse" 4 When a poly-T region at the 5' end of the sequence is found and removed, it is likely that the sequence is in the reverse sense. This option will change the sequence to the forward sense when it is written out. If this option is not set, then the sense will not be changed. fiveprime Remove poly-T tails at the 5' end of the sequence. (-fiveprime) 1 perl ($value)? "" : " -nofiveprime" 5 If this is set true, then the 5' end of teh sequence is inspected for poly-T tails. These will be removed if they are longer than any 3' poly-A tails. If this is false, then the 5' end is ignored. output Output section outseq outseq (-outseq) outseq.out perl " -outseq=$value" 6 seqsfile perl 1 outseq_sformat Output format for: outseq fasta fasta gcg gcg phylip phylip embl embl swiss swiss ncbi ncbi nbrf nbrf genbank genbank ig ig codata codata strider strider acedb acedb staden staden text text fitch fitch msf msf clustal clustal phylip phylip phylip3 phylip3 asn1 asn1 fasta perl " -osformat=$value" 7 auto perl " -auto -stdout" 8