Closed PeteHaitch closed 10 years ago
Reasons that reads can fail QC:
XR
or XG
-tag (this is checked multiple times, can probably remove some of these checks)0x01
, 0x40
or 0x80
FLAG
bit (paired-read FLAG
information)XG
or XR
-tags for are inconsistent with OT-strand or OB-strand (checked multiple times)comethylation_state
--overlappingPairedEndFilter
This will be simplified to the following categories:
comethylation_state
--overlappingPairedEndFilter
XR
or XG
tag Reads that fail QC will be written to the file <sampleName>.reads_that_failed_QC.txt
along with a message describing why the read failed. It might be a good idea to gzip this file to save space.
Additionally, a tally of the number of reads that failed each QC step will be reported to STDOUT
.
Rather than just output these to
STDOUT
, create a file of theQNAME
and the reason a read was filtered out during QC steps. Add the option to not create this file which instead would direct this output to/dev/null