Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
bios4biol
bioinfo-utils
Commits
92fbc971
Commit
92fbc971
authored
Mar 02, 2015
by
no_author
Browse files
No commit message
No commit message
parent
6d815d0c
Changes
1
Hide whitespace changes
Inline
Side-by-side
bin/splitbc.pl
View file @
92fbc971
...
...
@@ -541,75 +541,6 @@ sub match_sequences {
$seq_qualities_r2
=
substr
$seq_qualities_r2
,
$start2
,
$cut_len2
;
}
}
# trim the barcode
# if ($trim) {
# if ($barcodes_at_bol){
# $start=$barcodes_length+$adapt;
# $stop=$seq_len;
# } elsif ($barcodes_at_eol) {
# $start=0;
# $stop=$seq_len-$adapt;
# }
#
# #~ warn $seq_len," ",$barcodes_length;
## if ($adapt && $barcodes_at_bol){
## $start=$barcodes_length;
## $stop=$seq_len;
## }elsif ($adapt && $barcodes_at_eol){
## $start=0;
## $stop=$seq_len;
## }elsif ($barcodes_at_bol){
## $start=$barcodes_length+1;
## $stop=$seq_len;
## }elsif ($barcodes_at_eol){
## $start=0;
## $stop=$seq_len-1;
## }
# warn "position of fragment to keep : $start, $stop \n" if $debug;
# if ($paired==1 && $barcodes_at_eol){
# $seq_bases_r2 = substr $seq_bases_r2, $start, $stop;
# } else {
# $seq_bases = substr $seq_bases, $start, $stop;
# }
#
# if ($fastq_format) {
# if ($paired==1 && $barcodes_at_eol){
# $seq_qualities_r2 = substr $seq_qualities_r2, $start, $stop;
# } else {
# $seq_qualities = substr $seq_qualities, $start, $stop;
# }
# }
# }
#
# # trim both read
# if ($trim2){
# my $start;
# my $stop;
# my $stop2;
# if ($adapt && $barcodes_at_bol){
# $start=$barcodes_length;
# $stop=$seq_len;
# $stop2=$seq_len_r2;
# }elsif ($adapt && $barcodes_at_eol){
# $start=0;
# $stop=$seq_len;
# $stop2=$seq_len_r2;
# }elsif ($barcodes_at_bol){
# $start=$barcodes_length+1;
# $stop=$seq_len;
# $stop2=$seq_len_r2;
# }elsif ($barcodes_at_eol){
# $start=0;
# $stop=$seq_len-1;
# $stop2=$seq_len_r2-1;
# }
# $seq_bases = substr $seq_bases, $start, $stop;
# $seq_bases_r2 = substr $seq_bases_r2, $start, $stop2;
# if ($fastq_format){
# $seq_qualities = substr $seq_qualities, $start, $stop;
# $seq_qualities_r2 = substr $seq_qualities_r2, $start, $stop2;
# }
# }
}
if
(
(
!
defined
$best_barcode_ident
)
||
(
$best_barcode_mismatches_count
>
$allowed_mismatches
)){
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment