mixtureS: a novel tool for bacterial strain reconstruction from reads
mixtureS that can de novo identify bacterial strains from shotgun reads of a clonal or metagenomic sample, without prior knowledge about the strains and their variations. Tested on 243 simulated datasets and 195 experimental datasets, mixtureS reliably identified the strains, their numbers and ...Tags: mixtureS, novel, tool, bacterial, strain, reconstruction, reads
1352 days ago
Hifiasm: a haplotype-resolved assembler for accurate Hifi reads
Hifiasm is a fast haplotype-resolved de novo assembler for PacBio Hifi reads. It can assemble a human genome in several hours and works with the California redwood genome, one of the most complex genomes sequenced so far. Hifiasm can produce primary/alternate assemblies of quality competitive wit...Tags: Hifiasm, haplotype-resolved, assembler, accurate, Hifi, reads
1226 days ago
FastProNGS: fast preprocessing of next-generation sequencing reads
FastProNGS to integrate the quality control process with automatic adapter removal. Parallel processing was implemented to speed up the process by allocating multiple threads. Compared with similar up-to-date preprocessing tools, FastProNGS is by far the fastest.Tags: FastProNGS, fast, preprocessing, next-generation, sequencing, reads
1225 days ago
AlignGraph2: similar genome-assisted reassembly pipeline for PacBio long reads
AlignGraph2 is the second version of AlignGraph for PacBio long reads. It extends and refines contigs assembled from the long reads with a published genome similar to the sequencing genome. More at https://academic.oup.com/bib/advance-article-abstract/doi/10.1093/bib/bbab022/61467...Tags: AlignGraph2, genome-assisted, reassembly, pipeline, PacBio, long, reads
1147 days ago
How to resolve samtools index error "Chromosome blocks not continuous" ?
samtools index mappedPairs.sorted.bam return the following error ... [bam_sort_core] merging from 0 files and 2 in-memory blocks...[E::hts_idx_push] Chromosome blocks not continuoussamtools index: failed to create index for "mappedPairs.sorted.bam"Tags: samtools, index, sort, bam, reads
996 days ago
QuasR: Quantification and annotation of short reads in R
The QuasR package (short for Quantify and annotate short reads in R) integrates the functionality of several R packages (such as IRanges (Lawrence et al. 2013) and Rsamtools) and external software (e.g. bowtie, through the Rbowtie package, and HISAT2, through the Rhisat2 package). The packag...Tags: QuasR, Quantification, annotation, short, reads, R, Tool, Reads
995 days ago
Bioinformatics tools for telomere to telomere assembly !
● Merfin – k-mer-based assembly and variant calling evaluation for improved consensus accuracy (Arang Rhie)● PanGenie – algorithm that leverages a pangenome reference built from haplotype-resolved genome assemblies in conjunction with k-mer count information from r...Tags: optimize, HiFi, reads, correct, tools
990 days ago
Tags: k-mer, counting, learn, teach, skill, genome, reads
990 days ago
Short-read assembly using Spades !
If we only had Illumina reads, we could also assemble these using the tool Spades. 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 reve...Tags: spades, assembly, short, reads, genome, tools
824 days ago
HairSplitter: assembling long reads in an unknown number of haplotypes
Pros and cons of HairSplitter Limitations of HairSplitter: Not very fast: it re-polishes the whole assembly Limited in the number of haplotypes Strengths of HairSplitter: Very modular, can be used with any assembler Naive: makes no assumption on ploidy, parameter-free Safe: won’...Tags: HairSplitter, assembling, long, reads, unknown, number, haplotypes
514 days ago