search

mgalardini / pyseer

SEER, reimplemented in python 🐍🔮
http://pyseer.readthedocs.io
Apache License 2.0
109 stars 27 forks source link

gene annotation #170

Closed annarerra closed 3 years ago

annarerra commented 3 years ago

Hello,

After using the --kmer mode of pyseer, i have some gene annotations like "penA|murE". I know that these features in a specific chr of N.meningitidis are close. There is any other reason for the "|" in the annotation?

Thank you in advance :)

johnlees commented 3 years ago

I think this just happens when multiple features are overlapping the mapped region (so when the unitig is long/genes are close): https://github.com/mgalardini/pyseer/blob/master/pyseer/kmer_mapping/annotate_hits.py#L69

annarerra commented 3 years ago

Great! Thank you for the answer :)