Skip to content
Snippets Groups Projects
Commit 64ed292d authored by Aaron's avatar Aaron
Browse files

add omitted test dirs for "make test"

parent 17053b8d
No related branches found
No related tags found
No related merge requests found
File added
@HD VN:1.0 GO:none SO:coordinate
@SQ SN:chr1 LN:1000
one_blocks 16 chr1 1 40 30M * 0 0 GAAGGCCACCGCCGCGGTTATTTTCCTTCA CCCDDB?=FJIIJIGFJIJHIJJJJJJJJI MD:Z:50
\ No newline at end of file
File added
@HD VN:1.0 GO:none SO:coordinate
@SQ SN:chr1 LN:1000
one_bp_del 16 chr1 1 40 10M1D10M * 0 0 GAAGGCCACCGCCGCGCCGC CCCDDB?=FJIIJIGIIJIG MD:Z:50
BT=../../bin/bedtools
check()
{
if diff $1 $2; then
echo ok
else
echo fail
fi
}
###########################################################
###########################################################
# BAM files #
###########################################################
###########################################################
samtools view -Sb one_block.sam > one_block.bam 2>/dev/null
samtools view -Sb two_blocks.sam > two_blocks.bam 2>/dev/null
samtools view -Sb three_blocks.sam > three_blocks.bam 2>/dev/null
samtools view -Sb sam-w-del.sam > sam-w-del.bam 2>/dev/null
##################################################################
# Test three blocks without -split
##################################################################
# echo " coverage.t1...\c"
# echo \
# "chr1 0 50 1" > exp
# $BT coverage -abam three_blocks.bam -b three_blocks_nomatch.bed > obs
# check obs exp
# rm obs exp
rm *.bam
\ No newline at end of file
File added
@HD VN:1.0 GO:none SO:coordinate
@SQ SN:chr1 LN:1000
three_blocks 16 chr1 1 40 10M10N10M10N10M * 0 0 GAAGGCCACCGCCGCGGTTATTTTCCTTCA CCCDDB?=FJIIJIGFJIJHIJJJJJJJJI MD:Z:50
chr1 0 50 three_blocks_match 0 + 0 0 0 3 10,10,10, 0,20,40,
chr1 10 60 three_blocks_nomatch 0 + 0 0 0 3 11,10,10, 0,20,40,
chr1 10 60 three_blocks_nomatch 0 + 0 0 0 3 10,10,10, 0,20,40,
File added
@HD VN:1.0 GO:none SO:coordinate
@SQ SN:chr1 LN:1000
two_blocks 16 chr1 1 40 15M10N15M * 0 0 GAAGGCCACCGCCGCGGTTATTTTCCTTCA CCCDDB?=FJIIJIGFJIJHIJJJJJJJJI MD:Z:50
\ No newline at end of file
chr1 10 20 a1 1 +
chr1 30 40 a2 2 +
chr1 40 50 a3 3 -
chr1 45 100 a4 4 +
chr2 10 20 a1 5 +
chr2 30 40 a2 6 +
chr2 42 50 a3 7 +
chr2 45 100 a4 8 -
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment