"Total" sequence length doesn't include organelles
Nobody has claimed this yet.
- Dominant language
- Jupyter Notebook
- Stars
- 559
- Forks
- 70
- Avg merge
- 1m
- Merged PRs (30d)
- 2
Description
Before opening an issue, please:
- Make sure you are using the latest version using
datasets --version - Review our documentation
Describe the bug
Hello NCBI !
The assembly GCA_964199945.1 is reported as having a "Total Sequence Length" of 1,327,610,284 bp, but the the Fasta file actually contains 1,328,070,353 bp. The difference is exactly the MT and the plastid.
| assmstats-total-sequence-len | Assembly Stats Total Sequence Length |
|---|
To Reproduce
$ datasets summary genome accession GCA_964199945.1 --as-json-lines | dataformat tsv genome --fields assmstats-total-sequence-len --elide-header
1327610284
Expected behavior
I would expect the "total" sequence length to include everything. I would otherwise call it the length of "nuclear" genome only.
Best regards,
Matthieu
Contributor guide
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start by running the datasets summary genome accession GCA_964199945.1 --as-json-lines | dataformat tsv genome --fields assmstats-total-sequence-len --elide-header command shown in the issue and compare its result with the FASTA contents. Trace how the total is calculated for this assembly, including its MT and plastid sequences. Done when the reported total matches the complete FASTA sequence length.
Written by the indexing model from the issue text.
Assessment
- Domain
- cli, data
- Issue type
- Bug
- Difficulty
- 3/5
- Estimated time
- 1-2 days
- Activity status
- Stale
- Clarity
- Clearly specified
- Newbie friendliness
- 40/100