Closed IreneGregori closed 2 years ago
Hi Irene,
You should be able to do it with obifiles_to_metabarlist
(https://metabarfactory.github.io/metabaR/reference/obifiles_to_metabarlist.html) using those files:
file_obitab
: path to the file produced by an export to tab format of the ecotag result (the command should look like obi export --tab dms/ecotag_result > output.tab
)file_ngsfilter
: path to the ngsfilter file (the same file that you imported in the dms if you used ngsfilter)file_samples
: path to a 'sample characteristics table'. This does not exist in the OBITools3, it should look like this file here, but I believe that it can simply be the list of sample names with one sample per line.I'm supposed to make a function to import OBI3 data seamlessly but haven't done it yet..... might look into it these days if I can.
Please tell me if that works and ask for more information if needed :) Cheers Celine
Hi Celine, thanks for the answer. Unfortunately, the command didn't work. The output in R is "_Error in file.exists(fileobitab) : invalid 'file' argument" and i can't fix this. Could you please send me an example of how the file_obitab should look like? I saw the one in "litiere" files and it looks a lot different from mine. Thanks
Hi Irene, I looked into it, I think I can get something working for next Tuesday. It will be an obi export
option to produce the files for obifiles_to_metabarlist
. If you send me (celine.mercier@metabarcoding.org) your ngsfilter file, and your sample characteristics file if you have one, I can test it on them, but I can do without. Cheers
Hi Celine, i send you my ngsfile and my sample file. Thank you very much!
Ciao, Irene
Il giorno ven 15 ott 2021 alle ore 05:24 celine-mercier < @.***> ha scritto:
Hi Irene, I looked into it, I think I can get something working for next Tuesday. It will be an obi export option to produce the files for obifiles_to_metabarlist. If you send me @.***) your ngsfilter file, and your sample characteristics file if you have one, I can test it on them, but I can do without. Cheers
— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/metabaRfactory/metabaR/issues/39#issuecomment-943955628, or unsubscribe https://github.com/notifications/unsubscribe-auth/AWAJKHUL3I5ECUXBGPCGXETUG6NHNANCNFSM5FYFRNPA . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.
-- Irene Gregori
1 01 acacacac:gactgatg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 02 acagcaca:gactgatg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 03 gtgtacat:gactgatg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 04 tatgtcag:gactgatg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 05 tagtcgca:gactgatg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 06 tactatac:gactgatg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 07 agactatg:gactgatg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 08 gatcgcga:gactgatg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 09 acacacac:gtcacgtc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 10 acagcaca:gtcacgtc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 11 gtgtacat:gtcacgtc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 12 tatgtcag:gtcacgtc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 13 tagtcgca:gtcacgtc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 14 tactatac:gtcacgtc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 15 agactatg:gtcacgtc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 16 gatcgcga:gtcacgtc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 17 acacacac:gtcgtaga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 18 acagcaca:gtcgtaga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 19 gtgtacat:gtcgtaga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 20 tatgtcag:gtcgtaga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 21 tagtcgca:gtcgtaga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 22 tactatac:gtcgtaga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 23 agactatg:gtcgtaga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 24 gatcgcga:gtcgtaga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 25 acacacac:cgctctcg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 26 acagcaca:cgctctcg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 27 gtgtacat:cgctctcg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 28 tatgtcag:cgctctcg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 29 tagtcgca:cgctctcg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 30 tactatac:cgctctcg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 31 agactatg:cgctctcg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 32 gatcgcga:cgctctcg TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 33 acacacac:gatcgcga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 34 acagcaca:gatcgcga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 35 gtgtacat:gatcgcga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 36 tatgtcag:gatcgcga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 37 tagtcgca:gatcgcga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 38 tactatac:gatcgcga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 39 agactatg:gatcgcga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 40 gatcgcga:gatcgcga TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 41 acacacac:actagatc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 42 acagcaca:actagatc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 43 gtgtacat:actagatc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 44 tatgtcag:actagatc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 45 tagtcgca:actagatc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 46 tactatac:actagatc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 47 agactatg:actagatc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 48 gatcgcga:actagatc TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 49 acacacac:tactatac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 50 acagcaca:tactatac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 51 gtgtacat:tactatac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 52 tatgtcag:tactatac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 53 tagtcgca:tactatac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 54 tactatac:tactatac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 55 agactatg:tactatac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 56 gatcgcga:tactatac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 57 acacacac:tagtcgca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 58 acagcaca:tagtcgca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 59 gtgtacat:tagtcgca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 60 tatgtcag:tagtcgca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 61 tagtcgca:tagtcgca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 62 tactatac:tagtcgca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 63 agactatg:tagtcgca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 64 gatcgcga:tagtcgca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 65 acacacac:tatgtcag TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 66 acagcaca:tatgtcag TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 67 gtgtacat:tatgtcag TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 68 tatgtcag:tatgtcag TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 69 tagtcgca:tatgtcag TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 70 tactatac:tatgtcag TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 71 agactatg:tatgtcag TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 72 gatcgcga:tatgtcag TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 73 acacacac:gtgtacat TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 74 acagcaca:gtgtacat TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 75 gtgtacat:gtgtacat TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 76 tatgtcag:gtgtacat TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 77 tagtcgca:gtgtacat TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 78 tactatac:gtgtacat TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 79 agactatg:gtgtacat TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 80 gatcgcga:gtgtacat TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 81 acacacac:acagcaca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 82 acagcaca:acagcaca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 83 gtgtacat:acagcaca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 84 tatgtcag:acagcaca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 85 tagtcgca:acagcaca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 86 tactatac:acagcaca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 87 agactatg:acagcaca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 88 gatcgcga:acagcaca TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 89 acacacac:acacacac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 90 acagcaca:acacacac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 91 gtgtacat:acacacac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 92 tatgtcag:acacacac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 93 tagtcgca:acacacac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 94 tactatac:acacacac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 95 agactatg:acacacac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG 1 96 gatcgcga:acacacac TCCAACATCGAGGTCGCAA ACGAGAAGACCCTRTYGAG
Hi Irene,
Could you try sending me the files attached in a new email to celine.mercier.bioinfo@gmail.com, and not as a reply to this? I can't really see them like this.
Thanks! Cheers
In order to create a metabarlist file, I must create 4 different files (motus, pcrs,..). How can I create those files starting from an Obitools3 output? In particullary, which files do I have to take from all the outputs made by Obitools3? Thank you