Trees | Indices | Help |
|
---|
|
Script for filtering alignments with deficient properties (rqs,gaps, mismatches). The script generates a new .SAM file from the input SAM file, where only alignments are contained that fulfill certain quality thresholds. nm = numberOfMisMatches rq = readqualityscore MQ = mappingquality
|
|||
|
|||
|
|||
|
|||
|
|||
|
|
|||
counter1 = 0
|
|||
counter2 = 0
|
|||
__package__ = None
|
|
Reads a .fastqfile and filters them for certain quality values. If an alignment is equally good or better as the threshold it will be written to the filtered .SAM file. input: fastq file output: fastq dictionary key = readid; value = qualstr
|
Calculates the number of gaps / mistmatches from an HTSeq alngt. object.
|
Calculates the number of gaps / mistmatches from an HTSeq alngt. object.
|
Function for filtering a given sam file. The alignments are evaluated individually in regard to RQS,mismatches and gaps. The input needs to be sorted. |
Trees | Indices | Help |
|
---|
Generated by Epydoc 3.0.1 on Fri Mar 15 09:41:18 2013 | http://epydoc.sourceforge.net |