Add some unmatched constructs

broadinstitute · Sep 5, 2023 · f404cdf · f404cdf
1 parent f18df8a
commit f404cdf
Show file tree

Hide file tree

Showing 2 changed files with 89 additions and 1 deletion.
diff --git a/test-data/Makefile b/test-data/Makefile
@@ -191,7 +191,7 @@ test-demultiplexed:
 	--col-reference $(conditions) \
 	--row-reference $(reference) \
 	--demultiplexed \
-	--reads TTGAACCG:../../test-data/demultiplexed/TTGAACCG.construct.fastq,CCGAGTTA:../../test-data/demultiplexed/CCGAGTTA.construct.fastq,TTGAGTAT:../../test-data/demultiplexed/TTGAGTAT.construct.fastq,CCTCCAAT:../../test-data/demultiplexed/CCTCCAAT.construct.fastq,GGTCACCG:../../test-data/demultiplexed/GGTCACCG.construct.fastq,TTGACAAT:../../test-data/demultiplexed/TTGACAAT.construct.fastq,AATCCAAT:../../test-data/demultiplexed/AATCCAAT.construct.fastq,TTCTCATA:../../test-data/demultiplexed/TTCTCATA.construct.fastq,AATCCACG:../../test-data/demultiplexed/AATCCACG.construct.fastq,AATCGTGC:../../test-data/demultiplexed/AATCGTGC.construct.fastq,AAGAACTA:../../test-data/demultiplexed/AAGAACTA.construct.fastq,CCAGTGAT:../../test-data/demultiplexed/CCAGTGAT.construct.fastq,GGTCGTGC:../../test-data/demultiplexed/GGTCGTGC.construct.fastq,TTAGACCG:../../test-data/demultiplexed/TTAGACCG.construct.fastq,GGTCCACG:../../test-data/demultiplexed/GGTCCACG.construct.fastq,CCGAACTA:../../test-data/demultiplexed/CCGAACTA.construct.fastq,AACTCACG:../../test-data/demultiplexed/AACTCACG.construct.fastq,AATCACTA:../../test-data/demultiplexed/AATCACTA.construct.fastq,GGTCCATA:../../test-data/demultiplexed/GGTCCATA.construct.fastq,GGTCTGCG:../../test-data/demultiplexed/GGTCTGCG.construct.fastq,CCAGTGGC:../../test-data/demultiplexed/CCAGTGGC.construct.fastq,AACTTGCG:../../test-data/demultiplexed/AACTTGCG.construct.fastq \
+	--reads TTGAACCG:../../test-data/demultiplexed/TTGAACCG.construct.fastq,CCGAGTTA:../../test-data/demultiplexed/CCGAGTTA.construct.fastq,TTGAGTAT:../../test-data/demultiplexed/TTGAGTAT.construct.fastq,CCTCCAAT:../../test-data/demultiplexed/CCTCCAAT.construct.fastq,GGTCACCG:../../test-data/demultiplexed/GGTCACCG.construct.fastq,TTGACAAT:../../test-data/demultiplexed/TTGACAAT.construct.fastq,AATCCAAT:../../test-data/demultiplexed/AATCCAAT.construct.fastq,TTCTCATA:../../test-data/demultiplexed/TTCTCATA.construct.fastq,AATCCACG:../../test-data/demultiplexed/AATCCACG.construct.fastq,AATCGTGC:../../test-data/demultiplexed/AATCGTGC.construct.fastq,AAGAACTA:../../test-data/demultiplexed/AAGAACTA.construct.fastq,CCAGTGAT:../../test-data/demultiplexed/CCAGTGAT.construct.fastq,GGTCGTGC:../../test-data/demultiplexed/GGTCGTGC.construct.fastq,TTAGACCG:../../test-data/demultiplexed/TTAGACCG.construct.fastq,GGTCCACG:../../test-data/demultiplexed/GGTCCACG.construct.fastq,CCGAACTA:../../test-data/demultiplexed/CCGAACTA.construct.fastq,AACTCACG:../../test-data/demultiplexed/AACTCACG.construct.fastq,AATCACTA:../../test-data/demultiplexed/AATCACTA.construct.fastq,GGTCCATA:../../test-data/demultiplexed/GGTCCATA.construct.fastq,GGTCTGCG:../../test-data/demultiplexed/GGTCTGCG.construct.fastq,CCAGTGGC:../../test-data/demultiplexed/CCAGTGGC.construct.fastq,AACTTGCG:../../test-data/demultiplexed/AACTTGCG.construct.fastq,../../test-data/demultiplexed/unmatched.construct.fastq \
 	--row-barcode-policy PREFIX:CACCG@18)
 	(cd $(wd); diff counts.txt ../../test-data/expected-counts.txt && \
      diff lognormalized-counts.txt ../../test-data/lognormalized-counts.txt && \

diff --git a/test-data/demultiplexed/unmatched.construct.fastq b/test-data/demultiplexed/unmatched.construct.fastq
@@ -0,0 +1,88 @@
+@HWUSI-EAS100R:6:23:398:3989#1
+TATGTAAAGGGTAAAAGTGCAGTGCACCGCCATAATACTAGGTGACAGA
++
+=87E45/08.93:497@487-@91/-6;2.:;?475:24887;/7/8;5
+@HWUSI-EAS100R:6:23:398:3989#43
+TATGTAAAGGGTAAAAGTGCAGTGCACCGAGAGCGGATGAGGATACAGA
++
+.=964=:3659:-9/2774:855760844>8<;?956><:=;2=5,616
+@HWUSI-EAS100R:6:23:398:3989#6
+TATGTAAAGGGTAAAAGTGCAGTGCACCGCCATGATGTCTAAGCTGAGA
++
+2894/7A;2928.<<8;67391C;<:=-546<6;25.;9=58@:219A:
+@HWUSI-EAS100R:6:23:398:3989#40
+TATGTAAAGGGTAAAAGTGCAGTGCACCGGCAGATCTTTGAGATGATGG
++
+@499@5=;-63/=742:24984585675088>99?A5;.76558807>:
+@HWUSI-EAS100R:6:23:398:3989#10
+TATGTAAAGGGTAAAAGTGCAGTGCACCGCCATAATACTAGGTGACAGA
++
+0/?.;<2A9:56<3565,4555<-8A86:8=:?4:.2-47@155787).
+@HWUSI-EAS100R:6:23:398:3989#16
+TATGTAAAGGGTAAAAGTGCAGTGCACCGCAACCTGGAGACGCAGCACA
++
+3B8E707:775869/>07917463B53227=266219/@9;3D=5>;95
+@HWUSI-EAS100R:6:23:398:3989#18
+TATGTAAAGGGTAAAAGTGCAGTGCACCGCAACCTGGAGACGCAGCACA
++
+@4E73<968=C2=581::13925367741=558@<<B?-59;:=;44?.
+@HWUSI-EAS100R:6:23:398:3989#22
+TATGTAAAGGGTAAAAGTGCAGTGCACCGGTATCGTCCAAGGAATTGAT
++
+851/4:7@7A;7<50=7>1;,:98826B26121:5:=.?2471=1//A:
+@HWUSI-EAS100R:6:23:398:3989#11
+TATGTAAAGGGTAAAAGTGCAGTGCACCGGACAATGACTTTACCTATAT
++
+7<;9.48.C86:<B<;=>/:4:[email protected]<9>:4=>3</97>60>6296
+@HWUSI-EAS100R:6:23:398:3989#12
+TATGTAAAGGGTAAAAGTGCAGTGCACCGGAGTTCTGGAACGGGACTTT
++
+4637:=8=>688=073<=/*87472=6=9853?3=;-4BA5878?/D:E
+@HWUSI-EAS100R:6:23:398:3989#32
+TATGTAAAGGGTAAAAGTGCAGTGCACCGGTCCATCGTCCCAGATCTTA
++
+488;19379.74=<4<069>::3A28=52<><626@:74/25>9</3:0
+@HWUSI-EAS100R:6:23:398:3989#28
+TATGTAAAGGGTAAAAGTGCAGTGCACCGCCTTTGTACCACCAGCATGT
++
+779C71;?4:14>4742<5/18766>64;3=>54?7757223:CD=:55
+@HWUSI-EAS100R:6:23:398:3989#4
+TATGTAAAGGGTAAAAGTGCAGTGCACCGTCCTGGTATGGGAGGTTATT
++
+71:@:;4D<3867C7163;<7?;896545@7>624;4<7?-1C03?A37
+@HWUSI-EAS100R:6:23:398:3989#23
+TATGTAAAGGGTAAAAGTGCAGTGCACCGCACTACCAAATTGATACAAA
++
+8;@9;<:3?475=;<=6927;2;94:345::898;957>857<D50971
+@HWUSI-EAS100R:6:23:398:3989#3
+TATGTAAAGGGTAAAAGTGCAGTGCACCGGCCTTTGCATATCGCAGTAA
++
+07<;10:838,@667842(>:7;,9:?337B><><3>.0==<6=50:?0
+@HWUSI-EAS100R:6:23:398:3989#44
+TATGTAAAGGGTAAAAGTGCAGTGCACCGAGGGAGTCACACTGACCACC
++
+58:?7686+6274<=15=;;@339516573<2992;2;4<6/19:5397
+@HWUSI-EAS100R:6:23:398:3989#8
+TATGTAAAGGGTAAAAGTGCAGTGCACCGTCCTGGTATGGGAGGTTATT
++
+:<;5;534/<48.2->7<<73;@028993;23?=2656>3:48767>5<
+@HWUSI-EAS100R:6:23:398:3989#92
+TATGTAAAGGGTAAAAGTGCAGTGCACCGCCTTCAAATGATGCTGTAGT
++
+*78;4;;4A':915+=8822745.35:,5@5A7>4996;.=93?1B183
+@HWUSI-EAS100R:6:23:398:3989#7
+TATGTAAAGGGTAAAAGTGCAGTGCACCGGCTGAGATGATTGAACTGAA
++
+>:;<+2-965774748666=8*967A3>7,:B9447760<.63<;?227
+@HWUSI-EAS100R:6:23:398:3989#2
+TATGTAAAGGGTAAAAGTGCAGTGCACCGCCTCCGTTCTGATACTCACA
++
+961=2?:876?467=68E81=1;6888-66725<6;7875?;7@7?551
+@HWUSI-EAS100R:6:23:398:3989#13
+TATGTAAAGGGTAAAAGTGCAGTGCACCGCCATAATACTAGGTGACAGA
++
+..55725@594=69;8<3234637<@=98<96:67675?4;;,31?5:6
+@HWUSI-EAS100R:6:23:398:3989#52
+TATGTAAAGGGTAAAAGTGCAGTGCACCGCTATGAAGCTCTCACGGTTA
++
+7;/9076C71:99;.69344;91<58<2>7;/<>1247740;05@/57=