#1 combine coverage with confidence

[Anonymous]

If I understand right, taxa_coverage.py predicts whether each sequence will be amplified with the given primers, and gives out a reading organized by taxa of the fraction of each category that will show amplification. And taxa_assignment_report.py looks at reads from the region you have chosen and tells you whether it will give good resolving power by indicating how many sequences at each taxonomic level have an assignment confidence >=0.8. What would be nice is to have a breakdown of these numbers by individual domains, phyla, class, etc. so we know what species are included in that number. I see that is present in the output file *_reads_assignments.txt, but is there a simple way to put it in graphical format (percentages) as is done in taxa_coverage.py?