Entering edit mode
8.7 years ago
Louis Kok
▴
30
Hi all, Just wish to get your opinions. If the Sanger forward and reverse reads were base called separately, performed BLASTN against NCBI nt database and the result was found to be different species belonging to the same genus, how would you decide which result to take if both reads have clean peaks shown in chromatogram? Thanks a lot.
Are they exactly the reverse complementary of each other?
Or the forward sequence contains mainly the left part of the sequence, and the reverse one mainly the right part of the sequence?
Hi Antonio, using CAP3 to assemble forward and reverse reads, even with very low overlap percent identity cutoff (65%), still couldn't produce merged sequence. I suppose they are not overlapped each other.
You mention 65% of identity, but not the extension of that identity.. That could eventually join unrelated sequences, and therefore different sequences are found by Blast.
If you have sequenced the same samples in forward and reverse orientations, and you mention you got clean peaks, I am also wondering why you have used such as low cutoff value..
My first approach would be to increase that cutoff value with a more stringent value. CAP3 will render you singletons and contigs accordingly.
If different results are obtained, evaluate the values of E (the lower the better) and the query extension (the higher the better).