Reports sequences stored in one or more BAM files that cover one of more specific nucleotide positions in a reference.
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,
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
Complete or incomplete links?
--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.
See this program in action
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.