Skip to content

Commit

Permalink
describe the filtering criteria in SAM/BAM/BAMPE parser
Browse files Browse the repository at this point in the history
  • Loading branch information
taoliu committed Oct 7, 2024
1 parent 43140cb commit 7bbcfa8
Showing 1 changed file with 16 additions and 0 deletions.
16 changes: 16 additions & 0 deletions docs/source/docs/SAMBAMBAMPE.md
Original file line number Diff line number Diff line change
Expand Up @@ -24,3 +24,19 @@ The only exception is that the `callvar` command in MACS3 will need to
read the actual sequences and mapping information such as the
mismatch, insertion, deletion, etc from BAM. Also, please make sure
that the BAM file for `callvar` has to be sorted and indexed.

## Filtering Criteria

MACS3 will apply the following rule to throw away bad alignments while
scanning the 'flag' field in SAM/BAM/BAMPE. A record will be
discarded:

1. If a record is an unmapped read, not a primary alignment, or for a
read that fails QC, or a supplementary alignment
2. If this is for paired reads:
1. this read is not paired
2. the mate of this read (the other read in the pair) is not
mapped
3. this read is the second read in the pair (NOTE: the first read
in the pair is already enough for MACS to get enough
information)

0 comments on commit 7bbcfa8

Please sign in to comment.