[go: nahoru, domu]
Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Added tests #85

Merged
merged 9 commits into from
Sep 11, 2019
Merged

Added tests #85

Changes from 1 commit
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
607a100
Working on tests.
bsipos Sep 10, 2019
e642cd4
Merge branch 'master' of github.com:bsipos/seqkit
bsipos Sep 10, 2019
26b1617
Added tests for bam and fish.
bsipos Sep 10, 2019
d01e24a
Extended tests for fx2tab
bsipos Sep 10, 2019
61e9762
Merge branch 'master' of https://github.com/shenwei356/seqkit
bsipos Sep 10, 2019
5b599c8
Trying to fix missing csvtk problem in tests.
bsipos Sep 10, 2019
6e9d6d9
Fix.
bsipos Sep 10, 2019
c4d061a
+
bsipos Sep 10, 2019
7af76b2
+
bsipos Sep 10, 2019
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Next Next commit
Working on tests.
  • Loading branch information
@bsipos
bsipos committed Sep 10, 2019
commit 607a1004ccd9b88583605a69eaac10138f22847f
78 changes: 78 additions & 0 deletions tests/test.sh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -371,6 +371,8 @@ rm t.sort.*
BAM=tests/pcs109_5k.bam
PRIM_BAM=tests/pcs109_5k_prim.bam
PRIM_NANOPLOT=tests/pcs109_5k_prim_bam_NanoPplot.tsv
PRIM_WUB=tests/pcs109_5k_bam_alignment_length.tsv
PRIM_WUB_CLIP=tests/pcs109_5k_bam_soft_clips_tab.tsv

# accuracy
fun(){
Expand All @@ -384,6 +386,82 @@ R=$(cut -f 3 corr.tsv)
assert_exit_code
rm corr.tsv seqkit_acc.tsv

# MeanQual
fun(){
$app bam -f Read,MeanQual $PRIM_BAM 2> seqkit.tsv
paste seqkit.tsv $PRIM_NANOPLOT > joint.tsv
csvtk corr -t -f MeanQual,quals joint.tsv 2> corr.tsv
}
run bam_mean_qual fun
R=$(cut -f 3 corr.tsv)
(( $(echo "$R > 0.99" | bc -l) ))
assert_exit_code
rm corr.tsv seqkit.tsv

# MapQual
fun(){
$app bam -f Read,MapQual $PRIM_BAM 2> seqkit.tsv
paste seqkit.tsv $PRIM_NANOPLOT > joint.tsv
csvtk corr -t -f MapQual,mapQ joint.tsv 2> corr.tsv
}
run bam_map_qual fun
R=$(cut -f 3 corr.tsv)
(( $(echo "$R > 0.99" | bc -l) ))
assert_exit_code
rm corr.tsv seqkit.tsv

# ReadLen
fun(){
$app bam -f Read,ReadLen $PRIM_BAM 2> seqkit.tsv
paste seqkit.tsv $PRIM_NANOPLOT > joint.tsv
csvtk corr -t -f ReadLen,lengths joint.tsv 2> corr.tsv
}
run bam_read_len fun
R=$(cut -f 3 corr.tsv)
(( $(echo "$R > 0.99" | bc -l) ))
assert_exit_code
rm corr.tsv seqkit.tsv

# ReadAln
fun(){
$app bam -f Read,ReadAln $PRIM_BAM 2> seqkit.tsv
paste seqkit.tsv $PRIM_NANOPLOT > joint.tsv
csvtk corr -t -f ReadAln,aligned_lengths joint.tsv 2> corr.tsv
}
run bam_read_aln fun
R=$(cut -f 3 corr.tsv)
(( $(echo "$R > 0.99" | bc -l) ))
assert_exit_code
rm corr.tsv seqkit.tsv

# RefAln
fun(){
$app bam -f Read,RefAln $PRIM_BAM 2> seqkit.tsv
paste seqkit.tsv $PRIM_WUB > joint.tsv
csvtk corr -t -f RefAln,aligned_ref_bases joint.tsv 2> corr.tsv
}
run bam_ref_aln fun
R=$(cut -f 3 corr.tsv)
echo $R
(( $(echo "$R > 0.99" | bc -l) ))
assert_exit_code
rm corr.tsv seqkit.tsv

# LeftClip
fun(){
$app bam -f Read,LeftClip $BAM 2> seqkit.tsv
paste seqkit.tsv $PRIM_WUB_CLIP > joint.tsv
csvtk corr -t -f LeftClip,ClipStart joint.tsv 2> corr.tsv
}
run bam_left_clip fun
R=$(cut -f 3 corr.tsv)
(( $(echo "$R > 0.99" | bc -l) ))
assert_exit_code
rm corr.tsv seqkit.tsv




# ------------------------------------------------------------
# faidx
# ------------------------------------------------------------
Expand Down