Pockets Identification

CASTp

Pocket-Finder

PocketPicker

##################################################
[Mon May 7 11:29:18 2018] Reduction...
#file name genome size contigs heterozygous size [%] heterozygous contigs [%] identity [%] possible joins homozygous size [%] homozygous contigs [%]
/usr/lib/python2.7/dist-packages/matplotlib/font_manager.py:273: UserWarning: Matplotlib is building the font cache using fc-list. This may take a moment.
warnings.warn('Matplotlib is building the font cache using fc-list. This may take a moment.')
test/run1/contigs.fa 163897 245 66377 40.50 221 90.20 94.854 0 97520 59.50 24 9.80

##################################################
[Mon May 7 11:29:29 2018] Estimating parameters of libraries...
Aligning 19504 mates per library...
Insert size statistics Mates orientation stats
FastQ files read length median mean stdev FF FR RF RR
test/5000_1.fq.gz test/5000_2.fq.gz 50 4998 4990.20 721.47 0 4674 0 0
test/600_1.fq.gz test/600_2.fq.gz 100 599 598.63 47.68 0 10000 0 0

##################################################
[Mon May 7 11:29:29 2018] Scaffolding...
iteration 1.1: test/run1/contigs.reduced.fa 24 97520 39.355 17 94157 7321 2195 0 29603
19505 pairs. 17302 passed filtering [88.71%]. 1627 in different contigs [8.34%].
1526 pairs. 558 in different contigs [36.57%].
iteration 1.2: test/run1/_sspace.1.1.fa 3 97626 39.344 3 97626 87536 6063 821 87536
19505 pairs. 17607 passed filtering [90.27%]. 182 in different contigs [0.93%].
1077 pairs. 124 in different contigs [11.51%].
iteration 2.1: test/run1/_sspace.1.2.fa 3 97626 39.344 3 97626 87536 6063 821 87536
19505 pairs. 15112 passed filtering [77.48%]. 1295 in different contigs [6.64%].
3417 pairs. 396 in different contigs [11.59%].
iteration 2.2: test/run1/_sspace.2.1.fa 1 99133 39.344 1 99133 99133 99133 2328 99133
19505 pairs. 15152 passed filtering [77.68%]. 0 in different contigs [0.00%].
3398 pairs. 0 in different contigs [0.00%].

##################################################
[Mon May 7 11:29:34 2018] Gap closing...
iteration 1.1: test/run1/scaffolds.fa 1 99133 39.344 1 99133 99133 99133 2328 99133

##################################################
[Mon May 7 11:29:35 2018] Final reduction...
#file name genome size contigs heterozygous size [%] heterozygous contigs [%] identity [%] possible joins homozygous size [%] homozygous contigs [%]
[WARNING] Nothing reduced!
test/run1/scaffolds.filled.fa 99390 1 0 0.00 0 0.00 0.000 0 99390 100.00 1 100.00

##################################################
[Mon May 7 11:29:35 2018] Reporting statistics...
#fname contigs bases GC [%] contigs >1kb bases in contigs >1kb N50 N90 Ns longest
test/contigs.fa 245 163897 40.298 24 117391 3975 233 0 29603
test/run1/contigs.fa 245 163897 40.298 24 117391 3975 233 0 29603
test/run1/contigs.reduced.fa 24 97520 39.355 17 94157 7321 2195 0 29603
test/run1/_sspace.1.1.fa 3 97626 39.344 3 97626 87536 6063 821 87536
test/run1/_sspace.1.2.fa 3 97626 39.344 3 97626 87536 6063 821 87536
test/run1/_sspace.2.1.fa 1 99133 39.344 1 99133 99133 99133 2328 99133
test/run1/_sspace.2.2.fa 1 99133 39.344 1 99133 99133 99133 2328 99133
test/run1/scaffolds.fa 1 99133 39.344 1 99133 99133 99133 2328 99133
test/run1/_gapcloser.1.1.fa 1 99390 39.689 1 99390 99390 99390 2 99390
test/run1/scaffolds.filled.fa 1 99390 39.689 1 99390 99390 99390 2 99390
test/run1/scaffolds.reduced.fa 1 99390 39.689 1 99390 99390 99390 2 99390

##################################################
[Mon May 7 11:29:35 2018] Cleaning-up...
#Time elapsed: 0:00:17.376924

The Mix algorithm, approach and results were published in BMC bioinformatics : http://www.biomedcentral.com/1471-2105/14/S15/S16.

Address of the bookmark: https://github.com/cbib/MIX

gapFinisher can fill gaps in draft genomes quickly and reliably.

Address of the bookmark: https://github.com/kammoji/gapFinisher

CPAT - https://github.com/liguowang/cpat, http://lilab.research.bcm.edu/ (web server)

CPC2 - https://github.com/gao-lab/CPC2_standalone, http://cpc2.gao-lab.org/ (web server)

Infernal - http://eddylab.org/infernal/, https://github.com/EddyRivasLab/infernal

NCBI RefSeq - https://www.ncbi.nlm.nih.gov/refseq/

Rfam - http://rfam.xfam.org/, https://docs.rfam.org/en/latest/index.html

SILVA - https://www.arb-silva.de/

RNAmmer - http://www.cbs.dtu.dk/services/RNAmmer/ (web server, standalone download link)

https://evirusbioinfc.notion.site/18e21bc49827484b8a2f84463cb40b8d?v=92e7eb6703be4720abf17a901bc9a947

Address of the bookmark: https://evirusbioinfc.notion.site/18e21bc49827484b8a2f84463cb40b8d?v=92e7eb6703be4720abf17a901bc9a947

To fully leverage these opportunities, bioinformaticians require specialized tools and databases. This blog highlights essential resources for mycology research, focusing on databases, tools, and platforms tailored for fungal biology.

1. Fungal Databases

1.1. MycoCosm

Website: MycoCosm
Developed by the DOE Joint Genome Institute, MycoCosm is a comprehensive portal for fungal genomics. It offers genomic and transcriptomic data for a wide range of fungi, including saprobes, pathogens, and symbionts.

• Key Features: Genome browsers, comparative genomics tools, and functional annotations.
• Best For: Large-scale studies on fungal evolution and ecology.

1.2. FungiDB

Website: FungiDB
FungiDB is an integrated genomic resource for fungal pathogens and non-pathogens. It provides access to genome sequences, transcriptomic data, and functional annotations.

• Key Features: Advanced search options, BLAST, and pathway analysis tools.
• Best For: Studying fungal pathogenesis and host-pathogen interactions.

1.3. Index Fungorum

Website: Index Fungorum
This nomenclatural database provides information on the scientific names of fungi. It’s an essential resource for taxonomists and researchers focused on fungal biodiversity.

• Key Features: Taxonomic hierarchy and synonymy tracking.
• Best For: Identifying and classifying fungal species.

1.4. UNITE

Website: UNITE
UNITE is a specialized database for fungal ITS (Internal Transcribed Spacer) sequences, often used in fungal identification and phylogenetics.

• Key Features: Curated reference datasets and community annotations.
• Best For: Environmental mycology and microbial ecology studies.

2. Analytical Tools

2.1. Funannotate

Repository: GitHub - Funannotate
Funannotate is a genome annotation tool designed for fungi. It supports tasks like gene prediction, functional annotation, and orthology analysis.

• Best For: Annotating newly sequenced fungal genomes.

2.2. BUSCO (Benchmarking Universal Single-Copy Orthologs)

Website: BUSCO
BUSCO evaluates genome assembly and annotation completeness using orthologs. It includes a fungal-specific dataset.

• Best For: Assessing the quality of fungal genome assemblies.

2.3. Pathogen-Host Interactions Database (PHI-base)

Website: PHI-base
PHI-base is a manually curated resource containing information on pathogen-host interactions, including fungal pathogens.

• Best For: Exploring virulence factors and host-pathogen relationships.

3. Visualization Platforms

3.1. Cytoscape

Website: Cytoscape
A powerful tool for visualizing molecular interaction networks, Cytoscape can be used to study protein-protein interactions, gene networks, and metabolic pathways in fungi.

• Best For: Network biology and functional genomics.

3.2. iTOL (Interactive Tree of Life)

Website: iTOL
iTOL is an interactive tool for visualizing phylogenetic trees.

• Best For: Displaying fungal phylogenies and comparing evolutionary relationships.

4. Community Resources

4.1. Mycological Society of America (MSA)

Website: MSA
The MSA promotes fungal research and provides access to resources, conferences, and publications.

• Best For: Networking with fungal researchers and accessing recent studies.

4.2. OpenFungi

Website: OpenFungi
OpenFungi is an open-source initiative providing fungal genomic and transcriptomic datasets for research and education.

• Best For: Sharing and accessing public fungal datasets.

5. Genomics Workflows

5.1. Galaxy

Website: Galaxy Project
Galaxy offers a web-based platform for reproducible bioinformatics workflows, including tools for fungal genome and transcriptome analysis.

• Best For: User-friendly analysis pipelines without requiring coding skills.

5.2. Snakemake

Repository: Snakemake
A flexible pipeline management tool that supports fungal data processing and analysis.

• Best For: Custom workflows for large-scale fungal datasets.

Conclusion

Fungal research is a rapidly growing field with vast implications for medicine, agriculture, and industry. For bioinformaticians, the availability of specialized resources—databases, tools, and community platforms—opens doors to innovative discoveries. Whether you are investigating fungal genomics, studying host-pathogen interactions, or exploring fungal biodiversity, the resources outlined above will empower your research journey.

Dive into these resources and help unravel the mysteries of the fungal kingdom!

• Bioconductor – A plethora of tools for analysis and comprehension of high-throughput genomic data, including 1500+ software packages. [ paper-2004 | web ]
• Biopython – Freely available tools for biological computing in Python, with included cookbook, packaging and thorough documentation. Part of the Open Bioinformatics Foundation. Contains the very useful Entrez package for API access to the NCBI databases. [ paper-2009 | web ]
• Bioconda – A channel for the conda package manager specializing in bioinformatics software. Includes a repository with 3000+ ready-to-install (with conda install) bioinformatics packages. [ paper-2018 | web ]
• BioJulia – Bioinformatics and computational biology infastructure for the Julia programming language. [ web ]
• Rust-Bio – Rust implementations of algorithms and data structures useful for bioinformatics. [ paper-2016 ]
• SeqAn – The modern C++ library for sequence analysis.

Following are the list of standalone applications for network analysis:

Arena 3D

3D visualization of multi-layer networks

http://www.arena3d.org

Biana

Data integration and network management

http://sbi.imim.es/web/BIANA.php

BioLayout Express 3D 

2D/3D network visualization

http://www.biolayout.org/

BiologicalNetworks 

Efficient integrated multi-level analysis of microarray, sequence, regulatory and other data

http://www.biologicalnetworks.org

BioMiner

Modeling, analyzing and visualizing biochemical pathways and networks

http://www.zbi.uni-saarland.de/chair/projects/BioMiner

Cell Illustrator 

Petri nets for modeling and simulating biological networks

http://www.cellillustrator.com

COPASI

Analysis of biochemical networks and their dynamics

http://www.copasi.org/

Cytoscape 

Network visualization and analysis. Over 200 plugins [60]

http://www.cytoscape.org/

Dizzy

Chemical kinetics stochastic simulation software

http://magnet.systemsbiology.net/software/Dizzy/

DyCoNet

Gephi plugin that can be used to identify dynamic communities in networks

https://github.com/juliemkauffman/DyCoNet

GENeVis 

Network and pathway visualization

http://tinyurl.com/genevis/

GEPHI 

Interactive visualization and exploration for any network and complex system, dynamic and hierarchical graph.

https://gephi.org

Igraph

Collection of network analysis tools with the emphasis on efficiency, portability and ease of use

http://igraph.sourceforge.net

Medusa

Semantic and multi-edged simple networks

https://sites.google.com/site/medusa3visualization/

NAViGaTOR

Visualizing and analyzing protein-protein interaction networks

http://tinyurl.com/navigator1/

N-Browse

Interactive graphical browser for biological networks

http://www.gnetbrowse.org/

NeAT

Topological and clustering analysis of networks

http://rsat.ulb.ac.be/neat/

Ondex 

Data integration and visualization of large networks

http://www.ondex.org/

Osprey

Visualization and annotation of biological networks

http://biodata.mshri.on.ca/osprey/servlet/Index

Pajek 

Analysis and visualization of large networks and social network analysis

http://vlado.fmf.uni-lj.si/pub/networks/pajek/

PathwayAssist 

Navigation and analysis of biological pathways, gene regulation networks and protein interaction maps.

http://www.ariadnegenomics.com/downloads/

PIVOT 

Layout algorithms for visualizing protein interactions and families

http://acgt.cs.tau.ac.il/pivot/

ProCope 

Prediction and evaluation of protein complexes from purification data experiments

http://www.bio.ifi.lmu.de/Complexes/ProCope/

ProViz 

Visualization and exploration of interaction networks. Gene Ontology and PSI-MI formats supported

http://cbi.labri.fr/eng/proviz.htm

SpectralNET 

Network analysis and visualizations. Scatter plots and dimensionality reduction algorithms

https://www.broadinstitute.org/software/spectralnet

Tulip 

Enables the development of algorithms, visual encodings, interaction techniques, data models and domain-specific visualizations

http://tulip.labri.fr/TulipDrupal/

VANESA 

Automatic reconstruction and analysis of biological networks and Petri nets based on life-science database information

http://agbi.techfak.uni-bielefeld.de/vanesa/

VANTED 

Network reconstruction, data visualization, integration of various data types, network simulation

http://tinyurl.com/vanted/

yEd

Creation of diagrams manually and import external data

http://tinyurl.com/yEdGraph/

Web tools for network analysis

APID 

Unified protein-protein interactions from BIND, BioGRID, DIP, HPRD, IntAct and MINT

http://bioinfow.dep.usal.es/apid/

Arcadia 

Translates text-based descriptions of biological networks (SBML files) into standardized diagrams (Systems Biology Graphical Notation Process Description maps)

http://arcadiapathways.sourceforge.net/

AVIS 

Viewer for signaling networks

http://actin.pharm.mssm.edu/AVIS2

bioPIXIE 

Discovery of biological networks from diverse functional genomic data

http://pixie.princeton.edu/pixie

CellPublisher

Interactive representations of biochemical processes

http://cellpublisher.gobics.de/

Graphle

Distributed network exploration and visualization of interactive large, dense graphs

http://tinyurl.com/graphle/

GraphWeb 

Web server for graph-based analysis of biological networks

http://biit.cs.ut.ee/graphweb/

Hubba

Web-based service to explore the essential nodes in a network

http://hub.iis.sinica.edu.tw/Hubba

NetworkBLAST 

Analysis of protein interaction networks across species to infer protein complexes that are conserved in evolution

http://www.cs.tau.ac.il/~bnet/networkblast.htm

Pathview 

Tool set for pathway-based data integration and visualization

http://Pathview.r-forge.r-project.org/

PINA 

Integrated platform for protein interaction network construction, filtering, analysis, visualization and management

http://cbg.garvan.unsw.edu.au/pina/home.do

ReMatch 

Web-based tool for integration of user-given stoichiometric metabolic models into a database collected from public data sources

http://www.cs.helsinki.fi/group/sysfys/software/rematch/

SNOW 

Gene mapping on a reference or human protein-protein interaction network that SNOW hosts

http://snow.bioinfo.cipf.es

STITCH 

Resource to explore known and predicted interactions of chemicals and proteins

http://stitch.embl.de/

STRING

Protein interaction networks and integration of data such as genomic context, high-throughput experiments, conserved coexpression and previous knowledge derived from the literature

http://string-db.org

TVNViewer 

An interactive visualization tool for exploring networks that change over time or space

http://www.sailing.cs.cmu.edu/main/?page_id=545

tYNA 

System for managing, comparing and mining multiple networks

http://tyna.gersteinlab.org/tyna/

VisANT 

Visualization, mining, analysis and modeling of biological networks, metabolic networks and ecosystems

http://visant.bu.edu/

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.