fix bugs while processing with cigar

build.sh

@@ -4,3 +4,4 @@ env GOOS=darwin GOARCH=amd64 go build -o REDItools2_darwin_amd64 .
 env GOOS=linux GOARCH=amd64 go build -o REDItools2_linux_amd64 .
 env GOOS=windows GOARCH=amd64 go build -o REDItools2_win_amd64.exe .
 
+# scp -P 1100 REDItools2_linux_amd64 [email protected]:/mnt/data8/zhangyiming/REDItools2

const.go

@@ -20,7 +20,7 @@ var conf config
 var region *Region
 
 const (
-	VERSION            = "0.0.3-beta"
+	VERSION            = "0.0.4-beta"
 	DefaultBaseQuality = 30
 )
 

fast.go

@@ -56,8 +56,10 @@ func workerFast(
 			}
 
 			start, index := 0, 0
+
 			for _, i := range record.Cigar {
 				if i.Type() == sam.CigarMatch {
+
 					for j := 1; j <= i.Len(); j++ {
 						at := index + j - 1
 
@@ -69,14 +71,23 @@ func workerFast(
 						genomic := start + record.Start
 
 						if _, ok := edits[genomic]; !ok {
+							// GL000220.1
+							if genomic - 1 >= len(chrRef) {
+								sugar.Error(record.Record)
+								sugar.Fatalf("%s: genomic - 1[%d] >= len(chrRef)[%d]", ref.Ref, genomic - 1, len(chrRef))
+							}
 							edits[genomic] = NewEditsInfo(ref.Ref, chrRef[genomic-1], genomic)
 						}
 
 						edits[genomic].AddReads(record, at)
 						start++
 					}
 					index += i.Len()
-				} else if i.Type() != sam.CigarDeletion || i.Type() != sam.CigarHardClipped || i.Type() != sam.CigarInsertion {
+				} else if i.Type() != sam.CigarDeletion &&
+					i.Type() != sam.CigarHardClipped &&
+					i.Type() != sam.CigarInsertion &&
+					i.Type() != sam.CigarSoftClipped {
+
 					start += i.Len()
 				}
 			}

slow.go

@@ -73,7 +73,10 @@ func worker(wg *sync.WaitGroup, refs chan *Region, w chan string, omopolymericPo
 						start++
 					}
 					index += i.Len()
-				} else if i.Type() != sam.CigarDeletion || i.Type() != sam.CigarHardClipped || i.Type() != sam.CigarInsertion {
+				} else if i.Type() != sam.CigarDeletion &&
+					i.Type() != sam.CigarHardClipped &&
+					i.Type() != sam.CigarInsertion &&
+					i.Type() != sam.CigarSoftClipped {
 					start += i.Len()
 				}
 			}