HTseq code error repeatedly

Error occured when reading beginning of SAM/BAM file. ('SAM line does not contain at least 11 tab-delimited fields.', 'line 1 of file Sorted_KKUGCTM5_-VE_14-4_aligned.bam') [Exception type: ValueError, raised in _HTSeq.pyx:1276]

Warning: Malformed SAM line: MRNM != '*' although flag bit &0x0008 set Warning: Read 700463F:369:CB93CANXX:5:1103:15415:8299 claims to have an aligned mate which could not be found in an adjacent line. Error occured when processing SAM input (line 57 of file Sorted_KKUGCTM5_-VE_14-4_aligned1.sam): 'pair_alignments' needs a sequence of paired-end alignments [Exception type: ValueError, raised in __init__.py:603]

9 days ago by

michael.ante • 3.0k

Austria/Vienna

michael.ante • 3.0k wrote:

Hi, Please use the coding button to format your text. This will help others to understand your problem more easily.

To my understanding, htseq is assuming that your bam file is sorted a certain way. This assumption is violated by your file.

You can use the parameter -r pos if your reads are positionally sorted -r name otherwise.

Cheers, Michael

ADD COMMENT • link written 9 days ago by michael.ante • 3.0k

Thanks for pointing out how I should format my questions, I have done that now. With regards to the sorting by position or name, I believe the conversion of sam file to bam file and its subsequent sorting using the samtools code leads to sorting by name by default, so I am assuming mine is sorted by name because I continued with default parameters itself. Then in the Htseq website said that the default in Htseq is also by name so I didn't add that argument in my code. I still tried it now but the same error is coming up :(

ADD REPLY • link written 9 days ago by bnayer26 • 0

subsequent sorting using the samtools code leads to sorting by name

No. Default is co-ordinate sort. Here is samtools sort help.

Sort alignments by leftmost coordinates, or by read name when -n is used. An appropriate @HD-SO sort order header tag will be added or an existing one updated if necessary.

Can you show what samtools view -H your.bam | head -5 looks like?

ADD REPLY • link modified 8 days ago • written 8 days ago by genomax ♦ 62k

@HD VN:1.0  SO:queryname
@SQ SN:chr1 LN:248956422
@SQ SN:chr10    LN:133797422
@SQ SN:chr11    LN:135086622
@SQ SN:chr12    LN:133275309

this is what it looks like. is this the expected outcome? Thanks. I have added the <-n> in my code as well. The above output is what it looks like after I ran the code with <-n>

ADD REPLY • link modified 6 days ago • written 6 days ago by bnayer26 • 0

Looks like you file is sorted based on the sequence headers (names). So you will need to use -r name option as noted above by @michael.ante with htseq.

ADD REPLY • link written 6 days ago by genomax ♦ 62k

Also note that the default input format for htseq-count is sam, which is why the software yelled at you for not providing the right format. Rather than converting your data to .sam, you should tell htseq-count to expect .bam format.

ADD REPLY • link written 8 days ago by swbarnes2 ♦ 4.8k

What is the parameter I should introduce to specify that it is a bam file in the input? Shall I add -f bam in my HTseq code?

ADD REPLY • link modified 6 days ago • written 6 days ago by bnayer26 • 0

Similar posts • Search »