Skip to content

Latest commit

 

History

History
34 lines (24 loc) · 2.32 KB

README.md

File metadata and controls

34 lines (24 loc) · 2.32 KB

primate_numt

all scripts used for primate numts

Insertion_flank_read_extract.pl is a perl script to extract flanking reads around the Numt insertion breakpoint and also its mate pair. This can then be used for doing a de novo assembly of the Numt insertion.

Usage: Insertion_flank_read_extract.pl (File with Numts and samples they are found in bedFormat.txt) (directory with BAM files)

Where you first enter the (name of the script) followed by the (input file in bed format) (the input file has chromosome breakpoint (for any insertion) start coordinate and end coordinate and all the samples the insertion was discovered in saperated by a comma(,)). Followed by the (directory path) where all the BAM files are stored. Please look at the example for input file <bonobo_polymorphic_numts.txt>

Example of input file - bonobo_polymorphic_numts.txt

Hotspot.R is an R script that is able map hotspot and do a permuation analysis to see if the hotspots are significant.

Usage: Open the script and change names of three files <Reference_events.txt> <Polymorphic_events.txt> <Human_karyotype.txt>. 
You can then run different parts of the script based on your requirements (commented in the script). 

Example of input files - 1) Reference_events.txt 2) Polymorphic_events.txt 3) Human_karyotype.txt

Citation