Skip to content

mills-lab/primate_numt

Repository files navigation

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

About

all scripts used for primate numts

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published