How to obtain .seq in BIC-seq2.1.1

Question

[seq-seek]

Well, first of all I have to admit that I'm a freshman to bioinformatics.

I'm exploring BIC-seq-perl_1.1.2. My question is that how to use the SAMgetUnique attached to obtain the uniquely mapped read. And does this format as chrosome #.seq? I even can't run the example data...

Thanks!

[John Wick]

hello,you can use samtools motified by BICseq2:http://compbio.med.harvard.edu/BIC-seq/BICseq2/samtools-0.1.7a_getUnique-0.1.3.tar.gz and then input the command like : [modifiedSamtools_dir]/samtools view -U bwa,output,N,N mybam.bam maybe this can help