|
0
|
1 #+TITLE: RepeatExplorer based Assembly Annotation Pipeline
|
|
|
2
|
|
|
3 * Tools in repository
|
|
|
4 ** Extract Repeat Library from RepeatExplorer Archive
|
|
|
5 (=extract_re_contigs.xml=)
|
|
|
6
|
|
|
7 This toll will extract library of repeats based on RepeatExplorer2 analysis. Library is available as fasta file. Tool also filter out all the contig parts which has read depth and length below threshold. Parts of contigs with read depth below threshold are hardmasker. Contigs with full hardmasking are removed completelly
|
|
|
8
|
|
|
9 ** Format repeat library
|
|
|
10 (=format_repeat_library.xml=)
|
|
|
11
|
|
|
12 This tool append classification of repeats to library of repeats. Type of repeat is then part of sequence name in format:
|
|
|
13
|
|
|
14 ~>sequence_id#classification_level1/classification_level2/...~ this enable to specify classification hierarchy
|
|
|
15 Classification of sequneces in library is provided using =CLUSTER_TABLE.csv= (part of RE2 output)
|
|
|
16
|
|
|
17 This file can then be used for annotation of repeat in your assembly:
|
|
|
18 ** Repeat Annotation
|
|
|
19 (=repeat_annotate_custom.xml=)
|
|
|
20
|
|
|
21 Internally annotation is performed using RepeatMasker search. Output from RepeatMasker is parsed to remove duplicated and overlaping annotations, Conflicts in annotations are resolved using hierarchical classification of repeats provided in custom database.
|
|
|
22 ** TODO Summarize Annotation
|
|
|
23 This tool will create summary table from GFF annotation.
|
|
|
24 * test data
|
|
|
25
|
|
|
26 - ~test_assembly_1.fasta~ with ~test_db_1_satellites.fasta~ (include CLASS followed by double underscore - syntax 1)
|
|
|
27 - ~test_assembly_2.fasta~ with ~test_db_2_RE_repeats.fasta~ (include full hierarchical classification)
|
|
|
28
|
|
|
29
|
|
|
30
|
|
|
31 #+begin_comment
|
|
3
|
32 # create tarball for toolshed:
|
|
2
|
33 tar -czvf ../repeat_annotation_pipeline.tar.gz --exclude test_data \
|
|
|
34 --exclude .git --exclude tmp --exclude hg_repository --exclude .idea --exclude .gitignore .
|
|
1
|
35 #+end_comment
|
|
3
|
36
|
|
|
37
|
|
|
38
|
|
|
39 - TODO RM take only short name of sequences - validate name / adjust |
