[med-svn] [Git][med-team/rna-star][upstream] 2 commits: New upstream version 2.6.1c+dfsg

[Steffen Möller]

Steffen Möller pushed to branch upstream at Debian Med / rna-star


Commits:
ced0e9fb by Steffen Moeller at 2018-10-23T11:28:21Z
New upstream version 2.6.1c+dfsg
- - - - -
830c815a by Steffen Moeller at 2018-10-23T11:31:56Z
New upstream version 2.6.1c+dfsg
- - - - -


5 changed files:

- CHANGES.md
- source/ReadAlign_chimericDetectionPEmerged.cpp
- source/SequenceFuns.cpp
- source/VERSION
- source/readLoad.cpp


Changes:

=====================================
CHANGES.md
=====================================
@@ -1,3 +1,8 @@
+STAR 2.6.1c 2018/10/17
+======================
+
+* Enforced the consistent choice of supplementary chimeric alignments for overlapping mates.
+
 STAR 2.6.1b 2018/09/06
 ======================
 


=====================================
source/ReadAlign_chimericDetectionPEmerged.cpp
=====================================
@@ -39,7 +39,7 @@ void ReadAlign::chimericDetectionPEmerged(ReadAlign &seRA) {
                 };
             };
             for (uint jj=0; jj<2; jj++) {
-                if (segLen[ii][jj]<segLmin) {
+                if (segLen[ii][jj]<segLmin || (segLen[ii][jj]==segLmin && trChim[ii].exons[0][EX_G]>trChim[ii-1].exons[0][EX_G])) {
                     segLmin=segLen[ii][jj];
                     i1=ii;//trChim of the shortest segment length
                     i2=jj;//mate of the shortest segment length


=====================================
source/SequenceFuns.cpp
=====================================
@@ -139,7 +139,7 @@ void convertNucleotidesToNumbers(const char* R0, char* R1, const uint Lread) {//
                         case (84): case(116): 
                             R1[jj]=char(3);break;//T
                         default:                              
-                            R1[jj]=char(9);//anything else
+                            R1[jj]=char(4);//anything else is converted to N
                     };
                 };
 };
@@ -160,7 +160,7 @@ uint convertNucleotidesToNumbersRemoveControls(const char* R0, char* R1, const u
                 if (int(R0[jj]) < 32) {//control characters are skipped
                     continue;
                 } else {//all non-control non-ACGT characters are convreted to N
-                    R1[jj]=char(9);//anything else
+                    R1[jj]=char(4);//anything else
                 };
         };
         ++iR1;


=====================================
source/VERSION
=====================================
@@ -1 +1 @@
-#define STAR_VERSION "STAR_2.6.1b"
+#define STAR_VERSION "STAR_2.6.1c"


=====================================
source/readLoad.cpp
=====================================
@@ -84,7 +84,7 @@ int readLoad(istream& readInStream, Parameters& P, uint iMate, uint& Lread, uint
     } else {
         Lread=0;
     };
-    convertNucleotidesToNumbers(Seq+P.clip5pNbases[iMate],SeqNum,Lread);
+    convertNucleotidesToNumbersRemoveControls(Seq+P.clip5pNbases[iMate],SeqNum,Lread);
 
     //clip the adapter
     if (P.clip3pAdapterSeq.at(iMate).length()>0) {



View it on GitLab: https://salsa.debian.org/med-team/rna-star/compare/f926e182047e5dcf97be248fa0d30b9850a0ca1a...830c815a780f79edcea5d2d2eeaef6f29228f75a

-- 
View it on GitLab: https://salsa.debian.org/med-team/rna-star/compare/f926e182047e5dcf97be248fa0d30b9850a0ca1a...830c815a780f79edcea5d2d2eeaef6f29228f75a
You're receiving this email because of your account on salsa.debian.org.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://alioth-lists.debian.net/pipermail/debian-med-commit/attachments/20181023/13463168/attachment-0001.html>