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)	28 (0.630)	25/482	0.6	96.6%	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

ATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAACATTGGTGTAAACGTCAAACTGGTTAACCAGGAGTGGAAAACGTTCCTCGACACCCGTCACCAGGGTACTTTTGATGTGGCCCGTGCAGGCTGGTGTGCTGACTACAACGAACCAACTTCCTTCCTGAACACCATGCTTTCGAACAGCTCGATGAATACCGCGCATTATAAGAGCCCGGCCTTTGACAGCATTATGGCGGAAACGCTGAAAGTGACTGACGAGGCGCAGCGCAC  >  NC_000913/1302163‑1302647                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    aTTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGCTCTGCATAAAAAGCTGGCGATTGCTGCCTCTTAATTGTGaaaaccatt                                                                                                                                                                                                                                            >  2:152380‑M1/1‑242 (MQ=255)               aaaagGGCAGGGCAACAGGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAAGTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCCTCATTGTGaaaaacattggtgtaaacgtca                                                                                                                                                                                                                               <  2:265951‑M1/249‑23 (MQ=255)                                    cGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtgga                                                                                                                                                                                                          >  2:159126‑M1/1‑207 (MQ=255)                                     gCCTATGGATACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAACAACGAAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggcgtggac                                                                                                                                                                                                         >  2:397129‑M1/1‑206 (MQ=255)                                        tATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacg                                                                                                                                                                                                     >  2:9368‑M1/1‑203 (MQ=255)                                                   tCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacac                                                                                                                                                                                          >  1:313533‑M1/1‑192 (MQ=255)                                                              cTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGACACACGTACCGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagg                                                                                                                                                                               >  1:234986‑M1/1‑181 (MQ=255)                                                                   ggCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtac                                                                                                                                                                           <  1:424803‑M1/250‑75 (MQ=255)                                                                   ggCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCAAGGAAAAACGTAACGACGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaaccttggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtact                                                                                                                                                                          >  2:274428‑M1/1‑176 (MQ=255)                                                                              ggctcAGCCAGAATGGTTTGTCTGGAGCCAGGAAAAACGTAACGAAGAATCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgttg                                                                                                                                                                <  2:253413‑M1/248‑86 (MQ=255)                                                                                 ccaGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccg                                                                                                                                                            <  2:313533‑M1/250‑90 (MQ=255)                                                                                             gTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCCTAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagaggaaaacgttcctcgccacccgtcaccagggtccttttgatgaggcccgtgccggctgggg                                                                                                                                                >  2:239589‑M1/1‑150 (MQ=255)                                                                                                     ggAGCCGGGAAAAACGTAACGAAGAAGCGAAAATACTGCTGGCTGAAGCGGGATATACCGCAGACAAACCGTTGACCATCAACATGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgact                                                                                                                                        >  1:368005‑M1/1‑142 (MQ=255)                                                                                                       aaccaataaaaacttcactaataagcaaaaaaactgctggctgaagcGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGAGCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTTTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgcttactac                                                                                                                                      <  2:229448‑M1/229‑112 (MQ=255)                                                                                                         ccAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggcacttttgatgtggccagagcaggctggtgtgctgaataca                                                                                                                                     >  2:244915‑M1/1‑138 (MQ=255)                                                                                                          cAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGTTGAAGCGGGTTATACCGCAGACAAATCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaac                                                                                                                                   <  2:368005‑M1/251‑115 (MQ=255)                                                                                                                             aagCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctg                                                                                                                 >  1:404223‑M1/1‑118 (MQ=255)                                                                                                                                          gctggctgAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgcttt                                                                                                    <  1:108828‑M1/250‑146 (MQ=255)                                                                                                                                                tgAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaac                                                                                               >  1:196316‑M1/1‑99 (MQ=255)                                                                                                                                                        gttatagcgtacacccaCCGTTGACCATCTGCCTGTTGTATAAAAACTCCGATCTGCCTAAAAAGCCGGCGATTGCTGCCTCTTCATTG‑‑‑‑AGCAAACCCATGGTGTAAACGTCACACTGGTTAACCAGGAGTGGAACACGTTCATCGACAACCGTCACCCGGGTACTTTTGATGTGGCCCGTGCAGGCTGGCGTGCTGCCTACAACGCACAAAATTCCTTCCTGCACACCATCCTTTCGAACAGCCCGATg                                                                                  >  2:69383/17‑250 (MQ=255)                                                                                                                                                             aCCGCAGACAAACCGTTGACCAGCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaatac                                                                                 <  2:382265‑M1/250‑165 (MQ=255)                                                                                                                                                                    aCAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcaaacagctcgatgaataccgcgcatt                                                                         <  1:310305‑M1/251‑173 (MQ=255)                                                                                                                                                                          cGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattata                                                                      >  2:114250‑M1/1‑73 (MQ=255)                                                                                                                                                                                ccATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattataagagcccgg                                                             <  1:170138‑M1/251‑185 (MQ=255)                                                                                                                                                                                  aTCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattataagagcccggcc                                                           <  1:26440‑M1/251‑187 (MQ=255)                                                                                                                                                                                                           aTCTGCATAAAAAGCTGGCGATCGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacctgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattataagagaccggcctttgacagcattatggcggaaacgc                                  >  1:366403‑M1/1‑40 (MQ=255)                                                                                                                                                                                                                               aTTGCTGCCTCTTCATTGTGaaaaacattggtgtaaacgtcaaactggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattataagagcccggcctttgacagcattatggcggaaacgctgaaagtgactgacgaggcg              <  2:404223‑M1/251‑232 (MQ=255)                                                                                                                                                                                                                                           tCATTGTGaacaacattggtgtcaaagtcaaactggttaaccaggagtggaaaacgatcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctggtgtgctgactacaacgaaccaccttccttcctgaacaccatgctttcgcacagctcgatgaataccgcgcattataagagcccggcctttgacagcattatggaggaaacgctgaaggtgactgacgaggcgcagcgcaccgct  >  2:313031‑M1/1‑8 (MQ=255)                                                                                                                                                                                                                                           tCATTGTGaaaaacattggtgtaaacgtcaaacaggttaaccaggagtggaaaacgttcctcgacacccgtcaccagggtacttttgatgtggcccgtgcaggctgttgtgctgactacaacgaaccaacttccttcctgaacaccatgctttcgaacagctcgatgaataccgcgcattataagagcccggcctttgaaagcattatggcggaaacactgaaagtgacggacgaggcgcagcgcaccact  >  1:21174‑M1/1‑8 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    ATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAACATTGGTGTAAACGTCAAACTGGTTAACCAGGAGTGGAAAACGTTCCTCGACACCCGTCACCAGGGTACTTTTGATGTGGCCCGTGCAGGCTGGTGTGCTGACTACAACGAACCAACTTCCTTCCTGAACACCATGCTTTCGAACAGCTCGATGAATACCGCGCATTATAAGAGCCCGGCCTTTGACAGCATTATGGCGGAAACGCTGAAAGTGACTGACGAGGCGCAGCGCAC  >  NC_000913/1302163‑1302647

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