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?
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
Latest Articles
Collapse
-
by seqadmin
The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...-
Channel: Articles
07-08-2024, 03:19 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 07-25-2024, 06:46 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
07-25-2024, 06:46 AM
|
||
Started by seqadmin, 07-24-2024, 11:09 AM
|
0 responses
26 views
0 likes
|
Last Post
by seqadmin
07-24-2024, 11:09 AM
|
||
Started by seqadmin, 07-19-2024, 07:20 AM
|
0 responses
160 views
0 likes
|
Last Post
by seqadmin
07-19-2024, 07:20 AM
|
||
Started by seqadmin, 07-16-2024, 05:49 AM
|
0 responses
127 views
0 likes
|
Last Post
by seqadmin
07-16-2024, 05:49 AM
|