|
7 | 7 | $ ls -l $CRAMTMP/split.*.bam | wc -l | tr -d ' '
|
8 | 8 | 3
|
9 | 9 |
|
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 ' ' |
12 | 12 | 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 ' ' |
14 | 14 | 9
|
15 | 15 |
|
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 ' ' |
18 | 18 | 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 ' ' |
20 | 20 | 11
|
21 | 21 |
|
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 ' ' |
24 | 24 | 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 ' ' |
26 | 26 | 10
|
27 | 27 |
|
28 | 28 | $ $__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
|
56 | 56 | $ $__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'
|
57 | 57 | *Options --split-by-sample and --sample are mutually exclusive. Option --sample will be applied and --split-by-sample is ignored! (glob)
|
58 | 58 | $ [[ -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 |
60 | 60 | SM:MySample
|
61 | 61 |
|
62 | 62 | $ $__PBTEST_PBMM2_EXE align --split-by-sample $NO_SM_BIOSAMPLES $REF $CRAMTMP/split-no-sm.bam
|
|
0 commit comments