GitLab

obiannotate -C option seems to work partially when the header contain python dictionaries.

Example of input:

HISEQ:204:C8E5RANXX:7:1308:3148:82868_CONS_SUB_SUB ali_length=31; seq_ab_match=31; species_name=None; family=46569; class_name=Insecta; phylum_name=Arthropoda; seq_a_deletion=0; rank=subfamily; cluster=HISEQ:204:C8E5RANXX:7:1308:3148:82868_CONS_SUB_SUB; best_identity={'order_filtered_embl_r136_noenv_COL': 0.9030303030303031}; phylum=6656; forward_match=acgctgttatcccttagg; forward_primer=acgctgttatccctwagg; reverse_primer=gacgataagaccctwtaga; species=None; merged_sample={'H20_Cs_r4': 1, 'H20_Cs_r3': 11}; forward_score=72.0; seq_a_mismatch=0; start=taatt; forward_tag=atatagcg; seq_b_mismatch=0; scientific_name=Apicotermitinae; experiment=litiere_colle; species_list={'order_filtered_embl_r136_noenv_COL': ['Anoplotermes group sp. 1 TB-2017', 'Anoplotermes group nr. E1 TB-2014', 'Anoplotermes schwarzi', 'Ruptitermes nr. xanthochiton FG-ND2-26', 'Patawatermes sp. A TB-2017', 'Patawatermes nigripunctatus', 'Aparatermes sp. A TB-2017', 'Patawatermes turricola', 'Anoplotermes banksi', 'Humutermes krishnai', 'Grigiotermes hageni', 'Longustitermes manni']}; seq_a_single=94; cluster_weight=12; reverse_score=76.0; count=12; seq_b_insertion=0; taxid_by_db={'order_filtered_embl_r136_noenv_COL': 92739}; id_status={'order_filtered_embl_r136_noenv_COL': True}; seq_b_deletion=0; genus_name=None; cluster_center=True; seq_a_insertion=0; seq_length_ori=219; reverse_tag=agactatg; class=50557; goodAli=Alignement; match_count={'order_filtered_embl_r136_noenv_COL': 13}; family_name=Termitidae; best_match={'order_filtered_embl_r136_noenv_COL': 'KY224594; order_name=Blattodea; taxid=92739; scientific_name_by_db={'order_filtered_embl_r136_noenv_COL': 'Apicotermitinae'}; cluster_score=1.0; genus=None; order=85823; '}; order_name=Blattodea; rank_by_db={'order_filtered_embl_r136_noenv_COL': 'subfamily'}; seq_length=162; seq_b_single=94; status=full; mode=alignment; position=08_01D; genus=None; order=85823; distance=0.0; taatttaatcttataatcaaaataaatggatcaaaaaactataaacaaatatatagcagt aaagaggagttaaataaattcctcccatcgccccaacaaaacacctaaatcacttaataa aacaaaacaaacaaaataataaaaagtaaataaaatgttaac

command: obiannotate -C toto.fasta

Example corresponding output:

HISEQ:204:C8E5RANXX:7:1308:3148:82868_CONS_SUB_SUB '}; order_name=Blattodea; rank_by_db={'order_filtered_embl_r136_noenv_COL': 'subfamily'}; seq_length=162; seq_b_single=94; status=full; mode=alignment; position=08_01D; genus=None; order=85823; distance=0.0; taatttaatcttataatcaaaataaatggatcaaaaaactataaacaaatatatagcagt aaagaggagttaaataaattcctcccatcgccccaacaaaacacctaaatcacttaataa aacaaaacaaacaaaataataaaaagtaaataaaatgttaac