failed to convert from cram to bam
I got a cram file, and want to convert it to bam file. Here's the command I use and error I got
samtools view -b -T human_g1k_v37.fasta lucky.cram > lucky.bam
[sam_header_read2] 84 sequences loaded.
[sam_read1] reference 'Wm�;�����e ��i���f�[�=���~�����U��/���>���Q���O ����n����O��' is recognized as '*'.
Parse error at line 1: invalid CIGAR character
I don't know what to do with this. Please help~
• 5,663 views
•
link
updated
by
GenoMax
•
written
by
kygppyou •
•
written
by
kygppyou •
1 answer
You are using a prehistoric version of samtools, one that predates the CRAM format by several years, so of course it has no support for CRAM files. (The error messages shown indicate that it is attempting to parse the input as a SAM file, and failing miserably.)
You need at least samtools 1.0 to read CRAM files, and many CRAM-related improvements were made in 1.2. Of course, what you really want to install is the most up-to-date version of samtools, which is 1.3.1 at the time of writing.
• 0 views
•
link
written
by
John Marshall
Log in to answer this question.
More posts like this
-
sam file error
written by Chironex •
-
Convert SAM to BAM error
written by BeeWork •
-
samtools sort: truncated file. Aborting
written by m.t.lorenc •
-
How to rectify error while sorting the sam file
written by sunnykevin97
-
Converting .sam to .bam with samtools and I get a new error with every repeated command
written by ryanzapotocny •
-
small CRAM file to test
written by lien
-
Samtools Have Trouble "View" Bwa Mem Result, But Good With Bwa Sampe Result?
written by Chen Sun
-
Error In Sam To Bam Conversion
written by SK •
-
Tophat-Fusion Sam To Bam Error
written by Nicolas Rosewick
-
Truncated File, When Converting Sam To Bam
written by Jianfengmao •
Convert a CRAM file to a BAM file
Is cram-tools even still supported?
This suggests that there's a mismatch between the fasta file and what the CRAM file was created with. I'd double check that before proceeding.