BOL: Related items

UniqueKmer: Generate unique KMERs for every contig in a FASTA file

Abhi — Fri, 17 Dec 2021 00:08:15 -0600

Generate unique k-mers for every contig in a FASTA file.

Unique k-mer is consisted of k-mer keys (i.e. ATCGATCCTTAAGG) that are only presented in one contig, but not presented in any other contigs (for both forward and reverse strands).

This tool accepts the input of a FASTA file consisting of many contigs, and extract unique k-mers for each contig.

The output unique k-mer file and Genome file can be used for fastv: https://github.com/OpenGene/fastv, which is an ultra-fast tool to identify and visualize microbial sequences from sequencing data.

https://github.com/OpenGene/UniqueKMER

Address of the bookmark: https://github.com/OpenGene/UniqueKMER

School of Life Sciences, Jawaharlal Nehru University vacancy of JRF / SRF / RA in CSIR funded Project

Wed, 29 Apr 2015 21:26:19 -0500

School of Life Sciences, Jawaharlal Nehru University has issued notification dated 27.04.2015 to fill the vacancy of JRF / SRF / RA in CSIR funded Projec entitled "Structural and functional characterization of serine biosynthetic pathway enzymes from entamoeba histolytica". It is good chance to get job with IITKGP and brighten your future. Learn eligibility criteria and apply on or before 08.05.2015.

Employer: Jawaharlal Nehru University
Address: Dr. S. Gourinath, Principal Investigator, School Of Life Sciences, Jawaharlal Nehru University, New Delhi-110067
Email: not mentioned / provided for this job post
URL: http://www.jnu.ac.in/Career/currentjobs.htm
Phone: 011 2674 2575
Skills: not mentioned / required for this job post
Experience: Experience in molecular biology, structural biology and bioinformatics is desired
Education: M.Sc. in any field of life sciences.
Job Location: New Delhi, Delhi, India (View Jobs in New Delhi, Jobs in Delhi, Jobs in India)

Job Description: School of Life Sciences, Jawaharlal Nehru University vacancy of JRF / SRF / RA in CSIR funded Projec

Name of the Post: JRF / SRF / RA

Salary: As per rules

Required Job Profile:

Candidate must possess M.Sc. in any field of life sciences.

Desired Job Profile:

Candidate having NET - CSIR or UGC and experience in molecular biology, structural biology and bioinformatics is desired and experience with publication is preferred.

How to apply:

Eligible and interested candidates should need to apply with complete details to the above mentioned address on or before 08.05.2015.

Refer to http://www.jnu.ac.in/Career/currentjobs.htm

KAST

Neel — Wed, 23 Feb 2022 08:28:36 -0600

Perform Alignment-free k-tuple frequency comparisons from sequences. This can be in the form of two input files (e.g. a reference and a query) or a single file for pairwise comparisons to be made.

Address of the bookmark: https://github.com/martinjvickers/KAST

Research Fellows at AIMSCS, Hyderabad

Wed, 06 May 2015 06:23:33 -0500

C.R.Rao Advanced Institute of Mathematics, Statistics and Computer Science (AIMSCS) - Hyderabad, Andhra Pradesh
Advertisement No.: 5/2015

Research Fellows Systems Biology job vacancy in C.R.Rao Advanced Institute of Mathematics, Statistics and Computer Science (AIMSCS)

JRF : Qualification - M. Sc in Bioinformatics, Systems Biology, M. Sc statistics, or M. Tech in Bioinformatics,

Pay Scale : Rs. 25,000

SRF : Qualification- Qualification prescribed for JRF with 2 years of research experience.

Pay Scale : Rs. 28,000*

No.of Post: 2

Desirable: Candidates should have strong background in Computational biology, bioinformatics, statistics and algorithmic development. In addition to that previous experience of working on Linux, bio-informatics, NGS data analysis and Basic knowledge of biology is desirable. Programming on any one of the programming languages (C, C++, perl, python) and statistical framework (e.g. R, matlab, etc.) is highly desirable.

More at http://www.crraoaimscs.org/jrf_application_form_2015.pdf

ALE: Assembly Likelihood Estimator

BioStar — Wed, 08 Mar 2023 01:39:33 -0600

Just import the assembly, bam and ALE scores. You can convert the .ale file to a set of .wig files with ale2wiggle.py and IGV can read those directly. Depending on your genome size you may want to convert the .wig files to the BigWig format.

Address of the bookmark: https://github.com/sc932/ALE

Perl One liner basics !!

Abhimanyu Singh — Sun, 24 May 2015 09:28:33 -0500

Perl has a ton of command line switches (see perldoc perlrun), but I'm just going to cover the ones you'll commonly need to debug code. The most important switch is -e, for execute (or maybe "engage" :) ). The -e switch takes a quoted string of Perl code and executes it. For example:

$ perl -e 'print "Hello, World!\n"'
Hello, World!

It's important that you use single-quotes to quote the code for -e. This usually means you can't use single-quotes within the one liner code. If you're using Windows cmd.exe or PowerShell, you must use double-quotes instead.

I'm always forgetting what Perl's predefined special variables do, and often test them at the command line with a one liner to see what they contain. For instance do you remember what $^O is?

$ perl -e 'print "$^O\n"'
linux

It's the operating system name. With that cleared up, let's see what else we can do. If you're using a relatively new Perl (5.10.0 or higher) you can use the -E switch instead of -e. This turns on some of Perl's newer features, like say, which prints a string and appends a newline to it. This saves typing and makes the code cleaner:

$ perl -E 'say "$^O"'
linux

Pretty handy! say is a nifty feature that you'll use again and again.

Steps to find all the repeats in the genome !

Neel — Thu, 31 Aug 2023 02:43:28 -0500

To find repeats in a genome from 2 to 9 length using a Perl script, you can use the RepeatMasker tool with the "--length" option[0]. Here's a step-by-step guide:

Install RepeatMasker: First, you need to install RepeatMasker on your system. You can download it from the RepeatMasker website[0].

Prepare the genome sequence: Make sure you have the genome sequence in a FASTA file format. Let's assume the file is named "genome.fasta".

./RepeatMasker -pa -nolow -norna -no_is -div -lib RepeatMaskerLib.embl -gff -xsmall -small -poly -species -dir -length - genome.fasta

Replace the following placeholders with appropriate values:

: The number of processors/threads you want to use for parallel processing.
: The divergence value for the species you are analyzing. You can find divergence values for different species in the RepeatMasker documentation[0].
: The name of the species you are analyzing.
: The directory where you want the output files to be saved.
and : The minimum and maximum lengths of the repeats you want to find (in this case, 2 and 9).

Analyze the output: RepeatMasker will generate several output files, including a .out file. You can parse this file to extract the information you need. There is a Perl tool called "one_code_to_find_them_all.pl" that can help you parse RepeatMasker output files[0]. You can download it from the source provided.

Use the provided Perl script: Once you have the "one_code_to_find_them_all.pl" script, you can run it to conveniently parse the RepeatMasker output files. Here's an example of how to use it:

perl one_code_to_find_them_all.pl --rm --length

Replace with the path to your RepeatMasker .out file, and with the path to a file containing the lengths of the reference elements.

This script will generate several output files, including .log.txt and .copynumber.csv, which contain quantitative information about the identified repeat elements.

Remember to adjust the parameters and options according to your specific needs and the characteristics of your genome.

Tools to access the quality of your assembled genome !

LEGE — Thu, 08 Aug 2024 23:31:18 -0500

FASTA VALIDATOR + SEQKIT RMDUP: FASTA validation
GENOMETOOLS GT GFF3VALIDATOR: GFF3 validation
ASSEMBLATHON STATS: Assembly statistics
GENOMETOOLS GT STAT: Annotation statistics
NCBI FCS ADAPTOR: Adaptor contamination pass/fail
NCBI FCS GX: Foreign organism contamination pass/fail
BUSCO: Gene-space completeness estimation
TIDK: Telomere repeat identification
LAI: Continuity of repetitive sequences
KRAKEN2: Taxonomy classification
HIC CONTACT MAP: Alignment and visualisation of HiC data
MUMMER → CIRCOS + DOTPLOT & MINIMAP2 → PLOTSR: Synteny analysis
MERQURY: K-mer completeness, consensus quality and phasing assessment

Rosenberg lab

Wed, 27 May 2015 17:52:24 -0500

Research. Research in the lab focuses on mathematical, statistical, and computational problems in evolutionary biology and human genetics. Long-term interests of the lab include topics such as:

Human genetic variation
Inference of human evolutionary history from genetic markers
Statistical analysis of population-genetic data
Mathematical models of gene genealogies
Theoretical population genetics
Combinatorics of evolutionary trees
The relationship between gene trees and species trees
The role of human evolutionary genetics in the search for genes that contribute to disease-susceptibility
More at https://web.stanford.edu/group/rosenberglab/index.html

Genome Simulation with SLiM and msprime

BioStar — Fri, 31 Jan 2025 12:47:43 -0600

Genome simulation is an essential tool in population genetics, enabling researchers to model evolutionary processes and study genetic variation. Two widely used simulation tools in this field are SLiM and msprime. While both serve different purposes, they can be used together with the slendr framework to compare simulation outputs effectively.

Overview of SLiM and msprime

SLiM: Forward Genetic Simulator

SLiM is a free, open-source tool designed for forward genetic simulations. It allows researchers to model complex evolutionary scenarios, including selection, recombination, and demographic events, making it particularly useful for studying adaptation and selection in populations.

Key Features of SLiM:

Simulates population evolution forward in time
Supports custom evolutionary models using an embedded scripting language
Allows modeling of spatial and ecological dynamics
Provides high flexibility and extensibility for user-defined scenarios
Available on GitHub as an open-source project

msprime: Ancestry and Mutation Simulator

msprime is an efficient, open-source tool that simulates ancestry and mutations using a coalescent framework. It is known for its high-speed performance and low memory requirements, making it a popular choice for large-scale genomic simulations.

Key Features of msprime:

Implements coalescent simulations for ancestry modeling
Efficiently simulates large population histories
Supports the addition of mutations to genealogies
Developed using an open-source community model
Often faster and more memory-efficient than alternative simulators

Using SLiM and msprime with slendr

Both SLiM and msprime can be integrated with slendr, a framework that facilitates structured population genetic simulations. This integration allows for seamless comparison of simulation outputs.

How They Work Together:

SLiM and msprime simulations can be analyzed within slendr.
The ts_read() function in slendr enables loading and comparing tree sequence outputs from both simulators.
This integration allows researchers to validate simulation results and gain deeper insights into evolutionary processes.

Performance Considerations

While SLiM offers powerful forward simulations with extensive customization, msprime is often preferred for its speed and memory efficiency when simulating ancestry and mutations. The choice between the two depends on the research goals:

For detailed evolutionary modeling with selection and recombination: Use SLiM.
For large-scale coalescent simulations with mutations: Use msprime.
For comparing different simulation models and their outputs: Use slendr to integrate SLiM and msprime results.

Conclusion

SLiM and msprime are valuable tools for genome simulation, each serving distinct but complementary purposes in population genetics research. By leveraging the strengths of both simulators with slendr, researchers can conduct robust and efficient evolutionary simulations, enhancing our understanding of genetic diversity and adaptation.

For more information, check out the official GitHub repositories for SLiM and msprime, and explore the slendr framework for streamlined simulation workflow