Tag: minimap2

Using Minimap2 with FMLRC2

Using Minimap2 with FMLRC2 1 Hello all, I am using FMLRC2 (github.com/HudsonAlpha/rust-fmlrc) to correct PacBio reads with Illumina reads for hybrid genome assembly. Since FMLRC2 only corrects reads (does not do any assembly) another program is needed. In the paper published on FMLRC minimap (now succeeded by minimap2, github.com/lh3/minimap2) was…

Continue Reading Using Minimap2 with FMLRC2

[lh3/minimap2] Memory leak when using Python and threads

The program align.py uses mappy to align reads in Python using multiple worker threads. After loading the index the memory usage jumps up quickly to >20Gb and then continues to climb steadily through 40Gb an beyond. This issue was first discovered in bonito and isolated to mappy. The data flow…

Continue Reading [lh3/minimap2] Memory leak when using Python and threads

Variant calls of published already assembled genomes

Variant calls of published already assembled genomes 0 I have a set of short read sequencing for the 172 KB Epstein-barr virus genome. We successfully called our variants using GATK to a reference genome. A publication linked below from a different population compared variants (also from short read sequencing) to…

Continue Reading Variant calls of published already assembled genomes

Overestimation of number of reads from nanopore data (flagstat)

Same issue as mentioned on the minimap2 tool: github.com/lh3/minimap2/issues/236#issue-361097444 For example nanopore reads aligned to the host transcriptome the flagstat output is: 5953480 + 0 in total (QC-passed reads + QC-failed reads) 2961480 + 0 secondary 22696 + 0 supplementary 0 + 0 duplicates 4195469 + 0 mapped (70.47% :…

Continue Reading Overestimation of number of reads from nanopore data (flagstat)

Samtools flagstat

Samtools flagstat 1 I aligned my ONT sequencing run with minimap2, subsequently I filtered the file using samtools view -b -F 256 aln_transcriptome_sorted_6.bam -o filtered_aln_transcriptome_6.bam to end up with primary alignments only. When I run samtools flagstat on the filtered file I get the following output: 3502608 + 0 in…

Continue Reading Samtools flagstat

minimap2 – Github Help

38 9 7 minimap2,SneakySnake:snake: is the first and the only pre-alignment filtering algorithm that works efficiently and fast on modern CPU, FPGA, and GPU architectures. It greatly (by more than two orders of magnitude) expedites sequence alignment calculation for both short and long reads. Described in the Bioinformatics (2020) by…

Continue Reading minimap2 – Github Help

About QUAST

About Short intro Further reading How to cite Brief project history Authors Funding License and Terms of use Short intro QUAST evaluates genome assemblies by computing various metrics, including N50, length for which the collection of all contigs of that length or longer covers at least 50% of assembly length…

Continue Reading About QUAST

A Fast, Memory-Efficient, and Accurate Mechanism to Find Fuzzy Seed Matches

BLEND is a mechanism that can efficiently find fuzzy seed matches between sequences to significantly improve the performance and accuracy while reducing the memory space usage of two important applications: 1) finding overlapping reads and 2) read mapping. Finding fuzzy seed matches enable BLEND to find both 1) exact-matching seeds…

Continue Reading A Fast, Memory-Efficient, and Accurate Mechanism to Find Fuzzy Seed Matches

sniffles failed detect SV on minimap2 aligments

When I use ngmlr the sniffles worked. The coverage it more than 90% The code I sent on the github is exactly what it generated, I don’t think there any error Xu Zhang PhD Postdoctoral Associate, Department of Microbiology and Immunology Weill Cornell Medicine 1300 York Avenue, Box 62 New…

Continue Reading sniffles failed detect SV on minimap2 aligments

#1001716 – src:minimap2: fails to migrate to testing for too long: autopkgtest regression

#1001716 – src:minimap2: fails to migrate to testing for too long: autopkgtest regression – Debian Bug report logs Reported by: Paul Gevers <elbrus@debian.org> Date: Tue, 14 Dec 2021 20:03:02 UTC Severity: serious Tags: bookworm, sid Found in version minimap2/2.17+dfsg-12 Fixed in version minimap2/2.22+dfsg-3 Done: Paul Gevers <elbrus@debian.org> Reply or subscribe…

Continue Reading #1001716 – src:minimap2: fails to migrate to testing for too long: autopkgtest regression

Towards the biogeography of prokaryotic genes

1. Sunagawa, S. et al. Structure and function of the global ocean microbiome. Science 348, 1261359 (2015). PubMed  Google Scholar  2. Zou, Y. et al. 1,520 reference genomes from cultivated human gut bacteria enable functional microbiome analyses. Nat. Biotechnol. 37, 179–185 (2019). CAS  PubMed  PubMed Central  Google Scholar  3. Mohammad,…

Continue Reading Towards the biogeography of prokaryotic genes

Bug#1001716: src:minimap2: fails to migrate to testing for too long: autopkgtest regression

Source: minimap2 Version: 2.17+dfsg-12 Severity: serious Control: close -1 2.22+dfsg-3 Tags: sid bookworm User: release.debian….@packages.debian.org Usertags: out-of-sync Dear maintainer(s), The Release Team considers packages that are out-of-sync between testing and unstable for more than 60 days as having a Release Critical bug in testing [1]. Your package src:minimap2 has been…

Continue Reading Bug#1001716: src:minimap2: fails to migrate to testing for too long: autopkgtest regression

The Biostar Herald for Tuesday, September 14, 2021

The Biostar Herald publishes user submitted links of bioinformatics relevance. It aims to provide a summary of interesting and relevant information you may have missed. You too can submit links here. This edition of the Herald was brought to you by contribution from Istvan Albert, and was edited by Istvan…

Continue Reading The Biostar Herald for Tuesday, September 14, 2021

How to plot the gap distribution of contigs wrt to reference genome

How to plot the gap distribution of contigs wrt to reference genome 0 I have a contigs file that I generated using Minia and then I have used Minimap2 to map these contigs again to the reference genome. Now I want to plot the gap distribution i.e. gap distance of…

Continue Reading How to plot the gap distribution of contigs wrt to reference genome

Pacific Biosciences hiring Bioinformatics Software Engineer in United States

PacBio’s Application Software Group focuses on building solid, strategic value around our core data type – highly accurate, long read sequencing – by producing innovative software that unlocks genomics in ways never seen before. We’re growing an interdisciplinary team of bioinformatic experts to tackle some of the most interesting problems…

Continue Reading Pacific Biosciences hiring Bioinformatics Software Engineer in United States

The Biostar Herald for Friday, September 03, 2021

The Biostar Herald publishes user submitted links of bioinformatics relevance. It aims to provide a summary of interesting and relevant information you may have missed. You too can submit links here. This edition of the Herald was brought to you by contribution from zx8754, Istvan Albert, and was edited by…

Continue Reading The Biostar Herald for Friday, September 03, 2021

Global phylogenomic analyses of Mycobacterium abscessus provide context for non cystic fibrosis infections and the evolution of antibiotic resistance

1. Lee, M.-R. et al. Mycobacterium abscessus complex infections in humans. Emerg. Infect. Dis. 21, 1638–1646 (2015). CAS  PubMed  PubMed Central  Google Scholar  2. Prince, D. S. et al. Infection with Mycobacterium avium complex in patients without predisposing conditions. N. Engl. J. Med. 321, 863–868 (1989). CAS  PubMed  Article  Google…

Continue Reading Global phylogenomic analyses of Mycobacterium abscessus provide context for non cystic fibrosis infections and the evolution of antibiotic resistance

unable to find chromosome in SAM header

featureCounts: unable to find chromosome in SAM header 0 I am using featureCounts to try and create a count table for some RNA-Seq data I collected using an Oxford Nanopore platform. I have .sam files aligned with minimap2, and am running the following command to try to get a count…

Continue Reading unable to find chromosome in SAM header

hisat2 compatibility for long read

hisat2 compatibility for long read 0 Hi, I am trying to align PacBio transcriptome reads against the genome to count the gene number. For pair end read i used the following workflow: # convert gff to gtf /home/software/cufflinks-2.2.1/gffread xxx.gff -T -o xxx.gtf # build index /home/software/hisat2-2.2.1/hisat2_extract_exons.py xxx.gtf > xxx.exon /home/software/hisat2-2.2.1/hisat2_extract_splice_sites.py…

Continue Reading hisat2 compatibility for long read

De novo genome assembly

Forum:De novo genome assembly 0 Howdy, I have recently been tasked as the ‘bioinformatics guy” in my lab and am having trouble with a de novo genome assembly of Mother of a Thousand. I am working with Nanopore reads and have ran my reads through CANU. I have all of…

Continue Reading De novo genome assembly