# HG changeset patch
# User jjohnson
# Date 1614171763 0
# Node ID cb58d1961fd3cd2f7959f9bae107b59cc7e3e3c5
# Parent 286e35cf3f1c43139211ee5d0b9c4e1401b3a5a0
"planemo upload commit 692ea558cbbefee93859dc2b005fab5ac4970eb8"
diff -r 286e35cf3f1c -r cb58d1961fd3 fgbio_sort_bam.xml
--- a/fgbio_sort_bam.xml Sun Feb 21 23:41:16 2021 +0000
+++ b/fgbio_sort_bam.xml Wed Feb 24 13:02:43 2021 +0000
@@ -6,17 +6,23 @@
fgbio --version
-
+
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+ pairs much like the CIGAR string in BAM files. Four kinds of operators are recognized:
- T identifies a template read
- B identifies a sample barcode read
- M identifies a unique molecular index read
- S identifies a set of bases that should be skipped or ignored
+ - T identifies a template read
+ - B identifies a sample barcode read
+ - M identifies a unique molecular index read
+ - S identifies a set of bases that should be skipped or ignored
The last pair may be specified using a + sign instead of number to denote “all remaining bases”. This is useful if, e.g., fastqs have been trimmed and contain reads of varying length. For example to convert a paired-end run with an index read and where the first 5 bases of R1 are a UMI and the second five bases are monotemplate you might specify:
diff -r 286e35cf3f1c -r cb58d1961fd3 macros.xml.bak
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macros.xml.bak Wed Feb 24 13:02:43 2021 +0000
@@ -0,0 +1,56 @@
+
+ 1.3.0
+ 0
+
+
+ fgbio
+
+
+
+ (([1-9][0-9]*[TBMS])*([+]|[1-9][0-9]*)[TBMS])
+ @READ_STRUCTURE_PATTERN@(\s@READ_STRUCTURE_PATTERN@)*
+
+ ^@READ_STRUCTURES_PATTERN@$
+
+
+ ^[A-Za-z][A-Za-z]$
+
+
+
+
+
+
+
+
+
+
+ pairs much like the CIGAR string in BAM files. Four kinds of operators are recognized:
+
+ - T identifies a template read
+ - B identifies a sample barcode read
+ - M identifies a unique molecular index read
+ - S identifies a set of bases that should be skipped or ignored
+
+The last pair may be specified using a + sign instead of number to denote “all remaining bases”. This is useful if, e.g., fastqs have been trimmed and contain reads of varying length. For example to convert a paired-end run with an index read and where the first 5 bases of R1 are a UMI and the second five bases are monotemplate you might specify:
+
+::
+
+ --input r1.fq r2.fq i1.fq --read-structures 5M5S+T +T +B
+
+Alternative if you know your reads are of fixed length you could specify:
+
+::
+
+ --input r1.fq r2.fq i1.fq --read-structures 5M5S65T 75T 8B
+
+
+]]>
+
+
+
+
+
+