1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
#include "ChimericSegment.h"
ChimericSegment::ChimericSegment(Parameters &Pin, Transcript &alignIn) : P(Pin), pCh(Pin.pCh), align(alignIn)
{
if ( (align.intronMotifs[1]==0 && align.intronMotifs[2]==0) || (align.intronMotifs[1]>0 && align.intronMotifs[2]>0)) {//strand is undefined
str=0;
} else if ( (align.Str==0) == (align.intronMotifs[1]>0)) {//strand the same as RNA.
//This assumes that the aligns have consistent strands, i.e. only intronMotifs[1]>0 OR intronMotifs[2]>0
str=1;
} else {//strand opposite to RNA
str=2;
};
roS=align.Str==0 ? align.exons[0][EX_R] : align.Lread - align.exons[align.nExons-1][EX_R] - align.exons[align.nExons-1][EX_L];
roE=align.Str==0 ? align.exons[align.nExons-1][EX_R] + align.exons[align.nExons-1][EX_L] - 1 : align.Lread - align.exons[0][EX_R] - 1;
if (roS>align.readLength[0]) roS--;
if (roE>align.readLength[0]) roE--;
};
bool ChimericSegment::segmentCheck()
{
bool segGood = true;
segGood = segGood && align.rLength >= pCh.segmentMin; //mapped length >= chim segmentMin
segGood = segGood && align.intronMotifs[0]==0; //no non-canonical unannotated juncions.
return segGood;
//this is already tested for each align with default --outFilterIntronStrands RemoveInconsistentStrands
//segGood = segGood && (align.intronMotifs[1]==0 || align.intronMotifs[2]==0); //consistent intron motifs.
//this is not requiered since seg2 is tested for length
// segGood = segGood && (align.exons[align.nExons-1][EX_R] + align.exons[align.nExons-1][EX_L] + P.pCh.segmentMin <= Lread
// || align.exons[0][EX_R] >= P.pCh.segmentMin); //uncovered by seg1 read length is <= segmentMin
};
|
