Please see attached taxonomy_metadata.tsv (tab-delimited) and sequences.fasta files.
These are output files from the tool QIIME2 (for further details, please see https://docs.qiime2.org).
What would be an easy/short script to extract sequences that correspond, for example, to taxon
k__Fungi;p__unidentified;c__unidentified;o__unidentified;f__unidentified;g__unidentified;s__unidentified? then run the sequences through BLAST to get the results. Please attach your results and source code
I tried using python and pandas packages but... I don't really have knowledge in programming and I'm quite lost bc of the file types and how they are structured. I don't think working with them as tables will give an easy answer.
Thanks for your help!
That's a really bad title. Please read this: https://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1002202
Using a proper title is critical to people taking the time to even read your post.