I am now using samtools (1.3.1-33) mpileup to create mpileup files, and would like to remove reads that have flags of 4, 8, 256, 512 and 1024. So should the --rf be set to 1804 (4+8+256+512+1024) or the --ff set to 1804? Also, what is the difference between options of -a and -a -a (-aa)? Thanks you!
•
written
by
jing.mengrabbit •
If you want to remove reads that have certain bits set, then you need to use --ff, --excl-flags. If you don't set anything, the default is:
--ff [UNMAP,SECONDARY,QCFAIL,DUP]
That is, don't look at unaligned, secondary, QC-failed, or duplicate reads.
You should look at the flag table to see which you require:
upload image
Kevin
Thanks Kevin for your answer. Do you know what is the difference between -a and -a -a (-aa)?
By using either of these options, you should get a very large file.
-a, I believe, will output a record for every base that's listed in your reference sequence and for which there has been some form of
alignment.-aa or -a -a will get you an even larger file because it will
output a record for every base in your reference, irrespective of
whether there has been alignment or not.You should check how they both behave, just to make sure.
I think that these options are mainly fr use in conjunction with a target BED file
Thanks for your time. The following is the reply from samtools developers about the difference between -a and -a -a (-aa):
-a and -aa is all vs absolutely all (bad term I know). The difference is to do with references with zero coverage.
For example if we have a file with @SQ headers for all chromosomes, but only data on chr 1, then -a would show all of chr1 irrespective of coverage, but none on chr2, 3, etc, while -aa would show the zero coverage on all the other chromosomes too.
Thanks for your time. The following is the reply from samtools developers about the difference between -a and -a -a (-aa):
-a and -aa is all vs absolutely all (bad term I know). The difference is to do with references with zero coverage.
For example if we have a file with @SQ headers for all chromosomes, but only data on chr 1, then -a would show all of chr1 irrespective of coverage, but none on chr2, 3, etc, while -aa would show the zero coverage on all the other chromosomes too.
Log in to answer this question.
