Scripts for simulation and detection of Transposons
Simulate transposon insertions and exsisions with BAMSurgeon and RSVSim and evaluate results
bash run_simulations.sh
bash analyze_simulations.sh
mkdir data
mkdir results
mkdir piRNA
mkdir tables
Run TEMP to detect insertion events relative to the reference and TE-Locate to detect absences or maintenance of transposons relative to the reference. In data directory:
bash run_TransposonCaller.sh
- Filter out redundant calls
- Apply read support and population frequency filters
- Convert transposon names
- Collapse transposons
- Produce consistently formatted output files
bash run_ProcessTransposonCallerOutput.sh
- Check for and resolve contradictory calls
- Merge transposons
- Generate matrix of calls
- Filter monomorphic sites
- Account for sites with no coverage
- Calculate trait metrics
- Check transposon classification and prevalence
- Find intersection with genes
bash CleanTransposons.sh
Obtain mapping results file and analyze GWAS results in context of genes of interest and piRNAs
bash PostMappings.sh