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

MSG: trunc start,end -- there was no end for 1 #12

Closed
tseemann opened this issue Apr 3, 2018 · 1 comment
Closed

MSG: trunc start,end -- there was no end for 1 #12

tseemann opened this issue Apr 3, 2018 · 1 comment
Assignees
Labels

Comments

@tseemann
Copy link
Owner

tseemann commented Apr 3, 2018

*** [5] tig00003742 ***
Using first 3959 bp to BLAST
Writing tig00003742 ( 3959 bp ) to 5.fa
Writing tig00003742 ( 3959 bp ) to 5.head.fa
Running: blastn -query 5.head.fa -subject 5.fa -out 5.bls -evalue 1E-6 -dust no
blastn: 2..3959/3959 aligns to 2..3959/3959 at 98.5 %id
tig00003742 keep 1..0/3959 (remove 3960 bp)

------------- EXCEPTION -------------
MSG: trunc start,end -- there was no end for 1
STACK Bio::PrimarySeqI::trunc /home/linuxbrew/.linuxbrew/Cellar/perl/5.26.1_1/lib/perl5/site_perl/5.26.1/Bio/PrimarySeqI.pm:447
STACK main::check_overhang /home/tseemann/git/berokka/bin/berokka:149
STACK toplevel /home/tseemann/git/berokka/bin/berokka:74
-------------------------------------

@tseemann tseemann added the bug label Apr 3, 2018
@tseemann tseemann self-assigned this Apr 3, 2018
@tseemann
Copy link
Owner Author

tseemann commented Apr 3, 2018

Turns out this is a pacbio control sequence ~ 4000 bp!
And it is a palindrome.

I now realise that to circularise the hit must be on strand + / +

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

No branches or pull requests

1 participant