I am novice to NGS data analysis. I use SHRiMP2 to align SOLiD mate pair reads. Because I don't need the mate-pair information but just need the mapping information, I did the alignment of F3 / R3 reads separately. Each read has only the best hit recorded, with the parameter '--strata -o 1'. The read name in the obtained sam files didnot have the postfix, i.e., '_F3' and '_R3'. Then I used 'samtools view -bS' to convert SAM to BAM format and used 'samtools merge' and 'samtools sort' to do the merge and sort process.
The question is, the size of the merged file is nearly two fold of the sorted file. I guess that, because the reads of the same mate pair have identical names, only one of them is left in the sorted BAM file, am I right? How can I keep both reads in the sorted file?
-
Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...09-23-2024, 06:35 AM -
During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.
Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...09-09-2024, 10:59 AM