Commit 6d621196 authored by TomKellyGenetics's avatar TomKellyGenetics
Browse files

correct processing hairpin barcodes for SCI-Seq

parent dddf7af5
Loading
Loading
Loading
Loading
+11 −2
Original line number Diff line number Diff line
@@ -2455,13 +2455,22 @@ else
                s/^.{18}//g
                }'  $convFile > ${crIN}/.temp
            mv ${crIN}/.temp $convFile
            #remove linker
            #remove linker (10 bp barcodes)
            sed -E '
                /^(.{10})CAGAGC/ {
                s/^(.{10})CAGAGC(.{18})/\1\2/g
                n
                n
                s/^(.{26})(.{10})(.{6})(.{18})/\2\4/g
                s/^(.{10})(.{6})(.{18})/\1\3/g
                }'  $convFile > ${crIN}/.temp
            mv ${crIN}/.temp $convFile
            #remove linker (9 bp barcodes)
            sed -E '
                /^(.{9})CAGAGC/ {
                s/^(.{9})CAGAGC(.{18})/E\1\2/g
                n
                n
                s/^(.{9})(.{6})(.{18})/F\1\3/g
                }' $convFile > ${crIN}/.temp
            mv ${crIN}/.temp $convFile
            #swap barcode and UMI