-
BWA and Velvet : Anamolies in contigs and read mapping to reference
Hi All
I believe I have an interesting problem to look at or may be its just me who is experiencing this for the first time.
I have reads from GA for a bacterial genome. I used BWA to map them to known reference. Only 4% reads map back. Seeing the low % I decided to de novo assemble the reads using Velvet .
1. Total of 2150 contigs
2. 601 align to reference using 52 % reads from the original set.
3. k-mer size I used was 29. (read length : 50 bp PE)
I am not sure what to infer from these very contrasting numbers according to me. With Velvet contigs I am able to use 52% of the reads that aligns well to reference whereas with BWA only 4% reads align.
Looking forward for your thoughts.
Thanks!
-Abhi
-
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...04-22-2024, 07:01 AM -
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...04-04-2024, 04:25 PM
Leave a comment: