Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Abundance file shows more genomes than are used #262

Open
Naturalist1986 opened this issue May 27, 2024 · 2 comments
Open

Abundance file shows more genomes than are used #262

Naturalist1986 opened this issue May 27, 2024 · 2 comments
Labels
need answer issues that need clarifications from OP

Comments

@Naturalist1986
Copy link

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.

@HadrienG
Copy link
Owner

Hi!

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

@HadrienG HadrienG added the need answer issues that need clarifications from OP label Jun 13, 2024
@Naturalist1986
Copy link
Author

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.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
need answer issues that need clarifications from OP
Projects
None yet
Development

No branches or pull requests

2 participants