Take a look at the PRICE software from the DeRisi lab. Its meant to do something very similar. http://derisilab.ucsf.edu/index.php?page=software

You could also look at SSPACE (http://www.baseclear.com/landingpages/basetools-a-wide-range-of-bioinformatics-solutions/sspacev12/), ATLAS tools (http://www.hgsc.bcm.tmc.edu/content/bcm-hgsc-software), and SCARPA (http://compbio.cs.toronto.edu/hapsembler/scarpa.html).

See the PAGIT protocol: http://www.sanger.ac.uk/resources/software/pagit/

In particular, take a look at the IMAGE tool: http://genomebiology.com/2010/11/4/R41

Also SOAPdenovo has ha function for scaffolding. Not sure about ABYSS

Here there is a useful explanation of several tools.

https://bioinformaticsonline.com/search?q=scaffolding&entity_type=object&entity_subtype=bookmarks&offset=0&search_type=entities

I could be wrong, but the above answers to your hypothetical scenario appear to miss the point that you aren't interested in assembling the full genome, just the 100 kb part you're interested in. I suggest the following algorithm:

1. Start with the initial assembly C0 of the contigs you have identified as overlapping your region of interest, and the set S of reads those contigs contain. Let C = C0.

2. Repeat:
a. Identify paired-end reads (not in C) for which one or both ends align within, or extending, contigs in C.
b. Identify unpaired reads that align extending these new paired-end reads.
c. Construct a new assembly C' from C and the new reads identified in (a) and (b).
d. Trim C' so it does not extend more than 100 kb to either end of C0. Set C = C'.
e. Let S' denote the reads that contribute to C'. If S' does not contain any reads not present in S, stop. Otherwise, Set S = S'.

3. If you don't have a complete assembly of the region of interest, generate an STS for each end of each contig, probe a library for clones including these STSes, subclone these clones into a paired-end sequencing vector, and generate paired-end reads for this library; then try steps (1) and (2) again, adding these new sequencing reads to what you had before.

4. If your average sequencing depth for the region of interest exceeds 25 or so without filling all gaps, it is likely that the remaining gaps represent sequences that are not getting cloned in your sequencing vectors. Try different sequencing vectors.

Current projects in the Morota lab include developing kernel-based whole-genome prediction and kernel-based genome-wide association models, polygenic modeling of binary traits, reexamining the results from quantitative genetics analysis in light of functional annotation, and extending kernel methods (such as GBLUP and RKHS) specifically tailored for diverse types of emerging omics data.

In addition, candidates will be expected to leverage opportunities to interact with faculty in animal genetics and biometrics at the UNL in the areas of bioinformatics, breeding, functional genomics, quantitative genetics, and molecular genetics.

Candidates should have a B.S. or M.S. degree in quantitative disciplines with strong background and interest in statistical computing.
The starting date is Fall 2015.
For more information about research in the Morota lab at the UNL, visit: http://www.morotalab.org

A letter of interest in the position, C.V., and contact information for
three references should be emailed to Gota Morota at .
Review of applications will begin immediately, and continue until the
positions are filled. Informal inquiries are also welcome.

Also, please see: http://animalscience.unl.edu/anscprospectivegraduatestudents

Address of the bookmark: https://www.molinspiration.com/

Essential qualifications: First class M. Sc. in any branch of life sciences and qualified CSIR-UGC NET/GATE Examination.

Desirable qualifications: Practical experience in biochemical and biophysical studies of proteins

Emoluments: as per DBT norms

The project is tenable for two years, initially for one year.

Age: Below 28 years (relaxable in the case of SC/ST/OBC/women candidates)

Candidates are requested to bring two sets of complete applications on plain paper furnishing bio-data and copies of attested certificates along with originals (for verification) on the date of interview.

No TA/DA is admissible for candidates appearing at the interview.

Dr. Rajat Banerjee
Assistant Professor
Department of Biotechnology and
Dr. B. C. Guha Centre for Genetic Engineering and Biotechnology
University College of Science
35, Ballygunge Circular Road
Kolkata 700019

Advertisement: www.caluniv.ac.in/news/jrf_biotech_2.pdf

More at https://dahak-metagenomics.github.io/dahak/

Address of the bookmark: https://github.com/dahak-metagenomics/dahak

Required Skills: no special skills required for this job post
Required Experience:
Experience in computer aided drug design and or biochemical testing of natural or synthetic compounds is desired
Required Education:
M.Sc. / M.Tech.

Required Job Profile:
Candidate must possess M.Sc. in bioinformatics or computational biology or biotechnology or any branch of life sciences or pharmacology or chemical sciences or M.Tech. in any branch of life sciences with at least fifty five percent marks with NET or GATE.

Desired Job Profile:
Candidate having experience in computer aided drug design and or biochemical testing of natural or synthetic compounds.

How to apply:
Eligible and interested candidates should need to appear for walk-in interview on 21.11.2014 at 1700 hrs at the above mentioned address.

Contact
Pondicherry University
Dr. S. Mohane Coumar, Assistant Professor & Project Investigator, Centre for Bioinformatics, Pondicherry University, Puducherry 605 014
Email:registrar@pondiuni.edu.in
Phone: 0413-2655175

More at http://www.pondiuni.edu.in/sites/default/files/JRF-bioinfor-041114.pdf

At present, some people have posted articles about the analysis process of WGD. I searched for the keyword "wgd pipeline" and found the following:

GenoDup: https:// github.com/MaoYafei/GenoDup-Pipeline
https://peerj.com/articles/6303/
WGDdetector: https:// github.com/yongzhiyang2 012/WGDdetector
https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-019-2670-3
wgd: https:// github.com/arzwa/wgd
https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-016-1142-2#Sec1
https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-017-0399-x
GeNoGAP https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-016-1142-2
https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-017-0399-x
https://github.com/dfguan/purge_dups
https://www.biorxiv.org/content/10.1101/2020.01.24.917997v1

This article introduces the usage of wgd.

Wgd cannot be installed directly with bioconda at present, so it is a little troublesome to install, because it depends on a lot of software. wgd depends on the following software

BLAST
MCL
MUSCLE/MAFFT/PRANK
PAML
PhyML/FastTree
i-ADHoRe

But the good news is that most of the software it depends on can be installed with bioconda

conda create -n wgd python=3.5 blast mcl muscle mafft prank paml fasttree cmake libpng mpi=1.0=mpich
conda activate wgd

Here mpi=1.0=mpich is selected, because i-adhore depends on mpich. If openmpi is installed, an error will appear while loading shared libraries: libmpi_cxx.so.40: cannot open shared object file: No such file or directory

After that, the installation is much simpler

git clone https://github.com/arzwa/wgd.git
cd wgd
pip install .
pip install git+https://github.com/arzwa/wgd.git
For i-ADHoRe, you need to register at http:// bioinformatics.psb.ugent.be /webtools/i-adhore/licensing/Agree to the license to download i-ADHoRe-3.0

Since my miniconda3 installed ~/opt/, the installation path is so~/opt/miniconda3/envs/wgd/

tar -zxvf i-adhore-3.0.01.tar.gz
cd i-adhore-3.0.01
mkdir -p build && cd build
cmake .. -DCMAKE_INSTALL_PREFIX=~/opt/miniconda3/envs/wgd/
make -j 4
make insatall

Take the sugarcane genome Saccharum spontaneum L as an example. The genome is 8-ploid with 32 chromosomes (2n = 4x8 = 32)

Download the tutorial for CDS and GFF annotation files

mkdir -p wgd_tutorial && cd wgd_tutorial
wget http://www.life.illinois.edu/ming/downloads/Spontaneum_genome/Sspon.v20190103.cds.fasta.gz
wget http://www.life.illinois.edu/ming/downloads/Spontaneum_genome/Sspon.v20190103.gff3.gz
gunzip *.gz

First conda activate wgdstart our analysis environment, and then start the analysis

Step 1 : Use to wgd mclidentify homologous genes in the genome

wgd mcl -n 20 --cds --mcl -s Sspon.v20190103.cds.fasta -o Sspon_cds.out

Step 2 : Use to wgd ksdbuild Ks distribution

wgd ksd --n_threads 80 Sspon_cds.out/Sspon.v20190103.cds.fasta.blast.tsv.mcl Sspon.v20190103.cds.fasta

Step 3 : If the quality of the genome is good, then wgd syncollinearity analysis can be used . It can help us find the collinearity block in the genome and the corresponding anchor point

wgd syn --feature gene --gene_attribute ID \
-ks wgd_ksd/Sspon.v20190103.cds.fasta.ks.tsv \
Sspon.v20190103.gff3 Sspon_cds.out/Sspon.v20190103.cds.fasta.blast.tsv.mcl

 For more reading - There are 9 sub-modules in WGD

• kde: KDE fitting to the Ks distribution
• ksd: Ks distribution construction
• mcl: BLASP comparison of All-vs-ALl + MCL classification analysis.
• mix: Hybrid modeling of Ks distribution.
• pre: preprocess the CDS file
• syn: Call I-ADHoRe 3.0 to use GFF files for collinearity analysis
• viz: draw histogram and density plot
• wf1: Ks standard analysis procedure of the whole genome paranome (paranome), call mcl, ksd and syn
• wf2: Ks standard analysis procedure of one-vs-one homologous gene (ortholog), call wcl and kSD

Applications are invited for the AiCADD fellowship (MHRD, Govt. of India) of University of Kerala.

The terms and conditions of the fellowship is given below:

Ø The AiCADD PhD Fellowship scheme will be available for the students registered for full-time research or intending to register and pursue full time research at SIUCEB in frontier areas of bioinformatics, computational biology, systems biology and closely allied areas with focus on Ayur-Informatics.

Ø The fellowships will be widely announced and open to students irrespective of geographical consideration.

Ø Candidates availing of this fellowship shall not be in receipt of any other fellowships concurrently.

Ø Researchers will be selected on the basis of research aptitude test and personal interview.

Ø Each selected student will be eligible for a monthly fellowship of Rs. 10,000/- for the 1st and 2nd year and Rs. 12,000/- for the 3rd year.

Ø Candidates must register for PhD within one year of joining, failing which the fellowship will have to be remitted back.

Ø Candidates receiving the fellowship shall submit bi-annual reports of progress and the continuation of the fellowship will be based on the evaluation of the same.

Ø Candidates are also required to take up academic duties including teaching upto a maximum of 6 hours per week, as directed by AiCADD Principal Investigator.

Interested candidates may please forward their application along with resume on or before 15th November 2014 in the following address. Principal Investigator, AiCADD Centre, Dept. of Computational Biology and Bioinformatics, University of Kerala, Thiruvananthapuram - 695581.

More at https://sites.google.com/site/centreforbioinformatics/announcements

You can try this here, or try it later on your own data.

Get data

We will use the same Illumina data as we used above:

• illumina_R1.fastq.gz: the Illumina forward reads
• illumina_R2.fastq.gz: the Illumina reverse reads

Assemble

Run Spades:

spades.py -1 illumina_R1.fastq.gz -2 illumina_R2.fastq.gz --careful --cov-cutoff auto -o spades_assembly_all_illumina

• -1 is input file of forward reads
• -2 is input file of reverse reads
• --careful minimizes mismatches and short indels
• --cov-cutoff auto computes the coverage threshold (rather than the default setting, “off”)
• -o is the output directory

Results

Move into the output directory and look at the contigs:

infoseq contigs.fasta

Experience Required: 3-5 years of experience

No of Positions : Multiple

Qualifications: Candidates with minimum qualification as M.Sc Bioinformatics with 3-5 years of experience in Life sciences R&D or Pharma Industry.

Ph.D candidates with research experience in Bioinformatics with publications in International journal and minimum 2 years of industry experience in clinical genomics will be preferred for this position.

Requirement:

1. Must have basic understanding of molecular biology and Genomics.

2. Experience in application development or must have expertise in programming using either of Perl/Python.

3. Experience in statistical programming using R/Bioconductor/Matlab.

4. Strong concept in statistical and mathematical modelling.

5. Experience in designing and developing the bioinformatics pipeline.

6. Must have minimum 2+ years of hands on experience in NSG data analysis such as RNA-Seq,Exome-Seq ,Chip-Seq and downstream analysis.

7. Knowledge in WGS ,WES, Targeted re-sequencing,GWAS and population genomics will be preferred.

8. Must have experience working on opensource software/Framework and commercial software for NGS data analysis and reporting.

9. Should be aware of handling big data and guiding team members on multiple projects simultaneously.

10. Should have experience coordinating with different groups of clinical research scientist for various project requirements.

11. Ability to work as team as well as independently with minimal support.

More at http://www.ocimumbio.com/careers1/