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

about the -O1 output #38

Open
ljl000 opened this issue Jan 19, 2021 · 1 comment
Open

about the -O1 output #38

ljl000 opened this issue Jan 19, 2021 · 1 comment

Comments

@ljl000
Copy link

ljl000 commented Jan 19, 2021

Hi, I used spaln to align protein sequence to genomic dna sequence and I got the result like this
<chr03:18371009-18394190 [1:23181] ( 23181 - 1 ) - >1.faa [1:70] ( 1 - 70 )
;C complement(join(22030..22242))
1.faa was my query and chr03:18371009-18394190 was the target sequence coordinate from gff3.
How can I explain the result above?
Did 23181 - 1 mean that my query align with the reverse target sequence?
what did C complement(join(22030..22242)) mean?

@ogotoh
Copy link
Owner

ogotoh commented Jan 20, 2021

Yes, your query appears to match the reverse complementary strand of the genomic sequence.

;C complement(join(22030..22242)) indicates the exonic range of the genomic sequence that matches the query. In this example, no intron exists, so that only one range (pair of coordinates that flank the exon) is shown. The format after ;C is the same as that of feature table of Genbank entries.

Osamu,

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

2 participants