Closed:How to find mismatches in an assembly
2
0
Entering edit mode
7.0 years ago
mary ▴ 10

Edited: I'm using Tablet software to view my sorted BAM file and trying to set a feature to show mismatch threshold at which bases on an alignment coverage track are marked. i. e. if the percent of nucleotide differs from the consensus sequence is greater than that % of reads, it would be marked as a feature.

for this, I need to create a BED or GFF3 file from the sorted BAM file, which has the place of every nucleotide in every contig where the nucleotides of mapped reads has variation greater than the defined percent.

for example, I want to have a bed file with a list of positions where 30% of reads has a nucleotide variation from what is assigned to the assembled contig.

sequencing error mismateches • 188 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 2452 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6