Skip to content

Commit 3e30a53

Browse files
committed
Restore portable grep in CRAM tests
1 parent e2b3183 commit 3e30a53

File tree

5 files changed

+46
-46
lines changed

5 files changed

+46
-46
lines changed

tests/cram/splitsample.t

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -7,22 +7,22 @@
77
$ ls -l $CRAMTMP/split.*.bam | wc -l | tr -d ' '
88
3
99

10-
$ ID=$(samtools view -F 4 -H $CRAMTMP/split.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep 3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494 | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
11-
$ samtools view -F 4 $CRAMTMP/split.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
10+
$ ID=$(samtools view -F 4 -H $CRAMTMP/split.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep 3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494 | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
11+
$ samtools view -F 4 $CRAMTMP/split.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
1212
0
13-
$ samtools view -F 4 $CRAMTMP/split.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
13+
$ samtools view -F 4 $CRAMTMP/split.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
1414
9
1515

16-
$ ID=$(samtools view -F 4 -H $CRAMTMP/split.test_test.bam | grep test_test | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
17-
$ samtools view -F 4 $CRAMTMP/split.test_test.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
16+
$ ID=$(samtools view -F 4 -H $CRAMTMP/split.test_test.bam | grep test_test | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
17+
$ samtools view -F 4 $CRAMTMP/split.test_test.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
1818
0
19-
$ samtools view -F 4 $CRAMTMP/split.test_test.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
19+
$ samtools view -F 4 $CRAMTMP/split.test_test.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
2020
11
2121

22-
$ ID=$(samtools view -F 4 -H $CRAMTMP/split.UCLA_1023.bam | grep UCLA_1023 | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
23-
$ samtools view -F 4 $CRAMTMP/split.UCLA_1023.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
22+
$ ID=$(samtools view -F 4 -H $CRAMTMP/split.UCLA_1023.bam | grep UCLA_1023 | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
23+
$ samtools view -F 4 $CRAMTMP/split.UCLA_1023.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
2424
0
25-
$ samtools view -F 4 $CRAMTMP/split.UCLA_1023.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
25+
$ samtools view -F 4 $CRAMTMP/split.UCLA_1023.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
2626
10
2727

2828
$ $__PBTEST_PBMM2_EXE align $MERGED $REF $CRAMTMP/split_dataset.alignmentset.xml --split-by-sample
@@ -56,7 +56,7 @@ When both --split-by-sample and --sample were set, expect to see only one bam fi
5656
$ $__PBTEST_PBMM2_EXE align $REF $IN $CRAMTMP/splitsampleoverride.consensusalignmentset.xml --sort -j 8 --split-by-sample --sample "MySample" 2>&1 | fgrep -v 'Requested more threads'
5757
*Options --split-by-sample and --sample are mutually exclusive. Option --sample will be applied and --split-by-sample is ignored! (glob)
5858
$ [[ -f $CRAMTMP/splitsampleoverride.bam ]] || echo "File does not exist!"
59-
$ samtools view -H $CRAMTMP/splitsampleoverride.bam | grep "@RG" | grep -vP "@PG\tID:samtools" | cut -f 6 | sort | uniq
59+
$ samtools view -H $CRAMTMP/splitsampleoverride.bam | grep "@RG" | grep -v "@PG ID:samtools" | cut -f 6 | sort | uniq
6060
SM:MySample
6161

6262
$ $__PBTEST_PBMM2_EXE align --split-by-sample $NO_SM_BIOSAMPLES $REF $CRAMTMP/split-no-sm.bam

tests/cram/splitsamplejson.t

Lines changed: 9 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -28,22 +28,22 @@
2828
$ grep -c split_dataset_sorted_json.UCLA_1023.alignmentset.xml $CRAMTMP/split_dataset_sorted_json.json
2929
1
3030

31-
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted_json.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep 3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494 | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
32-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
31+
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted_json.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep 3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494 | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
32+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
3333
0
34-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
34+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
3535
9
3636

37-
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted_json.test_test.bam | grep test_test | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
38-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.test_test.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
37+
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted_json.test_test.bam | grep test_test | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
38+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.test_test.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
3939
0
40-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.test_test.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
40+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.test_test.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
4141
11
4242

43-
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted_json.UCLA_1023.bam | grep UCLA_1023 | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
44-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.UCLA_1023.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
43+
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted_json.UCLA_1023.bam | grep UCLA_1023 | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
44+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.UCLA_1023.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
4545
0
46-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.UCLA_1023.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
46+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.UCLA_1023.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
4747
10
4848

4949
$ grep -c split_dataset_sorted_json.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.alignmentset.xml $CRAMTMP/split_dataset_sorted_json.json | tr -d ' '

tests/cram/splitsamplejsonsorted.t

Lines changed: 9 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -28,20 +28,20 @@
2828
$ grep -c split_dataset_sorted_json.UCLA_1023.alignmentset.xml $CRAMTMP/split_dataset_sorted_json.json
2929
1
3030

31-
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted_json.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep 3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494 | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
32-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
31+
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted_json.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep 3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494 | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
32+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
3333
0
34-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
34+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
3535
9
3636

37-
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted_json.test_test.bam | grep test_test | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
38-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.test_test.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
37+
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted_json.test_test.bam | grep test_test | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
38+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.test_test.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
3939
0
40-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.test_test.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
40+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.test_test.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
4141
11
4242

43-
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted_json.UCLA_1023.bam | grep UCLA_1023 | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
44-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.UCLA_1023.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
43+
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted_json.UCLA_1023.bam | grep UCLA_1023 | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
44+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.UCLA_1023.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
4545
0
46-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.UCLA_1023.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
46+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted_json.UCLA_1023.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
4747
10

tests/cram/splitsamplesamename.t

Lines changed: 9 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -6,22 +6,22 @@
66
$ ls -l $CRAMTMP/splitname.*.bam | wc -l | tr -d ' '
77
3
88

9-
$ ID=$(samtools view -F 4 -H $CRAMTMP/splitname.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep 3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494 | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
10-
$ samtools view -F 4 $CRAMTMP/splitname.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
9+
$ ID=$(samtools view -F 4 -H $CRAMTMP/splitname.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep 3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494 | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
10+
$ samtools view -F 4 $CRAMTMP/splitname.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
1111
0
12-
$ samtools view -F 4 $CRAMTMP/splitname.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
12+
$ samtools view -F 4 $CRAMTMP/splitname.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
1313
9
1414

15-
$ ID=$(samtools view -F 4 -H $CRAMTMP/splitname.test-0.bam | grep "test(" | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
16-
$ samtools view -F 4 $CRAMTMP/splitname.test-0.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
15+
$ ID=$(samtools view -F 4 -H $CRAMTMP/splitname.test-0.bam | grep "test(" | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
16+
$ samtools view -F 4 $CRAMTMP/splitname.test-0.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
1717
0
18-
$ samtools view -F 4 $CRAMTMP/splitname.test-0.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
18+
$ samtools view -F 4 $CRAMTMP/splitname.test-0.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
1919
10
2020

21-
$ ID=$(samtools view -F 4 -H $CRAMTMP/splitname.test-1.bam | grep "test)" | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
22-
$ samtools view -F 4 $CRAMTMP/splitname.test-1.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
21+
$ ID=$(samtools view -F 4 -H $CRAMTMP/splitname.test-1.bam | grep "test)" | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
22+
$ samtools view -F 4 $CRAMTMP/splitname.test-1.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
2323
0
24-
$ samtools view -F 4 $CRAMTMP/splitname.test-1.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
24+
$ samtools view -F 4 $CRAMTMP/splitname.test-1.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
2525
11
2626

2727
$ $__PBTEST_PBMM2_EXE align $MERGED $REF $CRAMTMP/split_dataset.alignmentset.xml --split-by-sample

tests/cram/splitsamplesorted.t

Lines changed: 9 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -27,20 +27,20 @@
2727
$ grep -c split_dataset_sorted.UCLA_1023.alignmentset.xml $CRAMTMP/split_dataset_sorted.json
2828
1
2929

30-
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep 3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494 | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
31-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
30+
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep 3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494 | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
31+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
3232
0
33-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
33+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted.3260208_188nM-GTAC_2xGCratio_LP7_100fps_15min_5kEColi_SP2p1_3uMSSB_BA243494.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
3434
9
3535

36-
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted.test_test.bam | grep test_test | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
37-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted.test_test.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
36+
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted.test_test.bam | grep test_test | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
37+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted.test_test.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
3838
0
39-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted.test_test.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
39+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted.test_test.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
4040
11
4141

42-
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted.UCLA_1023.bam | grep UCLA_1023 | grep -vP "@PG\tID:samtools" | cut -f 2 | cut -f 2 -d ':')
43-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted.UCLA_1023.bam | grep -vc ${ID} | grep -vP "@PG\tID:samtools" | tr -d ' '
42+
$ ID=$(samtools view -F 4 -H $CRAMTMP/split_dataset_sorted.UCLA_1023.bam | grep UCLA_1023 | grep -v "@PG ID:samtools" | cut -f 2 | cut -f 2 -d ':')
43+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted.UCLA_1023.bam | grep -vc ${ID} | grep -v "@PG ID:samtools" | tr -d ' '
4444
0
45-
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted.UCLA_1023.bam | grep -vP "@PG\tID:samtools" | wc -l | tr -d ' '
45+
$ samtools view -F 4 $CRAMTMP/split_dataset_sorted.UCLA_1023.bam | grep -v "@PG ID:samtools" | wc -l | tr -d ' '
4646
10

0 commit comments

Comments
 (0)