merge bam
I am trying to create a combined bam file, from 8 files, to view on IGV
I get this error
What I can do ?
base) javier@iMac-de-JAVIER BWA % samtools merge -o merged.bam V350019555_L03_B5GHUMqcnrRAABA-549.bam V350019555_L03_B5GHUMqcnrRAABA-550.bam V350019555_L03_B5GHUMqcnrRAABA-551.bam V350019555_L03_B5GHUMqcnrRAAB(base) javier@iMac-de-JAVIER BWA % samtools merge -o merged.bam V350019555_L03_B5GHUMqcnrRAABA-549.bam V350019555_L03_B5GHUMqcnrRAABA-550.bam V350019555_L03_B5GHUMqcnrR(base) javier@iMac-de-JAVIER BWA % samtools merge -o merged.bam V350019555_L03_B5GHUMqcnrRAABA-549.bam V350019555_L03_B5GHUMqcnrRAABA-550.(base) javier@iMac-de-JAVIER BWA % samtools merge -o merged.bam V350019555_L03_B5GHUMqcnrRAABA-549.bam V350019555_L03_B5GHUMqcnrRAABA-550.bam V350019555_L03_B5GHUMqcnrRAABA-551.bam V350019555_L03_B5GHUMqcnrRAABA-552.bam V350019555_L03_B5GHUMqcnrRAABA-553.bam V350019555_L03_B5GHUMqcnrRAABA-554.bam V350019555_L03_B5GHUMqcnrRAABA-555.bam V350019555_L03_B5GHUMqcnrRAABA-556.bam
samtools: invalid option -- o
the version es 1.3.1
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
I guess the size increase because the output qualities are less homogenous, ie the number of different letters coded as qualities increased and hence it's "harder" to compress.
I guess such sequence is easier to compress:
",:,:,F,::F,,FFFF,,F,,,F:F::,F,F,F,F,F,FFF:,FF,F,:,:FF,F,F:FFFFFFFFF,FFFF:FFFFFF:FFFFFFFFFFFFFFFFFF,FF"
than
"/>-=-H'@>;/-BGCF-/B0.-F>B@@.J-B0C-F/B0CFF?'GC-F/70>FF/B0C=FHBHCFFHB0CGFH7GDGGHC@CFFHBGCFFHBGCFFGDI0GE" like in my case after recalibration.
Leave a comment:
-
Perhaps this is applicable in your case as well: http://gatkforums.broadinstitute.org...-recalibration
Leave a comment:
-
I haven't used it in forever but it may be keeping the old quality scores as well as the recalibrated ones. If that is default now there is an option to have it not do that (at least there was over a year ago when I last looked at this).
Leave a comment:
-
bam file increases in size following base recalibration
I am doing some whole genomes and after the alignments, I am doing the indel-realign and base reaclibration steps. I am seeing an odd result. The bam file after the base recalibration increases in size as below:
After the indel re-align, which is the same size as the originally mapped file:
114G Aug 23 08:51 Nrml_IndelRealin.bam
96G Aug 23 21:59 umr_IndelRealin.bam
After base quality recalibration:
126G Aug 25 20:40 Nrml_IndelRealin.BQSR.bam
127G Aug 25 15:44 Tumr_IndelRealin.BQSR.bam
Is this normal? Should the file size increase?
ThanksTags: None
Latest Articles
Collapse
-
by seqadmin
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...-
Channel: Articles
04-04-2024, 04:25 PM -
-
by seqadmin
Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...-
Channel: Articles
03-22-2024, 06:39 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
22 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
24 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
52 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Leave a comment: