This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Problem in converting SAM to BAM file by DROP SEQ
Exception in thread "main" htsjdk.samtools.SAMFormatException: Error parsing SAM header. @RG line missing SM tag. Line:
@RG ID:staralign1; ; Line number 10
    at htsjdk.samtools.SAMTextHeaderCodec.reportErrorParsingLine(SAMTextHeaderCodec.java:258)
    at htsjdk.samtools.SAMTextHeaderCodec.access$200(SAMTextHeaderCodec.java:46)
    at htsjdk.samtools.SAMTextHeaderCodec$ParsedHeaderLine.requireTag(SAMTextHeaderCodec.java:358)
    at htsjdk.samtools.SAMTextHeaderCodec.parseRGLine(SAMTextHeaderCodec.java:168)
    at htsjdk.samtools.SAMTextHeaderCodec.decode(SAMTextHeaderCodec.java:110)
    at htsjdk.samtools.SAMTextReader.readHeader(SAMTextReader.java:216)
    at htsjdk.samtools.SAMTextReader.<init>(SAMTextReader.java:63)
    at htsjdk.samtools.SAMTextReader.<init>(SAMTextReader.java:73)
    at htsjdk.samtools.SamReaderFactory$SamReaderFactoryImpl.open(SamReaderFactory.java:444)
    at htsjdk.samtools.SamReaderFactory$SamReaderFactoryImpl.open(SamReaderFactory.java:209)
    at picard.sam.SortSam.doWork(SortSam.java:154)
    at picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:308)
    at picard.cmdline.PicardCommandLine.instanceMain(PicardCommandLine.java:103)
    at picard.cmdline.PicardCommandLine.main(PicardCommandLine.java:113)
samtobam

Hi VANSHIKA,

welcome.

it might be good to add some more info on your "question". Now it only shows an error message which, without context, is rather meaningless.

What exactly are you stuck on? what are you trying to achieve? The "error" already gives you an indication what might go wrong, did you look into that?

It looks like RG lines are present, it just needs sample name in the line.

0 answers

No answers yet.

Log in to answer this question.