I MiSeqed heterogeneous templates amplified via rolling cycle amplification (RCA) by a specific 6-nt primer. PE reads were merged into single ones around 500 nt. As a consequence of RCA, some of them consist of 100-200 nt repeats from short-size templates. Is anyone aware of tools that can process, extract, and consensus these kinds of reads? Many thanks.
-
-
You could take an example sequence and then use "bbduk.sh" from BBMap suite in filter mode to extract reads containing that sequence with "literal=sequence_you_want" option. Guide to use bbduk is here. -
I believe that is what I want. Thanks!Comment
-
GenoMax: Can I extract/separate reads with 0, 1, or 2 motifs (the motif is the same)? Could BB scripts do this job?Comment
-
Possibly. You will have to play with parameters for bbduk.sh. "0" should be easy to get.Comment
-
Okay, I will try. ThanksComment
-
The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.
This convergence of genetics, immunology, and computation...Yesterday, 05:41 AM
Comment