Hi Scott,

The tool "Metagenomic analyses -> Find diagnostic hits" can be used to isolate the conserved sequences. Then, you use the tool "Join, Subtract and Group -> Compare" to find "Non Matching rows of 1st dataset" to filter out anything that you think is spurious for your analysis (put in original file first, output of diagnostic hits second) before moving forward with the other summary tools.


You will probably want to run the "Find diagnostic hits" tool more than once. The choice is yours whether to do the "Compare" after each, or to "Text Manipulation -> Concatenate" all the results together first, then "Compare". The first might work faster, it just depends on the size of your datasets (how much filtering occurred before this step, etc).

The "Compare" tool sorts and holds data in memory. Even if you need to break the data up and run in smaller chunks, the results should be the same in the end. None of these jobs require the data to be in one lump.

Others are welcome to add to this with their own strategies, I am sure there are others ways to do this. Some of the public servers specializing in Metagenomics may also have tools for this, or options, and some of those may have donated to the Tool Shed, for local or cloud use. May be worth a look.
http://wiki.galaxyproject.org/PublicGalaxyServers

Good question!

Jen
Galaxy team



On 9/18/13 7:03 AM, Scott W. Tighe wrote:
Dear Galaxy

When running HiSeq shot metagenomics sample from the environment against megablast and taxonomic representation, How do I filter/remove all the 16s and other conserved sequences.

The problem if blasting a single organism that has a fraction of conserved sequence, the results will align with E.coli 10,000 times more then the possible target organism. This data would be wrong and misleading. For example a 100mg sample that was negative for e coli using MUG test, give thousands of hits with galaxy.

1) Is there a "filter conserved sequences" setting?



2) Is there a "remove model organisms" setting?


Scott Tighe

--
Jennifer Hillman-Jackson
http://galaxyproject.org

___________________________________________________________
The Galaxy User list should be used for the discussion of
Galaxy analysis and other features on the public server
at usegalaxy.org.  Please keep all replies on the list by
using "reply all" in your mail client.  For discussion of
local Galaxy instances and the Galaxy source code, please
use the Galaxy Development list:

 http://lists.bx.psu.edu/listinfo/galaxy-dev

To manage your subscriptions to this and other Galaxy lists,
please use the interface at:

 http://lists.bx.psu.edu/

To search Galaxy mailing lists use the unified search at:

 http://galaxyproject.org/search/mailinglists/

Reply via email to