Add stats on chimeric transcripts
Ya and Smith did a detailed analysis of chimeras generated during transcriptome assembly:
We should incorporate their script for identifying chimeras from blastx results into
postassemble.py. It would run right after the
nr_annotate stage. The fraction of each type of chimera would be recorded and displayed in the transcriptome and tabular reports.
We could optionally exclude transcripts that are found to be chimeric from future analyses.