When aligning reads to a reference to generate a config for SNP/InDel detection, won't the alignment necessarily remove the reads that have the desired information?
I see this as a commonly-listed step before searching for mutations.
Other threads I have read here (and papers referred to in them) refer to software that does the alignments (I have used Bowtie and Ugene), but I clearly see the absence in the resulting contigs of known antibiotic resistance markers that are in the genomes (large insertions). My filtered Illumina sets are paired (which don't seem to work in the alignments), but the individual libraries provide about 10 fold coverage with most reads >200.
I see this as a commonly-listed step before searching for mutations.
Other threads I have read here (and papers referred to in them) refer to software that does the alignments (I have used Bowtie and Ugene), but I clearly see the absence in the resulting contigs of known antibiotic resistance markers that are in the genomes (large insertions). My filtered Illumina sets are paired (which don't seem to work in the alignments), but the individual libraries provide about 10 fold coverage with most reads >200.
Comment