#HISEQ1:20:H9V1RADXX:2:2103:20115:38796	163	1	535899	44	148M	=	536387	636	TGTAAACATGATTTCCTACTGGCATATCTGACTATAACTGACCACCTCAGGGTTCATTCCGATCTGTATATAAGTATCATGTAAACACGAGTTCCTGCTGGCATATCTGACTGTAACCGACCTCCTCAGGGTCCATTCTGATCTGTAT	CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJJJJJJJJJIJJIJJIJJIJJIJGEHJJJJGGHIJJJJJJJIJAHEHHEHHEHFFFFEDCCD=?A>ACDDDD<A@CDDCDCCD@@BCC>9>BB<CDCCDDCDCCBCCDEEEEEDEDDED	NM:i:0	MD:Z:148	MC:Z:148M	AS:i:148	XS:i:143	XA:Z:6,-171031411,148M,1;8,+35285,148M,1;5,-180873899,148M,1;1,-447345,148M,1;17,-81161206,148M,2;
#(base) root@4c386ca7f7bc:~/DeepMEI/data_cluster/split_softclipped_sort_666666# head tmp_map_clip.txt
#D00360:18:H8VC6ADXX:1:1109:10649:47125	1	536070
open FP,"tmp_map_clip.txt";
while(<FP>)
{
	chomp();
	$soft_read{$_}=1;
}
while(<STDIN>)
{
	if($_=~/^@/)
	{
		print "$_";
		next;
	}
	chomp();
	@line=split(/\t/,$_);
	if(exists $soft_read{"$line[0]\t$line[2]\t$line[3]"})
	{
		print "$_\n";
	}
}
