InSilicoSeq icon indicating copy to clipboard operation
InSilicoSeq copied to clipboard

Abundance file shows more genomes than are used

Open Naturalist1986 opened this issue 1 year ago • 2 comments

Hi,

I'm using this command to generate metagenomes from a set of genome fasta files:

iss generate --draft (list of *.fna files) --output --model nextseq --coverage uniform -p 128 -n 50M

Although the abundance file shows it uses 10 genomes for example, when I count the reads in the fastqs I only see 5-6 genomes.

Naturalist1986 avatar May 27 '24 09:05 Naturalist1986

Hi!

Are all your contig names unique across draft genomes? InSilicoSeq currently uses the contig names as read headers.

HadrienG avatar Jun 13 '24 13:06 HadrienG

Hi,

Yes, these are downloaded from NCBI Refseq. A portion of the genome contigs is missing from the generated fastqs, even though it is listed in the abundance file.

Naturalist1986 avatar Jun 13 '24 14:06 Naturalist1986