Hi Yan,

Yes, "Slice BAM" is like a filter, so the output should be smaller if some of the data mapped to regions outside of the intervals (entirely - no overlap). You can run flagstat on both files and compare the two numbers to see the difference/what was filtered.

For QC-passed/failed - this is a great question without a simple answer (as it has some history). Or so I thought. Then, I found a very simple post by Prof Albert on Biostar addressing the same question, that I think is about as good as it can get! (Main point = can safely ignore)
http://www.biostars.org/p/16100/

Hope this helps!

Jen
Galaxy team

On 9/12/13 10:11 PM, Yan He wrote:

Hi Jen and other galaxy-users,

 

I was using “Slice BAM” tool on Galaxy to get the alignment overlap with the targeted intervals. After I got the output BAM file, I used “flagstat” to get the detailed information of the output BAM file. What I got from “flagstat” is as following.

 

13704486 + 0 in total (QC-passed reads + QC-failed reads)

0 + 0 duplicates

2989995 + 0 mapped (21.82%:-nan%)

13704486 + 0 paired in sequencing

 

What’s the QC-passed reads? What’s the mapped reads? Should I only get the mapped reads to the targeted intervals? I am very confused. Any help is highly appreciated! Thanks a lot!

 

Best wishes,

Yan

 

 

 

 

 

 



___________________________________________________________
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/

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