How to find eversions with Pysam code?

[Martino]

I would like to write code to find eversion.

To find inversion, the code for this if statement is as follows:

elif read.is_reverse == read.mate_is_reverse:

But I cannot figure out what the code is for finding eversion.

Does somebode know this code?

I tried:

elif read.is_reverse == True and read.mate_is_reverse == False:

and

elif read.is_reverse == False and read.mate_is_reverse == True:

But both won't work. How can that be?

Can somebody help me?

[maubp]

Are you checking both the read and its mate are mapped?

P.S. From a style point of view it is more common to do "if statement:" and "if not statement:" rather than "if statement==True:" and "if statement=False:"

P.P.S. I think you mean "inversions" not "eversions". Places where the sequence has (somehow) been inverted (turned round).