Skip to content

Commit 484d358

Browse files
committed
some stupid GCC BAM doesn't have readgroup lines
1 parent a53f524 commit 484d358

File tree

1 file changed

+4
-4
lines changed

1 file changed

+4
-4
lines changed

wdl/pipelines/TechAgnostic/Utility/SplitBamByReadgroup.wdl

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -69,10 +69,10 @@ workflow SplitBamByReadgroup {
6969

7070
# basecall_model only applies to ONT, so PacBio data will always get 'None'
7171
Array[String] readgroup_attrs_to_get = ['ID', 'LB', 'PU']
72-
call BU.GetReadGroupInfo { input: bam = bam, keys = readgroup_attrs_to_get, null_value_representation = 'None' }
73-
String rgid = GetReadGroupInfo.read_group_info['ID']
74-
String library = GetReadGroupInfo.read_group_info['LB']
75-
String platform_unit = GetReadGroupInfo.read_group_info['PU']
72+
# call BU.GetReadGroupInfo { input: bam = bam, keys = readgroup_attrs_to_get, null_value_representation = 'None' }
73+
String rgid = "null" # GetReadGroupInfo.read_group_info['ID']
74+
String library = "null" # GetReadGroupInfo.read_group_info['LB']
75+
String platform_unit = "null" # GetReadGroupInfo.read_group_info['PU']
7676
7777
if (debug_mode) {
7878
call Utils.CountBamRecords { input: bam = bam }

0 commit comments

Comments
 (0)