Chimeric.out.junction file issues in segment 2 co-ordinates

18 views
Skip to first unread message

User_new

unread,
Oct 17, 2019, 1:46:41 PM10/17/19
to rna-star
I have single end reads and did 

STAR --runThreadN 8 --genomeDir $TRANS_DATA --readFilesIn $TRANS_DATA/5_S5_R1_001_trimmed.fq --outSAMtype BAM Unsorted --outFileNamePrefix chim --alignIntronMin 1 --chimSegmentMin 20 and get 

splash.PNG



When I look at the first segment say for example from line I get 841 - 891 , however for the second segment I always get 1971 - 1972 (meaning just a difference of 1). What am I doing wrong? I would like to get the start and end of both segments of the chimeras.

Alexander Dobin

unread,
Oct 19, 2019, 2:52:47 PM10/19/19
to rna-star
Hi @User_new

col 5 (e.g. 1971) is the last base of the intron of the chimeric junction, so the second segment starts at the next base (i.e. 1972).

Cheers
Alex
Reply all
Reply to author
Forward
0 new messages