diff --git a/src/bedtools/bedtools_bedtobam/test.sh b/src/bedtools/bedtools_bedtobam/test.sh index c5b24162..2d9b8d5c 100644 --- a/src/bedtools/bedtools_bedtobam/test.sh +++ b/src/bedtools/bedtools_bedtobam/test.sh @@ -37,12 +37,6 @@ printf "chr2\t.\texon\t1500\t1700\t.\t+\t.\tID=exon2;Parent=transcript1\n" >> "t printf "chr1\t.\tCDS\t1000\t1200\t.\t+\t0\tID=cds1;Parent=transcript1\n" >> "test_data/example.gff" printf "chr1\t.\tCDS\t1500\t1700\t.\t+\t2\tID=cds2;Parent=transcript1\n" >> "test_data/example.gff" -# Expected umcompressed SAM files -# printf "@HD\tVN:1.0\tSO:unsorted\n@PG\tID:BEDTools_bedToBam\tVN:Vv2.30.0\n@SQ\tSN:chr1\tAS:genome.txt\tLN:248956422\n@SQ\tSN:chr3\tAS:genome.txt\tLN:242193529\n@SQ\tSN:chr2\tAS:genome.txt\tLN:198295559\nmy_read/1\t0\tchr1\t129\t255\t100M\t*\t0\t0\t*\t*\nmy_read/2\t16\tchr1\t429\t255\t100M\t*\t0\t0\t*\t*\n" > "test_data/expected.sam" -# printf "" > "test_data/expected12.sam" -# printf "" > "test_data/expected_mapquality.sam" -# printf "" > "test_data/expected_gff.sam" - # Test 1: Default conversion BED to BAM mkdir test1 cd test1 @@ -60,7 +54,7 @@ echo "- test1 succeeded -" cd .. -# assert_identical_content # (bam file is compressed and -ubam option does not seem to work on this version.) +# I could not assert_identical_content (bam file is compressed and -ubam option does not seem to work on this version). # could use samtools view to check the content of the bam file, would need to add samtools as a dependency. # TODO: check if older versions works with -ubam option!!! # Version 2.27.1 (bedtools) -ubam outputs better than version 2.30.0 (bedtools). @@ -135,6 +129,7 @@ echo "> Run bedtools_bedtobam on GFF file" --input "../test_data/example.gff" \ --genome "../test_data/genome.txt" \ --output "output.bam" +# --uncompress_bam # checks assert_file_exists "output.bam" diff --git a/src/bedtools/bedtools_bedtobam/test_data/example.bed b/src/bedtools/bedtools_bedtobam/test_data/example.bed deleted file mode 100644 index dbcc2a97..00000000 --- a/src/bedtools/bedtools_bedtobam/test_data/example.bed +++ /dev/null @@ -1,2 +0,0 @@ -chr2:172936693-172938111 128 228 my_read/1 37 + -chr2:172936693-172938111 428 528 my_read/2 37 - diff --git a/src/bedtools/bedtools_bedtobam/test_data/example.gff b/src/bedtools/bedtools_bedtobam/test_data/example.gff deleted file mode 100644 index 7ddbb26f..00000000 --- a/src/bedtools/bedtools_bedtobam/test_data/example.gff +++ /dev/null @@ -1,7 +0,0 @@ -##gff-version 3 -chr1 . gene 1000 2000 . + . ID=gene1;Name=Gene1 -chr3 . mRNA 1000 2000 . + . ID=transcript1;Parent=gene1 -chr1 . exon 1000 1200 . + . ID=exon1;Parent=transcript1 -chr2 . exon 1500 1700 . + . ID=exon2;Parent=transcript1 -chr1 . CDS 1000 1200 . + 0 ID=cds1;Parent=transcript1 -chr1 . CDS 1500 1700 . + 2 ID=cds2;Parent=transcript1 diff --git a/src/bedtools/bedtools_bedtobam/test_data/expected_gff.sam b/src/bedtools/bedtools_bedtobam/test_data/expected_gff.sam deleted file mode 100644 index 72049867..00000000 --- a/src/bedtools/bedtools_bedtobam/test_data/expected_gff.sam +++ /dev/null @@ -1,11 +0,0 @@ -@HD VN:1.0 SO:unsorted -@PG ID:BEDTools_bedToBam VN:Vv2.30.0 -@SQ SN:chr1 AS:/viash_automount/Users/theo/Desktop/biobox/src/bedtools/bedtools_bedtobam/test_data/genome.txt LN:248956422 -@SQ SN:chr3 AS:/viash_automount/Users/theo/Desktop/biobox/src/bedtools/bedtools_bedtobam/test_data/genome.txt LN:242193529 -@SQ SN:chr2 AS:/viash_automount/Users/theo/Desktop/biobox/src/bedtools/bedtools_bedtobam/test_data/genome.txt LN:198295559 -gene 0 chr1 1000 255 1001M * 0 0 * * -mRNA 0 chr3 1000 255 1001M * 0 0 * * -exon 0 chr1 1000 255 201M * 0 0 * * -exon 0 chr2 1500 255 201M * 0 0 * * -CDS 0 chr1 1000 255 201M * 0 0 * * -CDS 0 chr1 1500 255 201M * 0 0 * * diff --git a/src/bedtools/bedtools_bedtobam/test_data/genome.txt b/src/bedtools/bedtools_bedtobam/test_data/genome.txt deleted file mode 100644 index bc3e0dfd..00000000 --- a/src/bedtools/bedtools_bedtobam/test_data/genome.txt +++ /dev/null @@ -1,3 +0,0 @@ -chr1 248956422 -chr3 242193529 -chr2 198295559 diff --git a/src/bedtools/bedtools_bedtobam/test_data/output.bam b/src/bedtools/bedtools_bedtobam/test_data/output.bam deleted file mode 100644 index 3d56a631..00000000 Binary files a/src/bedtools/bedtools_bedtobam/test_data/output.bam and /dev/null differ diff --git a/src/bedtools/bedtools_bedtobam/test_data/output.sam b/src/bedtools/bedtools_bedtobam/test_data/output.sam deleted file mode 100644 index f974874f..00000000 --- a/src/bedtools/bedtools_bedtobam/test_data/output.sam +++ /dev/null @@ -1,7 +0,0 @@ -@HD VN:1.0 SO:unsorted -@PG ID:BEDTools_bedToBam VN:Vv2.30.0 -@SQ SN:chr1 AS:genome.txt LN:248956422 -@SQ SN:chr3 AS:genome.txt LN:242193529 -@SQ SN:chr2 AS:genome.txt LN:198295559 -my_read/1 0 chr1 129 255 100M * 0 0 * * -my_read/2 16 chr1 429 255 100M * 0 0 * *