hi folks, I have 5 illumina sequencing fastq files. 1 is the parent and 4 are mutants. I alligned the sequence to the refenrence genome I have all the BAM bed etc etc files. I cleaned repetitive seq in the bam file with picard. what I would like to find are the variants (SNv and indels) that are present in each of the mutant strain but not in the parent. also I would like to have the variants just in the coding sequences. what program can I use?
-
-
You need just step 1 from this tutorial since you only want SNv. -
isnt this command comparing a bam file to the reference genome? what I want are the SNV/indels present in the mutant strain but not in the control strain.Comment
-
You are identifying changes w.r.t. a reference and then comparing those results to get the difference between your control and mutant results.Comment
-
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...09-01-2026, 05:41 AM
Comment