Hi,
I am planning to use the --very-sensitive option for aligning my ATAC-seq reads. However, after reading the manual section below, I am unclear as to whether this option should be combined with --no-1mm-upfront or not?
"By default, Bowtie 2 will attempt to find either an exact or a 1-mismatch end-to-end alignment for the read before trying the multiseed heuristic. Such alignments can be found very quickly, and many short read alignments have exact or near-exact end-to-end alignments. However, this can lead to unexpected alignments when the user also sets options governing the multiseed heuristic, like -L and -N."
Many thanks,
Lucy
•
written
by
Lucy •
1 answer
--very-sensitive doesn't change the permitted number of mismatches to 0, so you don't need to specify --no-1mm-upfront.
Log in to answer this question.
-
Goodness-of-fit test - observed vs. expected (some expected values < 5)
written by Lucy •
-
Chip-seq H3k27 Bowtie2 Alignement
written by zeeshan.fazal1986 •
-
Should batch correction be performed before or after feature selection in scRNA-seq
written by Lucy •
-
Aligning RNA-seq and ATAC-seq gene annotations
written by Lucy •
-
HOMER de novo motif enrichment
written by Lucy •
-
ATAC-seq transcription factor binding prediction
written by Lucy •
-
CENTIPEDE - values from fitCentipede
written by Lucy •
-
featureCounts allowMultiOverlap - fractional or not?
written by Lucy •
-
Determine peak summits after IDR analysis
written by Lucy •
-
When using MACS2 with -f BAM, it is possible to have a shift other than 0?
written by Lucy •
