SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Oxford Nanopore



Similar Threads
Thread Thread Starter Forum Replies Last Post
Calling Ambiguous Bases with Read Mapping roliwilhelm Illumina/Solexa 0 01-15-2018 01:40 PM
Nanopolish without reference genome jordi Bioinformatics 2 11-08-2017 04:27 AM
Nanopolish 1D anotherSAM Oxford Nanopore 4 11-03-2017 06:45 PM
NEWBIE: What is relation between SNP calling, INDEL calling and Genotype calling? arkilis Bioinformatics 4 10-01-2013 05:57 PM
Calling for help in using samtools mpileup: all Ns in reference bases ferris Bioinformatics 1 07-16-2012 05:10 AM

Reply
 
Thread Tools
Old 06-26-2018, 04:20 AM   #1
iltisanni
Member
 
Location: Mainz, Germany

Join Date: Mar 2017
Posts: 21
Default Nanopolish and pilon are calling wrong bases

Hi,

we sequenced and assembled a pichia genome with canu and afterwards polished it with nanopolish and with pilon+illumina reads.

We noticed that the nanopolish Assembly grows against the canu assembly and that it grows again after using pilon and illumina reads for polishing.

We then viewed our Assembly together with mapped illumina and nanopore reads in IGV and noticed that canu called some bases which nanopolish and pilon removed while polishing but those bases are correct since we can see them in the reads in IGV.

Why do nanopolish and pilon remove those bases although they are visible in both, illumina and nanopore reads...

commands we used for nanopolish and pilon:

Nanopolish:
nanopolish index -d PATH/TO/FAST5FILES -s PATH/TO/SEQUENCINGSUMMARY.TXT /PATH/TO/COMBINEDFASTQ

~/minimap2-2.10_x64-linux/minimap2 -d draft.mmi trimmed_contig.fasta
~/minimap2-2.10_x64-linux/minimap2 -ax map-ont -t 16 draft.mmi /PATH/TO/COMBINED.FASTQ | samtools sort [email protected] 16 -o reads_vs_trimmed_contig.sorted.bam
samtools index reads_vs_trimmed_contig.sorted.bam

python3 /usr/lib/nanopolish/nanopolish_makerange.py /PATH/TO/TRIMMED_CONTIG.FASTA | parallel --results nanopolish.results -P 3 nanopolish variants --consensus polished.{1}.fa -w {1} -r /PATH/TO/COMBINED.FASTQ -b /PATH/TO/READS_VS_TRIMMED_CONTIG.SORTED.BAM -g /PATH/TO/TRIMMED_CONTIG.FASTA -t 4 --min-candidate-frequency 0.1 --faster --fix-homopolymers

python3 /usr/lib/nanopolish/nanopolish_merge.py polished*.fa > polished_GENOME.fa

pilon:
java -jar pilon-1.22.jar \
--genome test.fasta \
--bam XXX.sorted.bam \
--fix "bases" \
--threads 16 \
--output test \
--outdir OUT
iltisanni is offline   Reply With Quote
Reply

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 05:49 AM.


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