diff --git a/chemistry.yaml b/chemistry.yaml index bfe2a4d..5aa8b39 100644 --- a/chemistry.yaml +++ b/chemistry.yaml @@ -1,7 +1,7 @@ chromiumV3_illumina: filter_R1_length: 58 bc_cut: "" - STAR_R1: "[WHITELIST_V3,\"--soloUMIlen 12 --clip5pNbases 58 0 --soloCBstart 1 --soloCBlen 16 --soloUMIstart 17\"]" + STAR_R1: "[WHITELIST_V3,\"--soloUMIlen 12 --clip5pNbases 58 0 --soloCBstart 1 --soloCBlen 16 --soloUMIstart 17 --outFilterMultimapNmax 1 --outFilterMismatchNmax 999 --outFilterMismatchNoverReadLmax 0.2\"]" STAR_R2: "[WHITELIST_V3,\"--soloUMIlen 12\"]" STAR_paired: "--alignEndsProtrude 58 ConcordantPair" chromiumV3_element: @@ -17,31 +17,31 @@ chromiumV3_ultima: chromiumV2_illumina: filter_R1_length: 56 bc_cut: "" - STAR_R1: "[WHITELIST_V2,\"--soloUMIlen 10 --clip5pNbases 56 0 --soloCBstart 1 --soloCBlen 16 --soloUMIstart 17\"]" + STAR_R1: "[WHITELIST_V2,\"--soloUMIlen 10 --clip5pNbases 56 0 --soloCBstart 1 --soloCBlen 16 --soloUMIstart 17 --outFilterMultimapNmax 1 --outFilterMismatchNmax 999 --outFilterMismatchNoverReadLmax 0.2\"]" STAR_R2: "[WHITELIST_V2,\"--soloUMIlen 10\"]" STAR_paired: "--alignEndsProtrude 56 ConcordantPair" dropseq_illumina: filter_R1_length: 50 bc_cut: "" - STAR_R1: "None --soloUMIlen 8 --clip5pNbases 50 0 --soloCBstart 1 --soloCBlen 12 --soloUMIstart 13" + STAR_R1: "None --soloUMIlen 8 --clip5pNbases 50 0 --soloCBstart 1 --soloCBlen 12 --soloUMIstart 13 --outFilterMultimapNmax 1 --outFilterMismatchNmax 999 --outFilterMismatchNoverReadLmax 0.2" STAR_R2: "None --soloUMIlen 8 --soloCBstart 1 --soloCBlen 12 --soloUMIstart 13" STAR_paired: "--alignEndsProtrude 50 ConcordantPair" microwellseq_illumina: filter_R1_length: 54 bc_cut: "CGACTCACTACAGGG...TCGGTGACACGATCG" - STAR_R1: "None --soloUMIlen 6 --clip5pNbases 54 0 --soloCBstart 1 --soloCBlen 18 --soloUMIstart 19" + STAR_R1: "None --soloUMIlen 6 --clip5pNbases 54 0 --soloCBstart 1 --soloCBlen 18 --soloUMIstart 19 --outFilterMultimapNmax 1 --outFilterMismatchNmax 999 --outFilterMismatchNoverReadLmax 0.2" STAR_R2: "None --soloUMIlen 6 --soloCBstart 1 --soloCBlen 18 --soloUMIstart 19" STAR_paired: "--alignEndsProtrude 54 ConcordantPair" bd_illumina: filter_R1_length: 53 bc_cut: "ACTGGCCTGCGA...GGTAGCGGTGACA" - STAR_R1: "None --soloUMIlen 8 --clip5pNbases 53 0 --soloCBstart 1 --soloCBlen 27 --soloUMIstart 28" + STAR_R1: "None --soloUMIlen 8 --clip5pNbases 53 0 --soloCBstart 1 --soloCBlen 27 --soloUMIstart 28 --outFilterMultimapNmax 1 --outFilterMismatchNmax 999 --outFilterMismatchNoverReadLmax 0.2" STAR_R2: "None --soloUMIlen 8 --soloCBstart 1 --soloCBlen 27 --soloUMIstart 28" STAR_paired: "--alignEndsProtrude 53 ConcordantPair" indrop_illumina: filter_R1_length: 32 bc_cut: "" - STAR_R1: "None --soloUMIlen 6 --clip5pNbases 32 0 --soloCBstart 1 --soloCBlen 8 --soloUMIstart 9" + STAR_R1: "None --soloUMIlen 6 --clip5pNbases 32 0 --soloCBstart 1 --soloCBlen 8 --soloUMIstart 9 --outFilterMultimapNmax 1 --outFilterMismatchNmax 999 --outFilterMismatchNoverReadLmax 0.2" STAR_R2: "None --soloUMIlen 6 --soloCBstart 1 --soloCBlen 8 --soloUMIstart 9" STAR_paired: "--alignEndsProtrude 32 ConcordantPair" \ No newline at end of file diff --git a/rules/cutadapt_star.snake b/rules/cutadapt_star.snake index 2364f9a..7351bc1 100644 --- a/rules/cutadapt_star.snake +++ b/rules/cutadapt_star.snake @@ -153,9 +153,6 @@ rule starsolo_R1: --soloCBwhitelist {params.chemistry} \ --soloBarcodeMate 1 \ --soloStrand Reverse \ - --outFilterMultimapNmax 1 \ - --outFilterMismatchNmax 999 \ - --outFilterMismatchNoverReadLmax 0.2 \ --genomeDir {STAR_INDEX} \ --soloFeatures Gene \ --alignMatesGapMax 100000 \ @@ -240,9 +237,6 @@ rule starsolo_paired: {params.star_args_paired} \ --soloBarcodeMate 1 \ --soloStrand Reverse \ - --outFilterMultimapNmax 1 \ - --outFilterMismatchNmax 999 \ - --outFilterMismatchNoverReadLmax 0.2 \ --genomeDir {STAR_INDEX} \ --soloFeatures Gene \ --alignMatesGapMax 100000 \