When running genomic_to_tx_segment, it will select the MANE Select transcript, and return this transcript even if it does not exist in UTA. We should add a check to see if the transcript exists in UTA before returning.
Steps to reproduce
An example where this occurs is selecting the MANE select transcript for the gene "NPIPB5"
Expected behavior
I expected the longest compatible transcript to be returned
Current behavior
The MANE select transcript, NM_001395849.1, is returned
Describe the bug
When running
genomic_to_tx_segment, it will select the MANE Select transcript, and return this transcript even if it does not exist in UTA. We should add a check to see if the transcript exists in UTA before returning.Steps to reproduce
An example where this occurs is selecting the MANE select transcript for the gene "NPIPB5"
Expected behavior
I expected the longest compatible transcript to be returned
Current behavior
The MANE select transcript, NM_001395849.1, is returned
Possible reason(s)
No response
Suggested fix
No response
Branch, commit, and/or version
0.8.0
Screenshots
No response
Environment details
Mac M1
Additional details
No response
Contribution
Yes, I can create a PR for this fix.