# Extracting WGS reads belonging to a particular organism **URL:** \ **Category:** Uncategorized\ **Tags:** gtn-tutorial, metagenomics, mapping\ **Created:** [January 11, 2019, 11:55am UTC](https://help.galaxyproject.org/t/extracting-wgs-reads-belonging-to-a-particular-organism/432 "2019-01-11T11:55:02Z")\ **Posts on this page:** 9\ **Page:** 1
**Author:** ![manju\_moorthy](https://avatars.discourse-cdn.com/v4/letter/m/3be4f8/32.png) [@manju\_moorthy](https://help.galaxyproject.org/u/manju_moorthy)\ **Post date:** [January 11, 2019, 11:55am UTC](https://help.galaxyproject.org/t/extracting-wgs-reads-belonging-to-a-particular-organism/432/1 "2019-01-11T11:55:02Z")
I have the sequenced data of an environmental culture. It has around 3 16srRNA belonging to 3 different organism. I want to extract the reads belonging to each different organism. How can I do that? Also when using Bowtie its not getting mapped to the reference genome of the 3 organisms at all. What to do? ---
**Author:** ![jennaj](https://sea2.discourse-cdn.com/flex020/user_avatar/help.galaxyproject.org/jennaj/32/27_2.png) [@jennaj](https://help.galaxyproject.org/u/jennaj)\ **Post date:** [January 11, 2019, 5:12pm UTC](https://help.galaxyproject.org/t/extracting-wgs-reads-belonging-to-a-particular-organism/432/2 "2019-01-11T17:12:47Z")
Hello, Mapping to a reference genome with Bowtie is unlikely to capture 16s hits. For how to classify 16s reads (and other analysis operations), please see the Metagenomics tutorials from the GTN. > [@](#): > > [Galaxy Training Network Tutorials](https://galaxyproject.github.io/training-material/): Some [GTN](https://galaxyproject.org/teach/gtn/) tutorials are appropriate for [Galaxy Main](https://usegalaxy.org) and some are not. Where you can run each is noted per tutorial – click on the _Galaxy instances_ gear icon ![galaxyserversGTN](https://us1.discourse-cdn.com/flex020/uploads/galaxy/original/1X/1287f50cfcdaa1aad750fe323f96995260ac4947.png) to review the [Public Galaxy](https://galaxyproject.org/use/) server choices. If a tutorial is supported by a pre-configured [Galaxy Docker](https://hub.docker.com/r/bgruening/galaxy-stable/) training image, instructions for how to get it will be listed below the tutorial listings, per category. Thanks! ---
**Author:** ![manju\_moorthy](https://avatars.discourse-cdn.com/v4/letter/m/3be4f8/32.png) [@manju\_moorthy](https://help.galaxyproject.org/u/manju_moorthy)\ **Post date:** [January 11, 2019, 5:29pm UTC](https://help.galaxyproject.org/t/extracting-wgs-reads-belonging-to-a-particular-organism/432/3 "2019-01-11T17:29:01Z")
OK. Thank you. Let me see the metagenomics tutorials ---
**Author:** ![manju\_moorthy](https://avatars.discourse-cdn.com/v4/letter/m/3be4f8/32.png) [@manju\_moorthy](https://help.galaxyproject.org/u/manju_moorthy)\ **Post date:** [January 11, 2019, 5:30pm UTC](https://help.galaxyproject.org/t/extracting-wgs-reads-belonging-to-a-particular-organism/432/4 "2019-01-11T17:30:31Z")
So do you mean to say that , even when the 16s sequences in our genome is similar to the 16s sequences in the reference genome, bowtie wont show hits? ---
**Author:** ![jennaj](https://sea2.discourse-cdn.com/flex020/user_avatar/help.galaxyproject.org/jennaj/32/27_2.png) [@jennaj](https://help.galaxyproject.org/u/jennaj)\ **Post date:** [January 11, 2019, 6:23pm UTC](https://help.galaxyproject.org/t/extracting-wgs-reads-belonging-to-a-particular-organism/432/5 "2019-01-11T18:23:20Z")
It is better to use a dedicated 16s reference database (example: [Silva](https://www.arb-silva.de/)). Reasons why can include: - Full reference genomes can have ribosomal regions masked. - 16s reads do not map with specificity against a full reference genome. - 16s reads do not meet the minimum mapping criteria set (eg: short length, too much variation). ---
**Author:** ![manju\_moorthy](https://avatars.discourse-cdn.com/v4/letter/m/3be4f8/32.png) [@manju\_moorthy](https://help.galaxyproject.org/u/manju_moorthy)\ **Post date:** [January 11, 2019, 6:55pm UTC](https://help.galaxyproject.org/t/extracting-wgs-reads-belonging-to-a-particular-organism/432/6 "2019-01-11T18:55:12Z")
But mine is not a 16s sequencing data. I am having whole genome sequence data. ---
**Author:** ![jennaj](https://sea2.discourse-cdn.com/flex020/user_avatar/help.galaxyproject.org/jennaj/32/27_2.png) [@jennaj](https://help.galaxyproject.org/u/jennaj)\ **Post date:** [January 11, 2019, 7:59pm UTC](https://help.galaxyproject.org/t/extracting-wgs-reads-belonging-to-a-particular-organism/432/7 "2019-01-11T19:59:51Z")
> [@manju\_moorthy](#): > > It has around 3 16srRNA belonging to 3 different organism. Thanks for clarifying, I thought this ^^ was what you had sequenced. > [@manju\_moorthy](#): > > Also when using Bowtie its not getting mapped to the reference genome of the 3 organisms at all. What to do? WGS reads can be mapped with Bowtie2 or BWA/BWA-MEM or BLASTN against a reference genome. Success involves doing some QA/QC on the reads, setting the correct parameters for the mapping tool used, and (possibly, if a built-in index is not available) making sure the custom reference genome is formatted correctly. See the tool forms for abbreviated help and links to the full documentation. Parameters available on Galaxy tool forms are for the most part the same as when using the tools line command. Many of the Galaxy tutorials include QA/QC and mapping steps. Help for Custom Genome formatting is [here](https://galaxyproject.org/learn/custom-genomes/). Also please see: [Troubleshooting resources for errors or unexpected results](https://help.galaxyproject.org/t/troubleshooting-resources-for-errors-or-unexpected-results/42) ---
**Author:** ![manju\_moorthy](https://avatars.discourse-cdn.com/v4/letter/m/3be4f8/32.png) [@manju\_moorthy](https://help.galaxyproject.org/u/manju_moorthy)\ **Post date:** [January 12, 2019, 7:07am UTC](https://help.galaxyproject.org/t/extracting-wgs-reads-belonging-to-a-particular-organism/432/8 "2019-01-12T07:07:40Z")
I have tried Bowtie, Bowtie2 and localblast too. But very very low hits or no mapping found with all the 3 reference genomes ---
**Author:** ![jennaj](https://sea2.discourse-cdn.com/flex020/user_avatar/help.galaxyproject.org/jennaj/32/27_2.png) [@jennaj](https://help.galaxyproject.org/u/jennaj)\ **Post date:** [March 17, 2026, 6:58am UTC](https://help.galaxyproject.org/t/extracting-wgs-reads-belonging-to-a-particular-organism/432/9 "2026-03-17T06:58:42Z")