anvi-report-linkmers
Table of Contents
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.
Authors
Can consume
Can provide
Usage
Reports sequences stored in a bam-file file that cover one of more specific nucleotide positions in a reference.
Basic mode of operation
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,
| contig_name_X | nucleotide_position_N |
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:
| contig_name_X | nucleotide_position_01,nucleotide_position_02,nucleotide_position_03 |
| contig_name_Y | nucleotide_position_04 |
| contig_name_Z | nucleotide_position_05,nucleotide_position_06 |
The resulting linkmers-txt would include all short reads that match any of these critera
Complete or incomplete links?
Using the --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.
See this program in action
http://merenlab.org/2015/12/09/musings-over-commamox/
Edit this file to update this information.
Additional Resources
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.