Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Nanopolish and pilon are calling wrong bases iltisanni Oxford Nanopore 0 06-26-2018 05:20 AM
Per base sequence quality and per base N content barthez95 Illumina/Solexa 1 05-22-2018 10:35 AM
Handling ambiguous nucleotides (IUPAC code) with MEGA software MLiorz Bioinformatics 0 04-06-2016 09:56 AM
Why is a seventh base needed with MiSeq 6-base indices? StopCodon Illumina/Solexa 4 03-06-2014 10:17 PM
denovo2 conversion pseudo-base for base jbdorr SOLiD 2 08-31-2012 10:28 AM

Thread Tools
Old 08-03-2018, 12:25 PM   #1
Location: Davis, CA

Join Date: Mar 2009
Posts: 82
Question Pilon ambiguous base handling?

Hi folks,

I'm using pilon to polishing a pacbio assembled genome with illumina data, and I'm confusing about how pilon treats single "N" characters for an ambigious base. There aren't gaps in the scaffolds, just single N characters representing an ambiguous base. For many of these there should be good support to correct this to an A, C, T, or G, but these aren't being touched by my current attempts with pilon. Ideas? Pilon is correcting other ambiguous bases (e.g. R, Y, K) to the correct base, but is ignoring all Ns.

The command I'm running is:
java -Xmx120g -jar ~/software/anaconda2/pkgs/pilon-1.22-1/share/pilon-1.22-1/pilon-1.22.jar \
    --genome ref.fasta \
    --frags aln.sorted.bam \
    --unpaired u.sorted.bam \
    --changes --vcf --tracks \
    --threads 16 \
    --fix bases,amb \
    --outdir pilon_02
ewilbanks is offline   Reply With Quote

genome assembly hybrid, illumina, pacbio correction, pilon, polish

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 06:28 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO