Does Minia discard singletons?

Hi everyone,

I am working with some Minia assemblies. These were run using the no-length-cutoff parameter. Normally I discard small contigs, but I like to see what comes out before i do.

My question: when using the no-length-cutoff parameter, does Minia keep everything, even the singletons, reads which were not part of any contig?

Just curious.

Thanks for any information.

Stacy

minia

• 1,939 views • link

updated 8.9 years ago by Rayan Chikhi

• written 9.5 years ago by stacy734 •

1 answer

Hi, yes Minia should keep singleton k-mers that were not filtered out by the abundance threshold. You can check this by looking at the length of the smallest contigs.

• 0 views • link

written 8.9 years ago by Rayan Chikhi

More posts like this

Statistics for bins
written 1.1 years ago by shevch2009 •

Hello everyone, I’m working with several dozen metagenomic bins and need to analyze contig/scaffold length distributions within each bin. Specifically, I’d like to calculate: Individual …
How do I filter the annotated genome ?
written 3.5 years ago by sunnykevin97

Hi, Sequencing data - Illumina NextSeq 500, 150 bp RL and PE reads. I had the masked fish genomes, using BRAKER2 by providing the vertebrate …
Binning fasta sequences by size?
written 8.6 years ago by stacy734 •

Hi everyone, I have a large fasta file and need to bin it by size: everything 200 nt and up in one file, and everything …
Binning Illumina reads by subread?
written 9.3 years ago by stacy734 •

Hi everyone, I have a very large file of genomic sequence, and need to bin them by the presence or absence of a specific 20-mer. …
Minia estimated genome size
written 9.6 years ago by stacy734 •

Hi everyone, I'm experimenting with Minia for genome assembly of a eukaryote with a 500 MB genome size. I've had good luck adjusting the parameters …
14 contigs recieved from SMRTportal using HGAP.3
written 10.5 years ago by komwit.s •

Hi, I try to assembly my Lactobacillus genome in HGAP.3 using SMRTPortal. I set genome size about 3000000 bp and use coverage at 15x. I …
What is "Total Length" in QUAST?
written 10.5 years ago by '

I have done genome assembly on an interleaved fastq file using many different assemblers (Velvet, ABySS, Minia, SPAdes, etc.) and have the "contigs.fasta" file from …
Forum: Behavior of short read assemblers regarding filtering out short contigs
written 11.3 years ago by Rayan Chikhi

Hi Biostar, This is not a question, just a small blog post. While developing a new version of [Minia][1], I was wondering how assemblers filter …
Cap3 Integration Of Velvet And Newbler Assemblies
written 13.7 years ago by tommivat

<p>I am conducting de novo assembly of ~33Mb genome using 454 and Illumina reads. I cannot use MIRA, since I have ~80M Illumina reads (would …
Very Few Reads Mapping Back To Contigs - Plant Transcriptome
written 14.5 years ago by Cerebralrust •

<p>I assembled plant transcriptome 454 data (non normalised) using trinity after the following</p> <p>1)pre processing (removal of adaptors, vector contamination) 2)removal of rRna sequences 3)removal …