Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	1,302,405	Δ4 bp	coding (1224‑1227/1632 nt)	oppA →	oligopeptide ABC transporter periplasmic binding protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1302404	1 (0.020)	31 (0.550)	30/490	0.7	96.9%	coding (1223/1632 nt)	oppA	oligopeptide ABC transporter periplasmic binding protein
	?	NC_000913	1302409 =	1 (0.020)					coding (1228/1632 nt)	oppA	oligopeptide ABC transporter periplasmic binding protein

GATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAACATTGGTGTAAACGTCAAACTGGTTAACCAGGAGTGGAAAACGTTCCTCGACACCCGTCACCAGGGTACTTTTGATGTGGCCCGTGCAGGCTGGTGTGCTGACTACAACGAACCAACTTCCTTCCTGAACACCATGCTTTCGAACAGCTCGATGAATACCGCGCATTATAAGAGCCCGGCCTTTGACAGCATTATGGCGGAAACGCTGAAAGTGACTGACGAGGCGCAGCGCA  >  NC_000913/1302157‑1302646                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   gATATCATTGTTAATCAAGAGAACGCGCAGGGCAACATGCCCGCCTATGGTTACACTACACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTACGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaa                                                                                                                                                                                                                                                  >  2:374628‑M1/1‑248 (MQ=255)  atatCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTGTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaa                                                                                                                                                                                                                                                <  1:75631‑M1/251‑5 (MQ=255)            taataaAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGCGaaaaacattggtgt                                                                                                                                                                                                                                      >  2:439834‑M1/1‑237 (MQ=255)              ataaAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTTTGaaaaacattggtgtaa                                                                                                                                                                                                                                    <  1:177245‑M1/251‑17 (MQ=255)                                                         tCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATACCACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacac                                                                                                                                                                                         >  1:526865‑M1/1‑192 (MQ=255)                                                               gTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtc                                                                                                                                                                                    >  1:421806‑M1/1‑186 (MQ=255)                                                                 ataCTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcacc                                                                                                                                                                                 >  1:158792‑M1/1‑184 (MQ=255)                                                                          gcgcAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtactt                                                                                                                                                                        <  1:167257‑M1/251‑77 (MQ=255)                                                                          gcgcAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgacaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtactt                                                                                                                                                                        <  1:337031‑M1/251‑77 (MQ=255)                                                                           cgcAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCCGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttacccaggagtggaaaacgttcctcgacacccgtcaccagggtacttt                                                                                                                                                                       >  2:488436‑M1/1‑174 (MQ=255)                                                                                        cAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggccc                                                                                                                                                            <  2:366621‑M1/249‑89 (MQ=255)                                                                                          gCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCGGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTTTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgc                                                                                                                                                        <  2:443778‑M1/251‑93 (MQ=255)                                                                                                           ggAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctga                                                                                                                                         <  1:498689‑M1/249‑108 (MQ=255)                                                                                                                   gAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaa                                                                                                                               <  2:99939‑M1/251‑118 (MQ=255)                                                                                                                                      cGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagaggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaac                                                                                                             >  1:400698‑M1/1‑115 (MQ=255)                                                                                                                                           aaaCTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatg                                                                                                       <  2:74177‑M1/251‑142 (MQ=255)                                                                                                                                                   ggctgAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcga                                                                                                >  1:261818‑M1/1‑102 (MQ=255)                                                                                                                                                           cGGGTTATACCGCAGACAAACCGTCGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcg                                                                                        >  2:173605‑M1/1‑94 (MQ=255)                                                                                                                                                                  taCCGCAGACAAAAAGTTGACCATCAACCTGTTGTATAACACCTCCGATATGCATAAAGAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaccgtcaaactggataaccaggagtggaaaacgttcctcgacacccgtcaccaggggaattttgatgcggcccgggcaggatggggtgctgactacaacgaaccaacttcattcctgaccaccatgctttcgaacagcgcgatcaatac                                                                                >  2:541804‑M1/1‑87 (MQ=255)                                                                                                                                                                                 gtcgaccATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattataagag                                                                 <  1:173605‑M1/248‑180 (MQ=255)                                                                                                                                                                                  ttGACCATCCACGTGTTGTATAACCCCTCCGTTCTGCCTAAAAAGCTGGCGATTGCTGCCTCTTCATTG‑‑‑‑AGAACAAAATTGGTGTAAACGTCACACTCGCTAACCTGGAGAGGAAAACGTTCCTCGACAACCGTCACCAGGGTACTTTTGATGTGGACCGTACAGGCTGATGTGCTGACTAAAACGAAACAACTTCCTTCCTGAACAACATGCTTTCCAACTGCTCGATGAATACCGCTCATTATAagcgc                                                            >  2:559518/1‑248 (MQ=255)                                                                                                                                                                                       cATCAACCTGGTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattataagagcccggc                                                           <  2:261818‑M1/251‑186 (MQ=255)                                                                                                                                                                                             ccTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattataagagcccggcctttg                                                      >  2:412407‑M1/1‑60 (MQ=255)                                                                                                                                                                                                                          aaaaGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattataagagcccggcctttgacagcattatggcggaaacgctgaaagtg                         >  1:196570‑M1/1‑31 (MQ=255)                                                                                                                                                                                                                                 ggCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgcacagctcgatgaataccgcgcattataagagcccggcctttgacagcattatggcggaaacgctgaaagtgactgacg                  >  2:216519‑M1/1‑24 (MQ=255)                                                                                                                                                                                                                                       tgctgcCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaatacagcgcattataagagcccggcctttgacagcattatggcggaaacgctgaaagtgactgacgaggcgcg           >  1:381709‑M1/1‑18 (MQ=255)                                                                                                                                                                                                                                           gcCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgttctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattataagagcccggcctttgacagcattatggcggaaacgctgaaagtgacttacgaggcgcagcg        <  2:199671‑M1/250‑237 (MQ=255)                                                                                                                                                                                                                                            ccTCTCCATTGTGaaaaacagtggtgtaaacgtcaaactggttaacaaggagtggaaaacgttcctcgacacccgtcaccagtgtactttggatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttactgaacaccatgctttcgaacagctcgatgaataccgcgcattataagagcccggcctgtgacagcattagggcggaaacgctgaaagtgactgacgaggctcagcgca      <  2:242500‑M1/251‑239 (MQ=255)                                                                                                                                                                                                                                             ctctTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattata                                                                     <  2:219993‑M1/187‑176 (MQ=255)                                                                                                                                                                                                                                             ctctTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattata                                                                     >  1:219993‑M1/1‑12 (MQ=255)                                                                                                                                                                                                                                                 tCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattataagagcccggcctttgacatcattatggcggaaacgctgaaagggactgacgaggcgcagcgcaca    <  1:216519‑M1/248‑241 (MQ=255)                                                                                                                                                                                                                                                 tCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattataagagcccggcctttgacagcattatggcggaaacgctgaaagtgactgacgaggcgcagcgcacagc  <  1:412407‑M1/250‑243 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   GATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAACATTGGTGTAAACGTCAAACTGGTTAACCAGGAGTGGAAAACGTTCCTCGACACCCGTCACCAGGGTACTTTTGATGTGGCCCGTGCAGGCTGGTGTGCTGACTACAACGAACCAACTTCCTTCCTGAACACCATGCTTTCGAACAGCTCGATGAATACCGCGCATTATAAGAGCCCGGCCTTTGACAGCATTATGGCGGAAACGCTGAAAGTGACTGACGAGGCGCAGCGCA  >  NC_000913/1302157‑1302646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ 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.