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

Is there a strong repair mode? #22

Open
Trandamere opened this issue Jul 19, 2020 · 1 comment
Open

Is there a strong repair mode? #22

Trandamere opened this issue Jul 19, 2020 · 1 comment

Comments

@Trandamere
Copy link

Hi:
The software you have developed is excellent.
I would like to know if there are any strong repair options including conventional SJ in the third generation including GTAG,GCAG.
Or remove any code in the software that can be completely fixed without considering the GTAG in PacBio

@Trandamere
Copy link
Author

I disabled these lines of code to achieve the goal

def cleanNoncanonical(transcript, refs, maxDist, logInfo):
    """ Iterate over each junction in the provided transcript object. If the 
        junction is noncanonical, attempt to correct it, and record the result.
        Returns a modified transcript object that originates from a copy of the 
        object passed into the function, along with the error log entries. 
    """ 
    logInfo.corrected_NC_SJs = 0
    logInfo.uncorrected_NC_SJs = 0 
    TE_entries = ""

    #if transcript.isCanonical == True:
    #    return transcript, TE_entries

    # Iterate over all junctions and attempt to correct
    for splice_jn_num in range(0,len(transcript.spliceJunctions)):
        junction = transcript.spliceJunctions[splice_jn_num]
        #if junction.isCanonical:
        #    continue
        #else:

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant