I found a strange read pair from IGV, so i want to extract the original information about the read pair.
I use " zcat file.fastq.gz | grep -B 3 "TCTCAAAAAGCACATGGATTCTACTCCTGATTTCTCCAAGGATCTACATATGATCAAAACATATCTTTACCAGATTCTCCGTGGAATTGCGTATTGCCACTCTCATAGGGTTCTCCATCGTGATC"
The strange thing is i can't find it in original fq file.
I use " zcat file.fastq.gz | grep -B 3 "TCTCAAAAAGCACATGGATTCTACTCCTGATTTCTCCAAGGATCTACATATGATCAAAACATATCTTTACCAGATTCTCCGTGGAATTGCGTATTGCCACTCTCATAGGGTTCTCCATCGTGATC"
The strange thing is i can't find it in original fq file.
Comment