From 83a0886e8e29387530ca365ace5f82298963ae74 Mon Sep 17 00:00:00 2001 From: Christopher Tomkins-Tinch Date: Mon, 5 Dec 2016 15:33:04 -0500 Subject: [PATCH] remove expected sam output file, do not redirect stderr to stdout remove expected sam output file, do not redirect stderr to stdout --- ...ected_no_leading_trailing_indel_cigars.sam | 42 ------------------- tools/samtools.py | 3 +- 2 files changed, 1 insertion(+), 44 deletions(-) delete mode 100644 test/input/TestToolSamtools/expected_no_leading_trailing_indel_cigars.sam diff --git a/test/input/TestToolSamtools/expected_no_leading_trailing_indel_cigars.sam b/test/input/TestToolSamtools/expected_no_leading_trailing_indel_cigars.sam deleted file mode 100644 index 4575810cb..000000000 --- a/test/input/TestToolSamtools/expected_no_leading_trailing_indel_cigars.sam +++ /dev/null @@ -1,42 +0,0 @@ -@HD VN:1.4 GO:none SO:coordinate -@SQ SN:Mock1-1 LN:15357 AS:assembly.nix -@RG ID:APCJ3.1 DT:2016-06-30T20:00:00-0400 PU:APCJ3.1.TAGGCATG-GCGTAAGA LB:Mock1 SM:Mock1 CN:M04004 PL:illumina -APCJ3:1:2108:7502:10921 163 Mock1-1 179 70 101M = 370 292 GGTGAGGAGGGTTCAATTCCGCCGGAGACTTTAAAGTCAGCAGTCAAAGTCTTCGTTATTAACACACCCAATCCCACCACACGCTATCAGATGCTAAACTT ABAAAEEECCCBGFGGGGFFGFGGGGGGGHHHHHHHHHHHGHHHGHHHHGHHHHGHHHHHHHGHHHGGGGHHGGHHHFHGHGGGGGGFHGHHGGHHHHGHH LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:31 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2103:11736:25696 99 Mock1-1 227 70 101M = 522 396 GTCTCCGTTCTTAACACACCCAATCCCACCACACGCTATCAGATGCTAAACTTTTGCCTAAGAATAATCTGCAGTCAAAATGCTAGGGCATCTCACAGGGT ABABADBCEFCFGGFGGGFFFGHHGGGHGFHGHGGGGGGGHHHHHHGHHGHHHHHGHHHHHHHHHHHHGHHGHHHHHHHHHHGHHGHHFGHHHGHGHHGGH LB:Z:Mock1 MC:Z:101M MD:Z:4T4A91 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:2 SM:i:70 MQ:i:70 PQ:i:65 UQ:i:60 AS:i:60 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2108:7502:10921 83 Mock1-1 370 70 101M = 179 -292 GCAAAACCATATTAGATTAGCAGATAGATCACCTGAAGCTCAGATAGAACGCTGTGAGATGGATGGTTTTGAGCCTGGTACATATAGGCTGATTCCAAATG GGGHHHHHHHHHHHHHHHHGHHGHHHHHHHHHHHHHHGHHHHHHHGGGGGHHHHHHHHHHGGHHGHHHHHGHGGHHGGHGGGGGGGFFGFEFFFEECCCCA LB:Z:Mock1 MC:Z:101M MD:Z:60T40 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:1 SM:i:70 MQ:i:70 PQ:i:31 UQ:i:30 AS:i:30 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2103:11736:25696 147 Mock1-1 522 70 101M = 227 -396 ATGACCTCCCTCCAACCATAAATAATGGAACTCCTTACGTACATGCAGATGTTGAAGGACAGCCATGTGATGAGATTGAGCAATTCCTGGATCGGTGTTAC HHFEEFCGGGBFFGGHGHGHHHHHGHHGHHGGGHHGHGGHHGHHHGHHHHHGHHHHGHHGHGHHHHGHGHHGHHHHHGHFGGGGGGGGFCCCCCEE@ABAA LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:65 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1108:7756:11224 99 Mock1-1 1469 70 101M = 1688 320 TTCCAAAATGGAGTACAAGTACCAGCTGTAAGACAGGGAGGTCAAACTGACTTTAGAGCACAGCCTTTACAAGATCCAATTCAAGCACAACTTTTCATGCC AA@ACEEDEFFFGGGFGGGGGFHHHHHHHHHHHHHHGGGGGGGHHHHHGHHHHHHGHGHHGHGGHHHHHGHHHHHHHGHHHGHHGHHFHHEHHHHGHHGGG LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:0 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1108:7756:11224 147 Mock1-1 1688 70 101M = 1469 -320 GCAAGGGGCGAACACGGAAACACTTTCCCAAACAATCCCAATCAAAACGCACAGTTGCAAGTGGGCGACTGGGATGAGTAAATCACTGACATGATCAAACT FGGHFFFFGHHHFHFGGHHHHGHHGGGGHHHHHHHEGGHHHGHGGGGFHGHGHHHHGHHHGGGGGGGHGHGGHHHHHHHGGGGGGGGFGFFFEEEEAAAA3 LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:0 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1103:11564:16446 147 Mock1-1 2653 70 101M = 2520 -234 TTGAACAAAAGGTGGACAAGGTGCTTGCACAGGGCAGCATGGTGACCCAAATAAAGAATGAATTATCAACAGTAAAGACAACATTAGCAACAATTGAAGGG GGGGGGHHHHFFHFGHGHHHFGFGHHGHGHHFE3HHHGHHFHHHGFFGHHGHGGHHHHGHHHGHHHGGGHHGHHHHHHGGGFGFGGFGGFFFEFFEBBAAA LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:1 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2107:14341:17176 163 Mock1-1 5018 70 50M1I50M = 5208 291 CAATACAAGACCATATCAATACTATTATGAACACCCAATTGAACAATATGTCTTGTCAGATCCTTGATAACCAGCTTGCAACCTCCCTAGGATTATACCTA ACCCBDEEFFEFGGGGGGGGGFHHHHHHHHHGHGGGGHHHHHHHHHHHGHHHHHGHHHHHHGGHIHHHHHFEHHHHHHGHHGFHHHGHHGFHHHHHHGGHH LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:1 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2107:14341:17176 83 Mock1-1 5208 70 101M = 5018 -291 TATGACACCTGCAGTGGTTCAAGCAACATTATCTACCTCAATTTCTGCTGCTGAAATACTAAGTGCCGGTCTAATGGAGGGTCAGATAGTTTCTGTTCTGC HHHGHHFHHHFGHFHHHHHHHHGHHHGHHHHHHHGHHHHHHHHHGHHGHHGHHGGHHHHHHHGGGGGHHHHHHHHGGHHFFGGGGGGGFFFFFFDEBCCBC LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:1 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1103:22054:22192 163 Mock1-1 5340 70 101M = 5604 365 CATTCCAACCATTGTCACACAATCAAATGCATTGGTGATTGACTTCTACTCAATTTCGAGTTTTATTAATAATCAAGAATCCATAATTCAATTGCCAGACA ACCCBEFFFDDFGGGGGGFFGGHGHHHHHGHHHGGHGHHHGGGHHHHHHHHHHHHHHGGGGFHHGHHGHHHHHGHGHHHHGHGHHHHHEHHHHGGGHGGFG LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:2 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1103:22054:22192 83 Mock1-1 5604 70 101M = 5340 -365 AGGGAGTTATATGAGGCGATTTGTAGCACTGGATGGAACAATTGTTGCAAACTGTCGAAGTCTAACGTGTCTATGCAAGAGTCCATCTTATCCTATATACC GGGGHGHHHHHHEGGGGHHHHHFGHHHHHHGHHHHGHHGHHGHHHHHGHHHGGFHGGHHGHGHHHGHHHHHHHHHHHHHGFGGGGGFGGFFEFFEEAAABA LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:2 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1113:17899:25240 163 Mock1-1 6364 70 101M = 6630 367 GCACTGTACCAGCGATCCTGCTCTCGCTGGGGTTTTGATCAATCACTCTAGAAAGATCCCCAGTTAGGACAAGTCCCGATCCGTCACGCTAGAACAAGCTG AA@@BEEFEEFFFGGGGGGGFGHHHGGGGGGGEGGGHGGHGHHHHGHHGHGHHHHHHGGGFHFHHHGHHHHHHGHHHGGGFGHGHGGGFGGFHHFHHHGHG LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:2 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1113:17899:25240 83 Mock1-1 6630 70 101M = 6364 -367 TTGCACCTGGACCTGTTGTTAATGCGGCTGGTAAGAAGACATTCCGAACCAGTTTCCGAATATTGGTCCTATCTGTACAAGCCGTTACCCTTATATTGGTT HFBF3G33GA3C2F?GF3G3E=0E@FBADF55GD5FEED31FAEA5F5G55BGFDGF5FDGCF2GB3GA@5D455FAGGGFGGGGGGGGFFFEFEBBAABB LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:2 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1103:15883:9960 129 Mock1-1 7517 70 101M = 7764 248 GGAATATATTTCGAAAATAAGTTGATCTTTCCTGCATATGGGGGTGTCTTGCCCAATAGTACACTAGGAGTTAAATCAGCAAGAGAATTTTTCCGGCCCGT CCBBBFEFFFFEGGGGGGFFGGHHGHHHHHGGHGHHHHHGGHGGGFGHHHGHGGHHHHGHHGHFHGGHHHHHHHHGHHHGHGHHGHHHHHHHGGFFFFFFF LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 NM:i:0 MQ:i:70 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1103:15883:9960 65 Mock1-1 7764 70 101M = 7517 -248 AGCTAGTTGTCCCCTCAAACAATCAGACACTGATGGGTGCAGAAGGAAGAGTTTTATTGATCAATAATCGGCTATTATATTATCAGAGGAGTACTAGCTGG BBBBBEEEFFFFFFGFGGGFGGHHHHHGHGHHGHHHHHGFHHHHHHHHHHGHHHHHHHHGHGHHHHHHHFGFGGHHGHHHHHHHHHGHFHFHHHHGHGHGG LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 NM:i:0 MQ:i:70 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2101:18746:24157 99 Mock1-1 8318 70 101M = 8503 286 AGATTGACCATCACTTGAGTTGTAGTGAATGTAGCAGGAAGCTTTATGGGCGTGTCTCATTTCTTATCGATTATTAAGAAAAAACAGGCCAGAATGGCGGG AA?@@EEEFFFEFEFGFGFBGFAGFFFHHHHHHHHFGEGHGGHFGHGFGHHEGFGGHHFHHHGFHHHFE1FEEHFHHHHHHHHCGHHHFFHFHHHFFFFB? LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:1 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2101:18746:24157 147 Mock1-1 8503 70 101M = 8318 -286 GTTACCAAATGATTTGGAGCCAGATGACTTGGGCCCATTAGCAAATCAGAATTGGAAGGCAATTCGAGCTGAAGAATCCCAGGTTCATGCACGGTTAAAAC FGHGHFHHHHGGHFFGHHGEHHGHHHHGAGFGGF5HHHHHGFFGHHGGAHHHHHHHHHGFHHFGEHHGGHCHHHEGFGFFGFFBGGGGFAADEEEEA@@>= LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:1 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1102:8829:24118 163 Mock1-1 10946 70 101M = 11078 233 GGTCGAATCTTGACTCAAGCGTTAAAGAACGTGAGCAAGATGTGCTTAACAGCTGATATACTGGGGGATTGTTCACAAGCATCATGCTCCAATTTAGCTAC AAA?@CBEBFFDFFFBGFGGFGGGGHHHHGDAFGEHFHGEHGFGHHEFHGHFGHGHGHHHFGHGFFFFGHHHGGHHHGGFGHGHHGFGFFHHHHHGGFGGE LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:1 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1102:8829:24118 83 Mock1-1 11078 70 101M = 10946 -233 GAGAAAGATTTGTGTTATTTCCTAAATGCATTCATGACAATTAGACAATTATGTTATGATCTAGTATTTCCCCAAACTAAATCTCTTAGTCAGGACATTAC GGFHGHFGGGHFHEHHHHHHGHHGGGHHHHHHGHHGHGHGGHFHHGHHHHHHHHHHHGGHHFFGHHEGGGGGHHHGHHFGGFGGFGGGEFFFFFFDAAAAA LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:1 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2110:27666:9079 177 Mock1-1 13147 70 3S98M = 13498 355 AGAAATGTGGGGAACCAAGCAGGTGTTGACCAACATTTCTCAAGGTATCGATTATGAGATAGTTGTTCCTTCTGAAAGCCAACTTACACTCAGTGATAGAG HFHFHGFHG@3EGFGGFFGFGDFHGGA51HHFHHHFHAHEBHHGGHHHGHHHHGGGEBGGHGGCFHGGGFHHHHHGHGGGGGGGGGGGGFFFEFEDBA@A@ LB:Z:Mock1 MC:Z:101M MD:Z:98 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 NM:i:0 MQ:i:70 UQ:i:30 AS:i:30 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2110:27666:9079 113 Mock1-1 13498 70 101M = 13147 -355 TCGAGACTCTGAAGAAGGACAATTCCTGATTGCATCCTATTATAACAGTTTTGGATATCTGGAACCGATATTAATGGAATCTAAAATATTCAATCTAAGTT HGGHHHHHHHFHHFHHGGHHHHHFFHHGHHHGHHHHHIIIIIGHHGHHHGHHHHGHHHHHFFGGGGHHHHHGHHHGGHHGGGGGGGGGGFFFEFEEBAACA LB:Z:Mock1 MC:Z:3S98M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 NM:i:0 MQ:i:70 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2111:11117:17348 99 Mock1-1 13725 70 101M = 13815 191 CTCAACCCCCACTTCATCACGTTCTCAGACCACTAGGTTTGTCATCCACCTCATGGTATAAAACAATCAGTGTTTTGAATTATATTAGCCATATGAAGATA ADBBCDDDEDCEGGFGGGGFGGHGGGHHHGGHHHHHHHHHGHHHHHHHHHHHHHGGHHHHHHHGHHHHHHHHGHHHHGHHHHHHHHHHFFHHHHHHHHHHH LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:2 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2111:11117:17348 147 Mock1-1 13815 70 101M = 13725 -191 ATATGAAGATATCTGACGGTGCCCATCTATACTTGGCAGAGGGAAGTGGAGCCTCTATGTCACTTATAGAGACTTTCTTGCCCGGGGAAACCATATGGTAC HHHGGHHGHHHHGGGHGGGHGHGHHHHHHHGHHHHHHHGHGGGHHGHGGHGGHHHHHGGHGGHHHHHHGHFHHHHHGGGFGGGGGFFGGFFFFFEE@ABAA LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:2 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1113:22664:16558 99 Mock1-1 14556 70 101M = 14687 232 AAGTTATATCAGAGTGTGTCAGAGATCAATATCTGGCAGACAACAACATTATCCTCCAGGCGGGAGGTACTCCAAGCACAAGAAAATGGTTGGATCTGCCT BBCABEEFEFFFGGGGGGGFGGHHHHGHHHHHHHHHHHGGGHHHHHHHGGHHHHHGHHHGGGGGGGGEHGHGGHHHGHGHHHHHHHHHGHHFGHHGHHGGH LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:0 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1107:13139:14438 163 Mock1-1 14602 70 101M = 14839 338 CATTATCCTCCAGGCGGGAGGTACTCCAAGCACAAGAAAATGGTTGGATCTGCCTGACTATTCTTCGTTCAATGAACTACAATCTGAAATGGCCAGACTCA @CBBBDEEFEEFGGFGGGGGGEHHHHHHHHGHHHHGHHHGHHHHHGGHHHHGHHHGHHHHHHGHHGGHHGHHHHHHHHHFHHHHHGGHHHGGGGHGHFGFH LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:1 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1113:22664:16558 147 Mock1-1 14687 70 101M = 14556 -232 GAAATGGCCAGACTCATAACAATTCATCTTAAAGAGGTAATAGAAATCCTAAAGGGCCAGGCATCAGATCATGACACCCTATTATTTACTTCATACAATGT FHHHHGGGGHFHGHGHHHHHHHHHHHHHHHHHHHHGGIHHHHHIHHHGGHHGGHGGHHHGGHHHHHHHHHHHHGFGHHHGGGGGGGGGGFFFFFEECBBAB LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:0 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1107:13139:14438 83 Mock1-1 14839 70 101M = 14602 -338 TCTTATGTATACTGTGAGGAACTGGTGTATCTTGCCTTCCCAAACTCGTCTCACCTTACGACAATCTATCGAGCTTGGAGAGTTTAGACTAAGAGATGTGA GGHHHHGHHHHHHFHHHGGHHHHGHHHHHHGHHGGGHHFFGHGGFGFGHGHHHHHHHGGGHHHHHGHHGGGHFHHHHHHFGFGGGGFGGFFFFEEDCBBAC LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:1 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2105:7518:11719 177 Mock1-1 14946 70 101M = 14977 32 CCATGGAGATTCTAAAACTATCCCCCAACAGGAAATATCTAAAGTCTGCATTAAATCAATCAACATTCAACCATCTAATGGGAGAAACATCTGACATATTG 3@G535F5GD33@GGGBFHEGGEGHHHHFHFG5HGGHGGHHHFFDEHEHHGHGGHFHGGBGGGFGBFEHHHFHHHHFHHFFFGGGGGGGFFFEEEDAABB@ LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 NM:i:0 MQ:i:70 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2105:7518:11719 113 Mock1-1 14977 70 101M = 14946 -32 GAAATATCTAAAGTCTGCATTAAATCAATCAACATTCAACCATCTAATGGGAGAAACATCTGACATATTGTTAAACCGAGCTTATCAGAAGAGAATTTGGA EFHHHHFFHHHHFHHHGGHHHHHHHHHHHHHHHHHHGHHGHGGGHHHHHHHGHGGHHHHHHGHHGHHGHGHHGGGGGGHFGGGGGGGFGFEFFEEEBBAAB LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 NM:i:0 MQ:i:70 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1106:11985:11287 129 Mock1-1 14981 70 101M = 15013 33 TATCTAAAGTCTGCATTAAATCAATCAACATTCAACCATCTAATGGGAGAAACATCTGACATATTGTTAAACCGAGCTTATCAGAAGAGAATTTGGAAAGC AAA?AEEFFFEFGFFGGGGFGFHHHHHHHHHHHHHHHGHHHGHHHFHHGGHHHFHGHHHHHHGHHGHHHHHGGGBGGGHHHGFGHHGHGHHFGFEFFHHGF LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 NM:i:0 MQ:i:70 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1106:11686:10415 99 Mock1-1 14982 70 101M = 15053 172 ATCTAAAGTCTGCATTAAATCAATCAACATTCAACCATCTAATGGGAGAAACATCTGACATATTGTTAAACCGAGCTTATCAGAAGAGAATTTGGAAAGCC BBBDDEEEEFFFFGGGGGGGFGHHHHHGHHHGHHHHGHGHHHHHHHHGGHHHHHHHGIIIIHHIHHHHHHGGGGGGGHHHFHHHHGHHHHHHHHHHHHGGG LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:2 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1106:11985:11287 65 Mock1-1 15013 70 101M = 14981 -33 CAACCATCTAATGGGAGAAACATCTGACATATTGTTAAACCGAGCTTATCAGAAGAGAATTTGGAAAGCCATCGGGTGTGTAATCTATTGCTTTGGTTTGC ABBA@BEFEFFFGGGGGGGGFGHHHHHGHHHHHHHHHHHGGGGGGGGHHHHGHHGHGHHHHHGGGHHHGGHHGGGGFEGGHHHHHHHHHHEHHHHHHHHFE LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 NM:i:0 MQ:i:70 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1105:10708:23204 99 Mock1-1 15033 70 101M = 15050 118 CATCTGACATATTGTTAAACCGAGCTTATCAGAAGAGAATTTGGAAAGCCATCGGGTGTGTAATCTATTGCTTTGGTTTGCTCACCCCGGATGTTGAAGAT ACCACEEFFFFFGGGGGGGFFGGGGGHHHGHHHHHHHHHHHHHHHHHHGHHHHFFFHFFFHHHHHHFHHHGHHHHHHHHGFHGHGHHEGGGGGHHHHHHHH LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:5 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2112:16035:6267 163 Mock1-1 15036 70 101M = 15092 157 CTGACATATTGTTAAACCGAGCTTATCAGAAGAGAATTTGGAAAGCCATCGGGTGTGTAATCTATTGCTTTGGTTTGCTCACCCCGGATGTTGAAGATTCT @AAB@EEFFFFFGBGGGGCAEFGHHHHHGBDFFA23GGHHFGFFEGHGGHG?EHAGAGGHGGHGGHGHHHHHHHHHEGHGHGGHGGGGFFHHHHGFHHHGH LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:3 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1105:13024:17729 163 Mock1-1 15044 70 101M = 15164 224 TTGTTAAACCGAGCTTATCAGAAGAGAATTTGGAAAGCCATCGGGTGTGTAATCTATTGCTTTGGTTTGCTCACCCCGGATGTTGAAGATTCTGAGCGCAT BBBBBFFFEFCCGFGGGGFGGGHHHHHHHHHHHHHHHHHGHHGGGHGGGHHHHHHHHHGHHGHGGHHHGHHGHGGHGGGGGGHHGHHHHHHHGGHGFFEGG LB:Z:Mock1 MC:Z:64M3D37M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:7 NM:i:0 SM:i:70 MQ:i:70 PQ:i:100 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1105:10708:23204 147 Mock1-1 15050 70 101M = 15033 -118 AACCGAGCTTATCAGAAGAGAATTTGGAAAGCCATCGGGTGTGTAATCTATTGCTTTGGTTTGCTCACCCCGGATGTTGAAGATTCTGAGCGCATTGATGT GGFFGHGGHHHHGHGHHHHHHHGHHHGHHHHGGGGGFEFHGHGHHHHGHHHHHHHHHGHGGGHHGGGGGGGHGHHHHHGGGFGGGGGFGCCFFEFFABBAB LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:5 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1106:11686:10415 147 Mock1-1 15053 70 101M = 14982 -172 CGAGCTTATCAGAAGAGAATTTGGAAAGCCATCGGGTGTGTAATCTATTGCTTTGGTTTGCTCACCCCGGATGTTGAAGATTCTGAGCGCATTGATGTTGA FFGFGHHHHGHFHHGHGHHHHHFGHHHGGGGGGGGGHGHGHHHHGHHGHHGHHHGGHHHHHGGGGGFFHHHHHHHHHHGGGGGGGGFGFFFFFEEEABBAB LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:2 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:2112:16035:6267 83 Mock1-1 15092 70 101M = 15036 -157 GTAATCTATTGCTTTGGTTTGCTCACCCCGGATGTTGAAGATTCTGAGCGCATTGATGTTGATAATGACATACCCGATTATGATATTCACGGGGACATAAT FHHHHHFHHHFGHHHGGHHHGHFGEGGFGGGHHHHHHHHGHHHHGGGGGHHHHHHHHGHHGHHHHHHHHHGGGGGGGDFGGFGFGGCGGBBFEEDE@AAAA LB:Z:Mock1 MC:Z:101M MD:Z:101 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:70 NM:i:0 SM:i:70 MQ:i:70 PQ:i:3 UQ:i:0 AS:i:0 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA -APCJ3:1:1105:13024:17729 83 Mock1-1 15164 70 64M3D37M = 15044 -224 CCCGATTATGATATTCACGGGGACATAATTTAAATCGACTAAAGACTCCTCTGGCATTACACATCAAAAAGTGCCAAACTAGCATCCAAATTCTTCTAAAC GGGGHHHHHHHHHGGFGGGGFFHHHHHHHHHGHGGGHHGIIHIGHFHFHHGHGHHHHHHHHHHGGHHHHHGHGGHHHHGGGFGGGGGGGFFFFFEEBBBBB LB:Z:Mock1 MC:Z:101M MD:Z:64^CAC37 PG:Z:MarkDuplicates.3 RG:Z:APCJ3.1 AM:i:7 NM:i:3 SM:i:7 MQ:i:70 PQ:i:100 UQ:i:100 AS:i:100 PU:Z:APCJ3.1.TAGGCATG-GCGTAAGA diff --git a/tools/samtools.py b/tools/samtools.py index 1f4373e4c..e42ecbbf8 100644 --- a/tools/samtools.py +++ b/tools/samtools.py @@ -154,8 +154,7 @@ def filterByCigarString(self, inBam, outBam, # an input samtools process to read a bam file # bufsize 0=unbuffered, 1=linebuffered, -1=system buffered (usually fully) in_process = subprocess.Popen(in_args, bufsize=1, stdout=subprocess.PIPE, - universal_newlines=True, - stderr=subprocess.STDOUT) + universal_newlines=True) # an output samtools process to write filtered output via stdin out_args = [self.install_and_get_path(), 'view', '-h', '-o', outBam]