Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913	3,383,196	Δ3 bp	2.6%	coding (997‑999/1068 nt)	degS →	serine endoprotease, periplasmic

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3383199 =	118 (1.020)	3 (0.030)	3/536	NT	2.6%	coding (1000/1068 nt)	degS	serine endoprotease, periplasmic
	?	NC_000913	= 3383195	112 (0.980)					coding (996/1068 nt)	degS	serine endoprotease, periplasmic

CGACGGTCAGTACGCCCGTTTCTTCTCTCAACCGCTGCTGCTGGGTAAAAACGGCGTGGAAGAGCGTAAATCTATCGGTACCCTGAGCGCATTTGAACAGAACGCGCTGGAAGGTATGCTGGATACGCTGAAGAAAGATATCGCCCTGGGCGAAGAGTTCGTTAATAAGTAATTGATTAGCGGATAATAAAAAACCGGAGCACAGACTCCGGTTTTTTGTTTTGAGCGCACGACTTAATTGGTTGCCGGATATTCCTGAATGGTGACCTGCAGCGTTAACTGCTTATCATCACGCATCACTAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3383501‑3383199 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑actacAGGGATCACCGAACCAGGGCGAATTTCCGCCACCTGATCCATCGTCTCCAGAGCAGAGATGGCCGGTTTGTTATCCACCGAAATAATCAGATCGTTGACCTGAATACCCGCATTCGCCGCCGGGCCGTCAGGTGACACTTCATTAACCACGATCCCTTGCAGTTGATCTATACCACCGCCCTGCGCGTGCAGTGGTGCGATCTCACGTCCGCCGATACCAATGTAGCCGCGGATCACGCGACCATCGCGGATCAGCTTATCCATAATTTTGGTTGCTAACTGGAAAGGAATCGCAAAG  <  NC_000913/3383195‑3382898                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           CGACGGCTAGGACGGCCGTTTCTTCTCTCAACCGCTGCTGCTGGGTAAAAACGGCGTGGAAGAGCGTAAATCGATCGGTAGCCGGAGCGAATTTGAACAGAACGCGCTGGAAGGTATGCTGGATACGCTGAAGAAAGATATCGCCCTGGGCGAAGAGTTCGTTAATAAGTAATTGATTAGCGGATAATAAAAAACCGGAGCACAGACTCCGGTTTTTTGTTTTGAGCGCACGACTTAATTGGTTGCCGGATATTCCTGAATGGTGACCTGCAGCGTTAACTGCTTATCATCACGCATCAC                                                                                                                                                                                                                                                                                                               <  1:208597/300‑1                                                                                                                                               GCCCTGGGCGAAGAGTTAGTTAATAAGTAATTGATTAGCGGATAATAAAAAACCGGAGCACAGACTCCGGTTTTTTGTTTTGAGCGCACGACTTAATTGGTTGCCGGATATTCCTGAATGGTGACCTGCAGCGTTAACTGCTTATCATCACGCATCACTACAGGGATCACCGAACAAGGGCGAATTTCCGCCACCTGATCCATCGTCACCAGAGCAGAGATGGCCGGTTGGTTATCCACCGAAATAATCAGATCGTAGACCTGAATACCCGCATTCGCAGCCGGGCCGTCAGGGGACACTT                                                                                                                                                                >  1:691937/1‑301                                                                                                                                                                                           ATAAAAAACCGGAGCACAGACTCCGGTTTTTTGTTTTGAGCGCACGACTTAATTGGTTGCCGGATATTCCTGAATGGTGACCTGCAGCGTTAACTGCTTATCATCACGCATCACTACAGGGATCACCGAACCAGGGCGAATTTCCGCCACCTGATCCATCGTCTCCAGAGCGGAGATGGCAGGTTTATTATCCACCGAAATAATCAGATCGTTGACCTGAATACCCGCATTCGCCGCCGGGCCGTCAGGTGACACTTAATTAAACACACACCATTTCAGTTGATCTATACCACCGCCCTGa                                                                                                                    >  1:398305/1‑300                                                                                                                                                                                                                                   CGCACTACTTAATTGGGTGCACGAGATTCCTGAATGGTGCCCTCCCGCGTTAACTGCTTATCATCACGCATCACTACAGGGATCACCGAACCAGGGCGAATTTCCGCCACCTGATCCATCGTCTCCAGAGCAGAGAGGGCCGGTTTGTTATCCACCGAAATAATCAGATCGTTGACCTGAATACCCGCATTCGCCGCCGGGCCGTCAGGTGACACTTCATTAACCACGATCCCTTGCAGTTGATCTATACCACCGCCCTGCGCGTGCAGTGGTGCGATCTCACGTCCGCCGATACCAATGT                                                                            <  2:398300/301‑1                                                                                                                                                                                                                                                                                                   aaagaacCACTAAAGGGATCACCGAACCAGGGCGAATTTCCGCCACCTGATACAGCGTCTCCAGAGCAGAGATGGCCGGTTTGTTATCCACCGAAATAATCAGATCGTTGACCTGAATACCCGCATTCGCCGCCGGGCCGTCAGGTGACACTTCATTAACCACGATCCCTTGCAGTTGATCTATACCACCGCCCTGCGCGTGCAGTGGTGCGATCTCACGTCCGCCGATACCAATGTAGCCGCGGATCACGCGACCATCGCGGATCAGCTTATCCATAATTTTGGTTGCTAACTGGAAAGG            <  2:1039500/294‑1                                                                                                                                                                                                                                                                                                            CTACAGGGATCACCGAACCAGGGCGAATTTCCGCCACCTGATCCATCGTCTCCAGAGCAGAGATGGCCGGTTTGTTATCCACCGAAATAATCAGATCGTTGACCTGAATACCCGCATTCGCCGCCGGGCCGTCAGGTGACACTTCATTAACCACGATCCCTTGCAGTTGATCTATCCCACCGCCCTTCGCGTGCAGTGGTGCGATCTCCCGTCCGCCGATCCCAATGTAGCCGCGGATCACGCGCCCATCCCGGCTCAGCTTATCCATAATTTTGGTTGCTACCTGGAACGGACTCGCAAA   >  2:742781/1‑301                                                                                                                                                                                                                                                                                                             TACAGGGATCACCGAACCAGGGCGAATTTCCGCCACCTGATCCATCGTCTCCAGAGCAGAGATGGCCGGTTTGTTATCCACCGAAATAATCAGATCGTTGACCTGAATACCCGCATTCGCCGCCGGGCCGTCAGGTGACACTTCATTAACCACGATCCCTTGCAGTTGATCTATACCACCGCCCTGCGCGTGCAGTGGTGCGATCTACCGTCCGCCGATACCAATGTAGCCGCGGATCACGCGACCATCGCGGATCAGCTTAACCACACTTTTGTTTGCTAACTGGAAAGAAATCGCAAAa  >  2:375779/1‑300                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           CGACGGTCAGTACGCCCGTTTCTTCTCTCAACCGCTGCTGCTGGGTAAAAACGGCGTGGAAGAGCGTAAATCTATCGGTACCCTGAGCGCATTTGAACAGAACGCGCTGGAAGGTATGCTGGATACGCTGAAGAAAGATATCGCCCTGGGCGAAGAGTTCGTTAATAAGTAATTGATTAGCGGATAATAAAAAACCGGAGCACAGACTCCGGTTTTTTGTTTTGAGCGCACGACTTAATTGGTTGCCGGATATTCCTGAATGGTGACCTGCAGCGTTAACTGCTTATCATCACGCATCACTAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3383501‑3383199 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑actacAGGGATCACCGAACCAGGGCGAATTTCCGCCACCTGATCCATCGTCTCCAGAGCAGAGATGGCCGGTTTGTTATCCACCGAAATAATCAGATCGTTGACCTGAATACCCGCATTCGCCGCCGGGCCGTCAGGTGACACTTCATTAACCACGATCCCTTGCAGTTGATCTATACCACCGCCCTGCGCGTGCAGTGGTGCGATCTCACGTCCGCCGATACCAATGTAGCCGCGGATCACGCGACCATCGCGGATCAGCTTATCCATAATTTTGGTTGCTAACTGGAAAGGAATCGCAAAG  <  NC_000913/3383195‑3382898

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.