JAL-2909 soft clip and handclap features (should use SO terms for ‘Clipped’)
authorJim Procter <jprocter@issues.jalview.org>
Sat, 3 Mar 2018 11:17:34 +0000 (11:17 +0000)
committerJim Procter <jprocter@issues.jalview.org>
Sat, 3 Mar 2018 11:17:34 +0000 (11:17 +0000)
src/jalview/datamodel/CigarParser.java

index cfbe842..a1124d4 100644 (file)
@@ -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;
     }