Entering edit mode
4 months ago
SaltedPork
▴
170
Hi,
We sequenced a plasmid on MiSeq, 150bp reads. I have a reference, and an expected sequence of the insert. I want to know at what position does the insert begin in the plasmid.
So far I have generated contigs using plasmid SPAdes, and bactopia. Neither align to the insert.
Some questions:
- how can i turn these contigs into a whole genom for the plasmid?
- Should I be using plasmid SPAdes if we didnt sequence the bacteria, only the plasmid?
Thanks,
Does this mean you can't find the expected insert when you try to align?
See Identification of the sequence insertion site in the genome for inspiration.
You could also use
bbduk.sh
in filter mode to see if you are able to find reads that align to the expected insert. Guide forbbduk.sh
here: https://jgi.doe.gov/data-and-tools/software-tools/bbtools/bb-tools-user-guide/bbduk-guide/