deinterleave with singletons
Hi!
I have a interleaved fastq containing unmapped reads produced by segemehl -u. I want to deinterleave it into the two mate pair files as well as removing/saving the singletons into a separate file.
Currently, reformat.sh cannot deal with it, even if I give outsingle= as parameter. The header contains the strand information (i. e. 2:N:0:2).
Is there some way to get at least the pairing reads extracted without singletons in between?
--
Kind regards,
Mathias
Hi!
I have a interleaved fastq containing unmapped reads produced by segemehl -u. I want to deinterleave it into the two mate pair files as well as removing/saving the singletons into a separate file.
Currently, reformat.sh cannot deal with it, even if I give outsingle= as parameter. The header contains the strand information (i. e. 2:N:0:2).
Is there some way to get at least the pairing reads extracted without singletons in between?
--
Kind regards,
Mathias
Comment