From: Jim Procter <jprocter@issues.jalview.org>
Date: Sat, 3 Mar 2018 11:17:34 +0000 (+0000)
Subject: JAL-2909 soft clip and handclap features (should use SO terms for ‘Clipped’)
X-Git-Url: http://source.jalview.org/gitweb/?p=jalview.git;a=commitdiff_plain;h=fe4666614489b0f7c98ba27c25a4432a5c10677d

JAL-2909 soft clip and handclap features (should use SO terms for ‘Clipped’)
---

diff --git a/src/jalview/datamodel/CigarParser.java b/src/jalview/datamodel/CigarParser.java
index cfbe842..a1124d4 100644
--- a/src/jalview/datamodel/CigarParser.java
+++ b/src/jalview/datamodel/CigarParser.java
@@ -170,8 +170,14 @@ public class CigarParser
       break;
     case S:
       // soft clipping - just skip this bit of the read
+      seq.addSequenceFeature(new SequenceFeature("SOFTCLIP", "",
+              nextPos + 1, nextPos + length, "bamfile"));
+      break;
     case H:
       // hard clipping - this stretch will not appear in the read
+      seq.addSequenceFeature(new SequenceFeature("HARDCLIP", "",
+              nextPos + 1, nextPos + length, "bamfile"));
+      break;
     default:
       break;
     }