| Index: generated/googleapis_beta/lib/genomics/v1beta2.dart
|
| diff --git a/generated/googleapis_beta/lib/genomics/v1beta2.dart b/generated/googleapis_beta/lib/genomics/v1beta2.dart
|
| index b0e678f313ba665d29754bc8a3a0a5f2cfec6461..7386460c551dbeb11d0dbfacc21ea9be7de013c8 100644
|
| --- a/generated/googleapis_beta/lib/genomics/v1beta2.dart
|
| +++ b/generated/googleapis_beta/lib/genomics/v1beta2.dart
|
| @@ -3940,7 +3940,25 @@ class Range {
|
| * reference sequence, in addition to metadata about the fragment (the molecule
|
| * of DNA sequenced) and the read (the bases which were read by the sequencer).
|
| * A read is equivalent to a line in a SAM file. A read belongs to exactly one
|
| - * read group and exactly one read group set.
|
| + * read group and exactly one read group set. Generating a reference-aligned
|
| + * sequence string When interacting with mapped reads, it's often useful to
|
| + * produce a string representing the local alignment of the read to reference.
|
| + * The following pseudocode demonstrates one way of doing this:
|
| + * out = "" offset = 0 for c in read.alignment.cigar { switch c.operation { case
|
| + * "ALIGNMENT_MATCH", "SEQUENCE_MATCH", "SEQUENCE_MISMATCH": out +=
|
| + * read.alignedSequence[offset:offset+c.operationLength] offset +=
|
| + * c.operationLength break case "CLIP_SOFT", "INSERT": offset +=
|
| + * c.operationLength break case "PAD": out += repeat("*", c.operationLength)
|
| + * break case "DELETE": out += repeat("-", c.operationLength) break case "SKIP":
|
| + * out += repeat(" ", c.operationLength) break case "CLIP_HARD": break } }
|
| + * return out
|
| + * Converting to SAM's CIGAR string The following pseudocode generates a SAM
|
| + * CIGAR string from the cigar field. Note that this is a lossy conversion
|
| + * (cigar.referenceSequence is lost).
|
| + * cigarMap = { "ALIGNMENT_MATCH": "M", "INSERT": "I", "DELETE": "D", "SKIP":
|
| + * "N", "CLIP_SOFT": "S", "CLIP_HARD": "H", "PAD": "P", "SEQUENCE_MATCH": "=",
|
| + * "SEQUENCE_MISMATCH": "X", } cigarStr = "" for c in read.alignment.cigar {
|
| + * cigarStr += c.operationLength + cigarMap[c.operation] } return cigarStr
|
| */
|
| class Read {
|
| /**
|
| @@ -3954,12 +3972,12 @@ class Read {
|
| core.List<core.int> alignedQuality;
|
|
|
| /**
|
| - * The bases of the read sequence contained in this alignment record.
|
| - * alignedSequence and alignedQuality may be shorter than the full read
|
| - * sequence and quality. This will occur if the alignment is part of a
|
| - * chimeric alignment, or if the read was trimmed. When this occurs, the CIGAR
|
| - * for this read will begin/end with a hard clip operator that will indicate
|
| - * the length of the excised sequence.
|
| + * The bases of the read sequence contained in this alignment record, without
|
| + * CIGAR operations applied. alignedSequence and alignedQuality may be shorter
|
| + * than the full read sequence and quality. This will occur if the alignment
|
| + * is part of a chimeric alignment, or if the read was trimmed. When this
|
| + * occurs, the CIGAR for this read will begin/end with a hard clip operator
|
| + * that will indicate the length of the excised sequence.
|
| */
|
| core.String alignedSequence;
|
|
|
|
|