Reports sequences stored in one or more BAM files that cover one of more specific nucleotide positions in a reference.
🔙 To the main page of anvi’o programs and artifacts.
Reports sequences stored in a bam-file file that cover one of more specific nucleotide positions in a reference.
Assume you wish to recover reads stored in one or more BAM files, where the matching reads contain at least one nucleotide position that align to a nucleotide position
nucleotide_position_N in a contig
contig_name_X. In that case, the user would first generate a two column TAB-delmited file, for example
positions_for_linkmers.txt with no header line,
And run the program this way to recover the short reads from this way:
anvi-report-linkmers --contigs-and-positions positions_for_linkmers.txt \ -i SAMPLE_01.bam SAMPLE_02.bam SAMPLE_03.bam (...) \ -o linkmers.txt
The user can define multiple contigs in the input file, and one or more nucleotide positions for each one of them:
The resulting linkmers-txt would include all short reads that match any of these critera
--only-complete-links flag, the user can enforce whether only complete links should be reported where each reported short read must cover each nucleotide position for a given contig.
Please note that if the nucleotide positions chosen for a given contig are too distant from each other given the short read length, zero reads may satisfy the complete links criterion.
Having complete links, however, will enable oligotyping analyses on metagenomic reads through the anvi’o program anvi-oligotype-linkmers.
Edit this file to update this information.
Are you aware of resources that may help users better understand the utility of this program? Please feel free to edit this file on GitHub. If you are not sure how to do that, find the
__resources__ tag in this file to see an example.