diff --git a/tests/cram/data/demultiplex.bc1009--bc1009.bam b/tests/cram/data/demultiplex.bc1009--bc1009.bam new file mode 100644 index 0000000..441a88e Binary files /dev/null and b/tests/cram/data/demultiplex.bc1009--bc1009.bam differ diff --git a/tests/cram/data/demultiplex.bc1009--bc1009.bam.bai b/tests/cram/data/demultiplex.bc1009--bc1009.bam.bai new file mode 100644 index 0000000..a2a94ee Binary files /dev/null and b/tests/cram/data/demultiplex.bc1009--bc1009.bam.bai differ diff --git a/tests/cram/data/demultiplex.bc1009--bc1009.bam.pbi b/tests/cram/data/demultiplex.bc1009--bc1009.bam.pbi new file mode 100644 index 0000000..db0e79e Binary files /dev/null and b/tests/cram/data/demultiplex.bc1009--bc1009.bam.pbi differ diff --git a/tests/cram/data/demultiplex.bc1009--bc1009.consensusreadset.xml b/tests/cram/data/demultiplex.bc1009--bc1009.consensusreadset.xml new file mode 100644 index 0000000..f065351 --- /dev/null +++ b/tests/cram/data/demultiplex.bc1009--bc1009.consensusreadset.xml @@ -0,0 +1,23 @@ + + + + + + + + + + + + + + + + + + + 21457684 + 5081 + + + diff --git a/tests/cram/data/demultiplex.bc1009--bc1009.ref.fasta b/tests/cram/data/demultiplex.bc1009--bc1009.ref.fasta new file mode 100644 index 0000000..8ec962e --- /dev/null +++ b/tests/cram/data/demultiplex.bc1009--bc1009.ref.fasta @@ -0,0 +1,2 @@ +>ref +TGCCTCCTAATTACTTTGAGGACCCTCCCAAGCATGCTAGATTTGCTGAGCCTGCAGCCAGCCCCATAAGGAACTCCTGCATCCTGTATGCTGAGCAGAGCACAAACACTCATAAACTGAGGCCCATGCATGGCGTGACACCACAGTTCCACACCAAAGAGCAAGAGCCATAAAAACAGAACAGATACTTTCTAAATTGCTTTCATCAATCTTACTTTGAGAAGAGGATTTTTTGGCTGAAGCACGTGAACAGAACTGCCAAAATCCAAACTTTGTGAGCAGTGTTTGTTACTCAAAATTAGTGAATACCTCCTTTCTCAACAGCAAAATACAGAGAAATATGACAAATATGCATTTTATTTAAAATAACTATAAAACTATAAATAAACCTAATGCAGTTGATTTCAGAGCATGCTAAGAACCTTTAAAAGTTATTATTCTAGCCAGGCATGGTGGCTCACACCTGTAATCCCAGCACTTTGGGAGGCCGAGGTGGACAGATCACAAGGTCAGGAGTTCAAGACCAGCCTGGCCAACATGGTGAAACCCCATCTCTACTAAAAATACAAAAATTAGCTGGACATGGTGGTGGGCGCCTGTAATCCCAGCTACTCAGAAGGCTGAGGCAGGAGAAGCATTTGAACCTGGGAGGCAGAGACTGCAGTGAGCCGAGATTGTATCACTGCACTCCAGCCTGGGTGACAATGTGAGACTACATCTCAAAAAATATATAGATATATATTATTCTAAACTACTTTTGTTTTTAGAGCATTTTAAATCTATCTTACAATATGAACAGGTTATATCTGGCAGGATACATGAGGAACTGAAATAGTGCTTGTTCTTGGGGACAAAGTTAGGTACTCAGAGGAAAGGGGTGGAAGGAAGATATTTTAATTATTGCCCTTTTGTAATTACAGAATTTGTATCATGTCAATGTATTACTTATGCAAAAATAATACATTAAAAAAAATTTTTTTTGAGATAGAGTCTCGCTCTGTCACCCAGGCTGGAGTGCAGTGGCGTGATCTCAACTCACTGCAGACTCTGCCTCCCGGGTTCAAGCGATTTTCCTGCCTCAGCCTCCCAGGTAGCTTGGATTACAGGCACGTGCCACCATGACCAGCTAATTTTTGTATTTTTAGTAGAGACAGGGTTTCACCATGTTAGCCAGGCTGGTCTCAAACTCCTGACCTCAGGTGATTCGTCCACCTAGGCCTCCCAAAGTGCTCGGATTAAAGGTGTAAGCCACCTTGCCAGTCCTGATTAAATTTTTTAAATAAGGTAATCATATATTACCTTAAAGTAATGAAACAAAATCTATCTCCTTTGTCAAGGCAGCAAAGGAAATTTCTTGATACCCATATTTTGAAGATAAGCAAATGCAGGTTCTGTCAAATCAGCATAAGAAAGAGAAGGCATATTTTTAGAATGATCAAGTGTCTGAAGCAGCCCAGCAGAAAGGTCTTCCTTCCCATCAGGAGCCCTGTTTAGTGTTTACACAAGATACATGCAATGTGTGAACTTTTATGACACCTACCTCCTAGAAAAGCACTACAGTAAATTTCTGGATGAGAAATTACTTACTGAAGAGATAAGCCAAATTGTTTTAAAATTTCTCCAAATATCCAAAAAGTTTTTCCATCCTTTAGAATAATTTATTCCTAAATTTTATTCCAACAAATAAGCTAACACTGCTCTCCGAATCATTTTCTAATTAAACTTTGAAGACGAGACAAGGATCATAACAGAAACAAAGAACAAGAAATATTACTAAAAATTAATCAGTCTTAACAGAGATATTCTTACCAAATGCCATCTGGGACACATCCAGATCGGAAATAGAATACTGAGAGGTAAGAGAAAGGCCTGAACTGCACCAGAGTACTTCCTAACTCCTACTGTACACTCGACACACAGACTTCCAGGGAAAGTCCTGGGGTTCCGCTAGATGTTTTAATGGTCTCTCTTTTCCTGTTTCAAAAATTCTCCCCAAGAACAAGCAACAGAAACCCCTAGCTTCCAAAAAGTAACAAGCTGTCAGCAATGAAGCAAACAGAAACCCTTTCACTGGCACATCTGAGATGCCCACGGAGTCTGCGGACCCCCCGTATTCTGGCAAGCCTGACTGGAATACTGGTAAAATGTACACCTCATTCATCCCTCGGCAATTAACCTCCCCATCAGCAACGTGTTGGTTTTATTCCCAACACACATAGGCAAAAAACTTCAATCAACTCTTCCTAATATTTAAGATAATTTCATTAATCACCACGTAAGTTCACCACTTTACTTGGCAACCACAACTCATTCAAACGCCTGCAGACCAACTTAGGCTTAGATGAGCAGATCTCTTAGGTGGTGCTTAAAATCGTGTTACTACCACCGTGATCATGAACTAAATGAAATGTCCTGACAATGTAGGAAGAGCCTCAAATCTCAGGCTGTTTTAAGTGCCACTAGGAAAAGAAATCTCCCATTGGACGCGCCCTCACTCCCCGCCCCCCCCCCTGCCCCGCCACTCCCGCCGAGGCGTCCTAAAAGCATTATGTCATCCACTGCCCCCACATCTCTCCAGCTCCAAACCCCGCCTCCTCTCGAAGTGTCCCGCGCCCCCAGCCCTCTTCCCTCTCAGACTAGAAGAGGGAAGTGGCACTGAGCAAATCTCGGCTCCTCCAAGTTCTCTTCACAACGCAGTTAAAAGAACCCCCGCCCTGGTTTCTCGCAAATAAAAAGAAACAAGTTCTCGCCCCAACTCTTCGCATCCTGCGAAGGGCGAGAGGGCTGAGGCCGTGACTCCCCTCGGTGAATTCAGGACAGGCCCCAACAAGGCTCTGCCTCCCCTCGCGAGAGGACAAGGGAAGACCCAAGTGAGGGAGCGGGGCTGAAGTGGGGGAAGGCCTCGCCCCAGGAGGGGGCGGGTGTCCCTCATGGGCTCTGGGTTGCTGGGTCACTCTGTCTCTGCGGGGCCGGGGGTTCGTGTCGCCGGCCCGCAGGCTGCAGGGTTACCGCCATCCCCGCCGTAGCCTGGGACCCGCCGGGACAGGGAGCTGCAGC diff --git a/tests/cram/demuxed.t b/tests/cram/demuxed.t new file mode 100644 index 0000000..abd1aec --- /dev/null +++ b/tests/cram/demuxed.t @@ -0,0 +1,8 @@ + $ IN=$TESTDIR/data/demultiplex.bc1009--bc1009.consensusreadset.xml + $ REF=$TESTDIR/data/demultiplex.bc1009--bc1009.ref.fasta + $ SAMTOOLS_VERSION=$(samtools --version | head -n 1 | cut -d '.' -f 2) + $ if [ $SAMTOOLS_VERSION -ge 10 ]; then SAMTOOLS_NOPG='--no-PG'; fi + + $ $__PBTEST_PBMM2_EXE align $IN $REF "${CRAMTMP}"/demultiplex.bc1009--bc1009.bam --log-level FATAL + $ samtools view ${SAMTOOLS_NOPG} "${CRAMTMP}"/demultiplex.bc1009--bc1009.bam | wc -l | tr -d ' ' + 1 diff --git a/tests/cram/meson.build b/tests/cram/meson.build index 843a931..8a95504 100644 --- a/tests/cram/meson.build +++ b/tests/cram/meson.build @@ -2,6 +2,7 @@ pbmm2_cram_test_data = files([ 'baminput.t', 'biosampleConsensus.t', 'biosampleSubreads.t', + 'demuxed.t', 'empty.t', 'fastxinput.t', 'idempotence.t',