Menu

trformat usage

Q:

.. how is trformat.pl used ?

A:
evigene/scripts/rnaseq/trformat.pl
is used BEFORE tr2aacds to regularize IDs in fasta of Velvet,Soap,Trinity, ensure unique IDs, add prefixes for parameter sets. It is is mainly to ensure unique IDs for input transcripts, where assemblers run with varying options (kmer,..) will produce subsets with same ID.

trformat.pl guesses from file names options like kmer and assembly type, then ensures combined transcript set has unique IDs. It does ID reformatting to make those of various assemblers look similar, from these assembler types: velvet, soaptr, trinity, idba, and pacbio (as of 2017 version).

When I run multi-kmer assemblers, each kmer output has subfolder "vel_k32" "vel_k61" .. then IDs of combined assembly look like "velk32g001t1", "velk61g001t1", ...

Here is a script to run trformat.pl on set of velvet, soap multi-kmer assemblies:
http://arthropods.eugenes.org/EvidentialGene/arthropods/mosquito/evg_scripts/trformgroup.sh
where calls are like this

  nap=MySppIDPrefix
  subd=assemblies_folder
  $evigene/scripts/rnaseq/trformat.pl -pre $nap -out $subd.tr -log -in  $subd/vel*/transcripts.fa
  $evigene/scripts/rnaseq/trformat.pl -pre $nap -out $subd.tr -log -in  $subd/sod*/so*.scafSeq.gz

trformat.pl is not essential to using tr2aacds, which checks for duplicate IDs, but you should somehow ensure your input transcript set has unique transcript IDs.

Related

Wiki: Home

Posted by Don Gilbert 2017-10-06

Log in to post a comment.

MongoDB Logo MongoDB