Skip to content

Extra scripts

Pavel V. Dimens edited this page Apr 28, 2021 · 20 revisions

extract_markers.sh

This simple script extracts contig and position information from the data_f.call.gz file generated from running the first steps of the pipeline (ParentCall2 and Filtering2) into the file snps.txt.

lepmap2anchor

This is a simple bash script that concatenates the distance intervals files in intervals into a single file necessary to feed into the Lep-Anchor wrapper lepanchor_wrapper.sh (part of Lep-Anchor, not included). It should be run in the intervals folder and the only argument is the markers file that can be obtained from extract_markers.sh.

lepmap2allmaps

This is a python script included in the repository but not used by LepMak3r. The script converts Lep-Map3 distances output into ALLMAPS format. Run it without arguments to see the usage text. Note: LepMak3r does not run ALLMAPS, so please explore the documentation for ALLMAPS found here.

popmap2pedigree

If you have a simple single-family cross, you can use popmap2pedigree to convert your Stacks or dDocent generated popmap file into a LepMap3-compliant pedigree file. Make sure your popmap is a tab-separated two-column text file with the first column of sample names and the second column of parent or progeny designations:

sample_1<tab>parent
sample_2     parent
f1_001     progeny
f1_002     progeny
f1_003     progeny
etc...

So long as the parents are designated parent, it doesn't matter what progeny are designated as. Go wild 😃

(... but like, not too wild)

Getting started

LepMap3 Workflow

LepAnchor Workflow

Other

Clone this wiki locally