Comments (4)
rsharris
commented on August 10, 2024
Looking at the resulting blocks, it appears that they are consistent with the fact that the reference sequence, PT.PT1S09221, appears on the reverse strand.
In other words, for the first block, I believe it *does* contain position 570 along the forward strand of PT.PT1S09221.
I was initially confused about why it reports 127 as component.start. But working this through:
s PT.PT1S09221 117 12 - 698 GCTGAAGTTGGC
is interval (117,129) on the reverse strand (using zero-based half-open coordinates). That corresponds to GCCAACTTCAGC on the forward strand in (569,581). The request for position 570 would be the second position in that interval. Converting that position as a starting position on the reverse strand is 127, right? I.e. the sub-interval starting at 127 is (127,129), which is the final two characters, e.g. the aligned sequence is split into GCTGAAGTTG/GC and GC is reported.
Bob H
…On Apr 17, 2020, at 3:56 AM, Daniel Lang ***@***.***> wrote:
Hi,
either I'm doing something wrong or I've stumbled over a peculiar bug in the MAF indexing code.
My setup:
Linux/Python 3.7.1/Anaconda
bx-python=0.8.8 installed from -c bioconda
I use a multi-species, multi sequence maf file indexed via:
maf_build_index.py test.maf test.idx
I want to study SNPs locations on the MSA. I have created a small example repo that contains a test set and an ipynb to demonstrate the bug:
https://github.com/dandaman/bxmafissue/blob/master/BugTracer.ipynb <https://github.com/dandaman/bxmafissue/blob/master/BugTracer.ipynb>
As you can see in the notebook, get (or probably the index) is somehow mixing up positions. In my full data set some positions are wrong some are not?!
Am I doing something wrong or is this a bug?
This old report https://www.biostars.org/p/122420/ <https://www.biostars.org/p/122420/> seems to have encountered the same issue.
All the best,
Daniel
—
You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub <#61>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AAFDOI6Y6O3H6COBOIT6GB3RNADUFANCNFSM4MKRBIFQ>.
from bx-python.
rsharris
commented on August 10, 2024
Continuing that thought ... from what you've reported it appears that the lookup always considers the position to be on the forward strand, and the resulting start position is reported along whatever strand the alignment is on.
I wonder if that is documented anywhere.
from bx-python.
dandaman
commented on August 10, 2024
You're right. I did forget about how MAF represents alignments on the negative strand :-(
It's mentioned in the documentation of slice_by_component and I overread it.
I solved my issue by checking the returned alignment strand for my target sequence and the reverse complementing it.
I've updated my example repo and keeping it in case anyone falls into the same trap I have.
Updated the master of my repo accordingly. https://github.com/dandaman/bxmafissue/blob/master/BugTracer.ipynb
If you're not very familiar with MAF it might be a common misconception. Dealing with GFF3 or VCF you'd expect different. Maybe also add this to the doc of get() and point to the MAF standard description at ucsc or so.
From my side we can close this issue, Bob. Thanks for your explanation and for your work on bx-python!
from bx-python.
rsharris
commented on August 10, 2024
Reverse-strand treatment, as well as origin-zero vs origin-one, pop up over and over on forums. Reverse strand in MAF bit me early on, probably in 2004.
It would be helpful if the situation you encountered was in some bx-python doc, e.g. if there is an example of how to use maf index lookup it ought to include an example where the 'hit' is on the reverse strand.
I haven't been a contributor to bx-python in over a decade. And my contribution was pretty small to begin with. I haven't been active in the project as it has matured and improved, and with James' death I don't know any of the people who have done that work. But I do occasionally respond to issues when I think I can understand them.
from bx-python.
Related Issues (20)
- Incompatible types str and byte under python3 HOT 1
- Wheel for Python 3.8 HOT 2
- StopIteration transformed to RuntimeError on python >= 3.7
- No source tarballs for the release 0.8.8 on the PyPI website HOT 1
- error can't found bx.intervals module
- Importing ABC directly from collections was deprecated and will be removed in Python 3.10. Use collections.abc
- time.clock has been removed in Python 3.8
- ImportError: undefined symbol: PyUnicodeUCS2_FromStringAndSize HOT 1
- binned_bitsets_from_list errors when chromosome size is larger than set MAX (512M) HOT 1
- Got wrong positions on chomped query maf subset HOT 2
- How to dealing with maf files compressed with bzip2 or lzop ? HOT 4
- Failed to load seekbzip2 module when dealing with bz2 files HOT 1
- undefined symbol: PyUnicodeUCS2_FromStringAndSize HOT 3
- 0.8.11 Doesn't install some built shared libraries HOT 8
- documentation? HOT 3
- Add support to release linux aarch64 wheels HOT 1
- multiple sequence alignment
- upgrading bx-python to newer versions of python HOT 1
- ERROR Missing dependencies: oldest-supported-numpy HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from bx-python.