Skip to content

An integrated retroduplication caller based on 1) exon-exon junction and 2) discordant reads

Notifications You must be signed in to change notification settings

gersteinlab/retrodup

Repository files navigation

Retrodup

An integrated retroduplication caller based on 1) exon-exon junction and 2) discordant reads

Please refer to /Landscape and Variation of Novel Retroduplications in 26 Human Populations/.

  1. Create junction libraries

create_junction_library.pl -m null_mode -g genome_dir annotation_file (gencode)

  1. Extract unmapped reads and re-map to junctions

run_bwa.sh -o output_suffix junction_libs bam_files >& out.txt

  1. Call retroduplications from exon-exon junction

call_retroduplications.pl annotation_file sam_file > result.txt

  1. Cluster the discordant reads to discover insertion sites

samtools view bamfile gene_coor | java Cluster gene_coor2 > insertionpoints.txt

About

An integrated retroduplication caller based on 1) exon-exon junction and 2) discordant reads

Topics

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published