Tweet
Dear all,
Using bwa, I would like to select only reads that I am very sure about ...
it is to say, reads that have unique match (mismatch would be allowed as SNP are too frequent to consider the no-mismatch option, however, I am considering putting a limit but I don't really know which one, my reads are 75 bp illumina for the moment.)
I am trying to think to some solution regarding paired-ends : should both paired-ends be unique match, only one ? to be kept as very sure mapped reads.
I am not quite sure if I should play with bwa option to restrict the search or if I should use the sam format (I will probably need some help with the flag thing, it has been quite a while since I had a look at this) to select afterwards (if we consider that calculation time is not a problem ...) .
I would really appreciate comments, ideas and help
thanks a lot
Using bwa, I would like to select only reads that I am very sure about ...
it is to say, reads that have unique match (mismatch would be allowed as SNP are too frequent to consider the no-mismatch option, however, I am considering putting a limit but I don't really know which one, my reads are 75 bp illumina for the moment.)
I am trying to think to some solution regarding paired-ends : should both paired-ends be unique match, only one ? to be kept as very sure mapped reads.
I am not quite sure if I should play with bwa option to restrict the search or if I should use the sam format (I will probably need some help with the flag thing, it has been quite a while since I had a look at this) to select afterwards (if we consider that calculation time is not a problem ...) .
I would really appreciate comments, ideas and help
thanks a lot
?
Comment