BOL: Related items

My commonly used commands in Bioinformatics

Rahul Nayak — Thu, 26 Jul 2018 04:58:45 -0500

FYI, I've found it useful to use MUMmer to extract the specific changes that Racon makes, so I can evaluate them individually:

minimap -t 24 assembly.fasta long_reads.fastq.gz | racon -t 24 long_reads.fastq.gz - assembly.fasta racon_assembly.fasta
nucmer -p nucmer assembly.fasta racon_assembly.fasta
show-snps -C -T -r nucmer.delta

This reports Racon's changes in a table. You can exclude indels with the -I option in show-snps.

This process (Racon -> MUMmer -> SNP table) solves the problem I originally raised in this issue. So as far as I'm concerned, you can close this issue (or keep it open if you still want to implement some kind of variant table).

Referee: Genome assembly quality scores

Jit — Sun, 04 Nov 2018 16:44:30 -0600

Modern genome sequencing technologies provide a succint measure of quality at each position in every read, however all of this information is lost in the assembly process. Referee summarizes the quality information from the reads that map to a site in an assembled genome to calculate a quality score for each position in the genome assembly.

We accomplish this by first calculating genotype likelihoods for every site. For a given site in a diploid genome, there are 10 possible genotypes (AA, AC, AG, AT, CC, CG, CT, GG, GT, TT). Referee takes as input the genotype likelihoods calculated for all 10 genotypes given the called reference base at each position.

Referee is a program to calculate a quality score for every position in a genome assembly. This allows for easy filtering of low quality sites for any downstream analysis.

https://github.com/gwct/referee

Address of the bookmark: https://gwct.github.io/referee/#

jackalope: A swift, versatile phylogenomic and high-throughput sequencing simulator

Abhimanyu Singh — Fri, 26 Jul 2019 00:58:12 -0500

jackalope simply and efficiently simulates (i) variants from reference genomes and (ii) reads from both Illumina and Pacific Biosciences (PacBio) platforms. It can either read reference genomes from FASTA files or simulate new ones. Genomic variants can be simulated using summary statistics, phylogenies, Variant Call Format (VCF) files, and coalescent simulations—the latter of which can include selection, recombination, and demographic fluctuations. jackalope can simulate single, paired-end, or mate-pair Illumina reads, as well as reads from Pacific Biosciences These simulations include sequencing errors, mapping qualities, multiplexing, and optical/PCR duplicates. All outputs can be written to standard file formats.

A swift, versatile phylogenomic and high-throughput sequencing simulator https://jackalope.lucasnell.com

Address of the bookmark: https://github.com/lucasnell/jackalope

Kevler: Reference-free variant discovery in large eukaryotic genomes

Jit — Tue, 28 Jan 2020 03:21:53 -0600

Welcome to kevlar, software for predicting de novo genetic variants without mapping reads to a reference genome! kevlar's k-mer abundance based method calls single nucleotide variants (SNVs), multinucleotide variants (MNVs), insertion/deletion variants (indels), and structural variants (SVs) simultaneously with a single simple model.

More at https://kevlar.readthedocs.io/en/latest/

https://www.cell.com/iscience/pdf/S2589-0042(19)30259-7.pdf

Address of the bookmark: https://github.com/kevlar-dev/kevlar

Caretta – A multiple protein structure alignment and feature extraction suite

Rahul Nayak — Fri, 18 Dec 2020 02:09:44 -0600

Caretta – a multiple protein structure alignment and feature extraction suite

Caretta, a multiple structure alignment suite meant for homologous but sequentially divergent protein families which consistently returns accurate alignments with a higher coverage than current state-of-the-art tools. Caretta is available as a GUI and command-line application and additionally outputs an aligned structure feature matrix for a given set of input structures, which can readily be used in downstream steps for supervised or unsupervised machine learning.

Address of the bookmark: http://www.bioinformatics.nl/caretta/

Genomic Scientist at UDSC

Thu, 28 May 2015 19:14:23 -0500

Centre for Genetic Manipulation of Crop Plants

Department of Genetics

University of Delhi South Campus

NEW DELHI – 110 021

WALK-IN-INTERVIEW FOR THE TEMPORARY POSITIONS OF RESEACH SCIENTIT & LAB / FIELD ATTENDANT

1 Research Scientist (RS) – 3

DBT, Ph. D.

Experience on DNA Markers, plant genome mapping and bioinformatics

Salary: 60,000 (Consolidated) + 5% annual increment

Date and time: 25.06.2015 at 10:30 AM

These temporary positions have been sanctioned in a DBT funded project for the Phase II on ‘Centre of Excellence on genome mapping and molecular breeding of Brassicas.’

The applicants are requested to register their names on the day of interview in the First Floor, Biotech Centre, Centre for Genetic Manipulation of Crop Plants, Department of Genetics before the stipulated time for the interview. Only the registered eligible candidates will be interviewed on the day in the Committee Room.

Applicants are requested to bring all related documents, in original and a set of photocopy, for verification.

No TA/DA will be paid for attending the interview.

www.du.ac.in/du/index.php?mact=News,cntnt01,detail,0&cntnt01articleid=5492&cntnt01returnid=83

GRIDSS: the Genomic Rearrangement IDentification Software Suite

Rahul Nayak — Sun, 17 May 2020 10:27:44 -0500

GRIDSS is a module software suite containing tools useful for the detection of genomic rearrangements. GRIDSS includes a genome-wide break-end assembler, as well as a structural variation caller for Illumina sequencing data. GRIDSS calls variants based on alignment-guided positional de Bruijn graph genome-wide break-end assembly, split read, and read pair evidence.

Address of the bookmark: https://github.com/PapenfussLab/gridss

PeGAS: A Comprehensive Bioinformatic Solution for Pathogenic Bacterial Genomic Analysis

LEGE — Mon, 01 Sep 2025 01:18:10 -0500

This is PeGAS, a powerful bioinformatic tool designed for the seamless quality control, assembly, and annotation of Illumina paired-end reads specific to pathogenic bacteria. This tool integrates state-of-the-art open-source software to provide a streamlined and efficient workflow, ensuring accurate insights into the genomic makeup of pathogenic microbial strains.

Address of the bookmark: https://github.com/liviurotiul/PeGAS

kraken: A universal genomic coordinate translator for comparative genomics

Jit — Thu, 07 Dec 2017 04:45:43 -0600

If you planning on conducting a study involving dozens of large genomes, then you do not have to run all pairwise synteny alignments .. simply try kraken: A universal genomic coordinate translator for comparative genomics

Address of the bookmark: https://github.com/nedaz/kraken

TACOA: Taxonomic classification of environmental genomic fragments using a kernelized nearest neighbor approach

Poonam Mahapatra — Tue, 15 May 2018 09:52:28 -0500

TACOA is a software that can accurately predict the taxonomic origin of genomic fragments from metagenomic data sets by combining the advantages of the k -NN approach with a smoothing kernel function. TACOA can be easily installed and run on a desktop computer, therefore allowing researchers to locally analyze their metagenomic sequence data or integrate it into their pipelines.

Address of the bookmark: http://www.cebitec.uni-bielefeld.de/index.php/2-uncategorised/99-tacoa