SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   samtools picard SamFormatConverter (http://seqanswers.com/forums/showthread.php?t=5503)

Bio.X2Y 06-14-2010 06:22 AM

samtools picard SamFormatConverter
 
Hi,

I'm trying to run SamFormatConverter to convert a BAM to a SAM, and I'm getting an error (see below).

The text suggests that samtools is trying to parse my BAM as a SAM. However, the file is definitely a BAM, has a '.bam' extension, and appears to have the necessary "magic number" at the top.

Any ideas if I might be doing something wrong?

Thanks!

Command:
java -jar <location>/SamFormatConverter.jar I=in.bam O=out.sam

Head in.bam:
BAM?@HD VN:1.0 GO:none SO:unsorted
@SQ SN:chrM LN:16571 AS:HG18 UR:/seq/references/Homo_sapiens_assembly18/v0/Homo_sapiens_assembly18.fasta M5:d2ed829b8a1628d16cbeee88e88e39eb SP:Homo sapiens

Exception:
Exception in thread "main" net.sf.samtools.SAMFormatException: Error parsing text SAM file. Not enough fields; Line 1
Line: BAM?@HD VN:1.0 GO:none SO:unsorted
at net.sf.samtools.SAMTextReader.reportFatalErrorParsingLine(SAMTextReader.java:169)
at net.sf.samtools.SAMTextReader.access$400(SAMTextReader.java:40)
at net.sf.samtools.SAMTextReader$RecordIterator.parseLine(SAMTextReader.java:261)
at net.sf.samtools.SAMTextReader$RecordIterator.next(SAMTextReader.java:224)
at net.sf.samtools.SAMTextReader$RecordIterator.next(SAMTextReader.java:196)
at net.sf.picard.sam.SamFormatConverter.doWork(SamFormatConverter.java:64)
at net.sf.picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:150)
at net.sf.picard.sam.SamFormatConverter.main(SamFormatConverter.java:72)

bh1 04-26-2011 11:16 AM

Hi,
Did you figure a solution for this problem? I'm getting exactly the same error message.
thanks.

Richard Finney 04-26-2011 01:19 PM

Validation_stringency=lenient
 
Try VALIDATION_STRINGENCY=LENIENT in the java invocation.

Example:

java-Xmx67T -jar wherever/SamTowhatevernameis.jar \
VALIDATION_STRINGENCY=LENIENT \ INPUT=uncle.sam O=nephew.sam


That might shut up the error message and get Picard to accept the input.

bw. 07-08-2013 07:46 AM

In this case the issue is probably the
"BAM?" on the 1st line in front of "@HD VN:1.0 GO:none SO:unsorted"
I would try deleting it.


All times are GMT -8. The time now is 12:20 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.