-
Notifications
You must be signed in to change notification settings - Fork 5
Description
For some proteins, I get the following funny offset in the output (for position 121). Not consistent. Must be a trivial bug somewhere. Volker
query
MPAWPAAAVAKAVPSPSTPPPPHSRGAGRRRLRPCGAKKGPGTDERGATAGGGGVVTRGA
LLRSGAALFALGFVDAGYSGDWSRIGAISKDTEEALKLAAYAVVPLCLAVVFSPSSEDGS
NNT*
gdna
AACCGTTTGCATCGGCGGGAGTGCCAGGTTTTGTCCGTACCTCGCACCCCGAGAGCGTGCAGCGCTCCACGTTCTGCTGCTAACGCACGTGTGGACGAGAAGCAAACACCCAAACACCCAACTACCACCGCGCACACAGCACAGCATGCCGGCGTGGCCCGCCGCCGCCGTCGCCAAGGCCGTCCCGTCCCCGTCCACGCCGCCGCCGCCGCACTCGCGCGGAGCAGTCCGCCGTCGCCTCCGCCCGTGCGGCGCAAAGAAGGGCCCCGGGACCGACGAGCGAGGTGCCACCGCCGGCGGCGGCGGCGTGGTCACCAGGGGCGCCCTCCTGCGGTCCGGCGCCGCTCTCTTCGCGCTCGGCTTCGTCGACGCCGGGTGAGTCGCACGCACGCGCGCGCGCGCGCGCGCTGATGACGGATTTCCATTGCGTTTTGCACCGGTGGTCGATCTGGTCGGTCTCACCCGTTTCGGGGCGGCTTTGCTGCAGGTACAGCGGCGACTGGTCACGCATCGGCGCCATCTCCAAGGACACCGAGGAGGCGCTCAAGCTGGCCGCCTACGCCGTCGTGCCTCTCTGCCTCGCCGTCGTGTTCTCGCCATCGTCGGAAGACGGCAGTAACAACACCTGAATCCTAGGAAGAAGAGCCATGCTGCAGATCATGTAGTCGTCGTGATACCATCATGCTTACGTGTTGTATACTTCGATCGTGATCGCAATGGGCATTTCTTTTGCCTAAGGGGTGTGTTATCAGCAGTGAAAAGTTTTGGCGTGTCACATCGAATATACGAATATAGATTGATAGCAAAACAAATTGCAGATTCCGTCTATAAATTACGAGACGAATTTATTATTTAATTAATCCATCATTAGCAAATATTTACTGTACCACCACATTATCAAATCATAGAACAATTAGGCTTAAAAGATTCATCTCGTAATTTACACACAATCTGTGTAATTAGTTATTTACATTTAATGTGACTGAGTGAAAATTTTTTG
gth -protein query -genomic gdna -species rice > out
more out
$ GenomeThreader 1.7.1
$ Date run: 2021-01-06 14:53:45
$ Arguments: -protein query -genomic gdna -species rice
Protein Sequence: file=query, description=query
1 MPAWPAAAVA KAVPSPSTPP PPHSRGAGRR RLRPCGAKKG PGTDERGATA GGGGVVTRGA
61 LLRSGAALFA LGFVDAGYSG DWSRIGAISK DTEEALKLAA YAVVPLCLAV VFSPSSEDGS
121 NNT*
Genomic Template: file=gdna, strand=+, from=1, to=926, description=gdna
....