Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	1,303,073	Δ12 bp	coding (175‑186/921 nt)	oppB →	murein tripeptide ABC transporter/oligopeptide ABC transporter inner membrane subunit OppB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1303072	1 (0.020)	36 (0.570)	34/482	0.9	97.3%	coding (174/921 nt)	oppB	murein tripeptide ABC transporter/oligopeptide ABC transporter inner membrane subunit OppB
	?	NC_000913	1303085 =	1 (0.020)					coding (187/921 nt)	oppB	murein tripeptide ABC transporter/oligopeptide ABC transporter inner membrane subunit OppB

GGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGATCATGACACAGTATTTCAGCTACCTGAAACAACTGGCGCACGGTGATTTCGGTCCATCGTTTAAATATAAAGATTATTCGGTCAATGACTTGGTTGCATCCAGTTTTCCCGTTTCTGCCAAACTGGGAGCCGCAGCATTTTTCCTTGCGGTAATACTGGGTGTTAGTGCTGGCGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACC  >  NC_000913/1302826‑1303309                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              gggCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACTTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCAGCTTAATGATCCGtatt                                                                                                                                                                                                                                           <  2:38668‑M1/251‑5 (MQ=255)           gtCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacct                                                                                                                                                                                                                                 <  1:84956‑M1/251‑15 (MQ=255)                   ccACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaac                                                                                                                                                                                                                         <  2:20367‑M1/251‑23 (MQ=255)                                  ttaatcgCATTACAGAAGGCACAGGCCAGAAGGTAGGTCAATGTTAAAATTTATTCTACGTCGCTGGCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCATTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATAATCTTAATGATCAGtatttcagctacctcaaacaagtggcgcacgcacatt                                                                                                                                                                                                          >  2:441020‑M1/4‑214 (MQ=255)                                         cATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtcc                                                                                                                                                                                                   <  2:97063‑M1/251‑45 (MQ=255)                                                                               aaaaTTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTGCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttgg                                                                                                                                                             <  1:300105‑M1/251‑83 (MQ=255)                                                                                         cTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggggatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatcc                                                                                                                                                     >  1:490080‑M1/1‑159 (MQ=255)                                                                                                   gTCTGGAAGCGATTACGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtt                                                                                                                                         >  1:575184‑M1/1‑149 (MQ=255)                                                                                                              aTTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGTCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaact                                                                                                                              <  1:258364‑M1/251‑114 (MQ=255)                                                                                                                     cGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCAGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagcc                                                                                                                       >  2:614121‑M1/1‑131 (MQ=255)                                                                                                                             attCTTATTACTATTTCGTTCTTTATGATGCGCTTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCTCCAGAAGGGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGATCATGACACAGTATTTCAGCTACCTTAAACAACTGTCGCACGGTGATTGCGGTCCATCGTTTAAATATAAAGATTATTCGGTCAATGACTTGGTTGAATCCAGTTTTCCCGTTTCTGCCaacatggg                                                                                                               >  1:522212/1‑245 (MQ=255)                                                                                                                              ttCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgc                                                                                                                                    <  2:291489‑M1/229‑108 (MQ=255)                                                                                                                              ttCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgc                                                                                                                                    >  1:291489‑M1/1‑122 (MQ=255)                                                                                                                                     ttntagttcgttctTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGATCGTACTTTACCGCCAGAAGTTATGGCCAATATCTAATCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaacttgcgcacggtgatttcgggccttcgtttaaatataaagattattcggtcaattacttggttgcatccagttttcccgtttctgccaaactgggagctgcagcattcttcctt                                                                                                        <  2:599869‑M1/244‑136 (MQ=255)                                                                                                                                     taCTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttcctt                                                                                                        <  2:348126‑M1/250‑136 (MQ=255)                                                                                                                                      aCTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgc                                                                                                      >  1:135290‑M1/1‑114 (MQ=255)                                                                                                                                        tATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaattttaaagattattcgttcaatgacttggttgcatccatttttcccgtttctgccaacctgggagccgcagcatttttccttccgg                                                                                                    >  2:246132‑M1/1‑112 (MQ=255)                                                                                                                                         aTTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggt                                                                                                   <  1:3938‑M1/251‑141 (MQ=255)                                                                                                                                             cGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccttttctgccaaacttggagccgcagcatttttccttgcggtaata                                                                                               <  1:614121‑M1/251‑145 (MQ=255)                                                                                                                                                tCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactg                                                                                            >  1:291740‑M1/1‑104 (MQ=255)                                                                                                                                                   ttATGATGCGCCTCGCGCCGGGAAGGCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGCGATGGCCAATATCGAAGCGAGATATCATCTTAATGATCCGtatttcagctacctgaaacaacttgcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttgggtgcatccagttttctcgtttctgccaaactgggagccgcagaatttttcattgcgttaatactgggt                                                                                         >  2:349656‑M1/1‑101 (MQ=255)                                                                                                                                                            gcCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgct                                                                                <  1:246132‑M1/251‑160 (MQ=255)                                                                                                                                                                cgcgCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATAGCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccttttctgtcaatctgtgagccgcagcatttttccttgcggtaatactgggtgttagtgctgg                                                                              <  1:79431‑M1/249‑162 (MQ=255)                                                                                                                                                                 gcgcCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcgt                                                                           <  2:537904‑M1/251‑165 (MQ=255)                                                                                                                                                                       ggANGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctgcctgaaacaactggcgcgcggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccaggtttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagggctggcgtgattgc                                                                     <  2:592498‑M1/251‑171 (MQ=255)                                                                                                                                                                            cccTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcgtt                                                                          <  2:121413‑M1/241‑166 (MQ=255)                                                                                                                                                                            cccTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcgtt                                                                          >  1:121413‑M1/1‑76 (MQ=255)                                                                                                                                                                                gttaCCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagcttccggaaacaactggcgcacggtgatttcgggccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcattgtgccttgcggtaatactgggtgttagtgctggctttattgccgcatta                                                              <  2:385689‑M1/248‑178 (MQ=255)                                                                                                                                                                                   aCCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttctcgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcgttattgccgcattaaaaca                                                         <  1:539869‑M1/251‑183 (MQ=255)                                                                                                                                                                                     cGGCAAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcgatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcgtgattgccgcattaaaacaaa                                                       <  2:317534‑M1/251‑185 (MQ=255)                                                                                                                                                                                                    aCCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcgttattgccgcattaaaacaaaacaccaaatgggac                                         <  2:490080‑M1/250‑199 (MQ=255)                                                                                                                                                                                                          aGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggttcatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcgttattgccgcattaaaacaaaacaccaaatgggactataccg                                  <  2:342147‑M1/251‑206 (MQ=255)                                                                                                                                                                                                                                      tatCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtggtagtgctggcgttattgccgcattaaaacaaaacaccacatgggactataccgtgatggggctggcaatgaccgggggtgt      >  2:353449‑M1/1‑18 (MQ=255)                                                                                                                                                                                                                                       atcatcTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcgttattgccgcattaaaacaaaacaccaaatgggactataccgtgatggggctggcaatgaccggggttgtt     >  1:51990‑M1/1‑17 (MQ=255)                                                                                                                                                                                                                                         catcTTAATGCTCCGtctttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttcc                                                                                                          >  2:52637‑M1/1‑15 (MQ=255)                                                                                                                                                                                                                                         catcTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttcc                                                                                                          <  1:52637‑M1/148‑134 (MQ=255)                                                                                                                                                                                                                                              tAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcgttattgccgcattaaaacaaaacaccaaatgggactataccgtgatggggctggcaatgaccggggttgttatc  <  2:223863‑M1/247‑238 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              GGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGATCATGACACAGTATTTCAGCTACCTGAAACAACTGGCGCACGGTGATTTCGGTCCATCGTTTAAATATAAAGATTATTCGGTCAATGACTTGGTTGCATCCAGTTTTCCCGTTTCTGCCAAACTGGGAGCCGCAGCATTTTTCCTTGCGGTAATACTGGGTGTTAGTGCTGGCGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACC  >  NC_000913/1302826‑1303309

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