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)	29 (0.720)	25/488	0.6	96.7%	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

GCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGATCATGACACAGTATTTCAGCTACCTGAAACAACTGGCGCACGGTGATTTCGGTCCATCGTTTAAATATAAAGATTATTCGGTCAATGACTTGGTTGCATCCAGTTTTCCCGTTTCTGCCAAACTGGGAGCCGCAGCATTTTTCCTTGCGGTAATACTGGGTGTTAGTGCTGGCGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTT  >  NC_000913/1302841‑1303318                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       gCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaac                                                                                                                                                                                                                                     <  1:120847‑M1/251‑20 (MQ=255)                                 gAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagcttcctgaaacaactggcgcacggtgatttcggtccatcgttt                                                                                                                                                                                                     <  1:334223‑M1/251‑52 (MQ=255)                                    ggCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaa                                                                                                                                                                                                   >  2:99008‑M1/1‑197 (MQ=255)                                     gCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaat                                                                                                                                                                                                 <  1:82921‑M1/251‑56 (MQ=255)                                         aGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataa                                                                                                                                                                                             <  2:12916‑M1/251‑60 (MQ=255)                                          ggCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaa                                                                                                                                                                                            <  1:362660‑M1/251‑61 (MQ=255)                                                 aGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtat                                                                                                                                                                                                                                                     >  2:8560‑M1/1‑184 (MQ=255)                                                 aGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtat                                                                                                                                                                                                                                                     <  1:8560‑M1/187‑4 (MQ=255)                                                          aaTGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatg                                                                                                                                                                            >  1:84993‑M1/1‑175 (MQ=255)                                                                    tttATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGACCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGCAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgc                                                                                                                                                                  >  1:31531‑M1/1‑165 (MQ=255)                                                                                 gCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATTATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttccc                                                                                                                                                     <  2:223250‑M1/251‑100 (MQ=255)                                                                                    gTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtt                                                                                                                                                  >  1:217665‑M1/1‑149 (MQ=255)                                                                                         gAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgc                                                                                                                                             >  2:364985‑M1/1‑144 (MQ=255)                                                                                             cGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaa                                                                                                                                          >  1:300727‑M1/1‑140 (MQ=255)                                                                                                ttacgacGCTATTTATTCTTATTCCTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAATCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGGGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggcgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaacctg                                                                                                                                      >  2:18375‑M1/4‑137 (MQ=255)                                                                                                   cGACGCTATTTCTTCTTATAACCATTTCGTTCTTTATGATCCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTCTACCGCCAGAAGTGATGGCCAATATCGAAGCGAACTATCATCTTAATGATCCGtctttcagctacctgaaacaactgccgcacggtgttttcggtccatcgtttaactatacagattattcggtcaatgacttggttgcatccagttttcccgtttttgccaaactggg                                                                                                                                    >  2:42387‑M1/1‑134 (MQ=255)                                                                                                    gACGCTATTTATTCTTATTACTATTTCGCTCTTTATGATGCGCCTCGCGCCGGGAATCCCTTTTACCGGCGAACGTACTTTACCGCCAGACGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGgatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataagaattattcggtcaatgacttggttgaatccagttttaccgtttctgccaaactggg                                                                                                                                    >  2:105292‑M1/1‑133 (MQ=255)                                                                                                               ttCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcattt                                                                                                                       >  2:310964‑M1/1‑122 (MQ=255)                                                                                                                    attaCTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttcc                                                                                                                   <  2:300727‑M1/250‑134 (MQ=255)                                                                                                                             tCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaa                                                                                                          <  1:99008‑M1/250‑143 (MQ=255)                                                                                                                                          ggggcctcGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGgatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaggtgggagccgcagcatgtttcctttcggtaatactgggtgttagt                                                                                            <  2:179286‑M1/248‑157 (MQ=255)                                                                                                                                                 cgcgCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcg                                                                                     <  1:104027‑M1/251‑164 (MQ=255)                                                                                                                                                    gcCGGTAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAATTGAGGGCCAATATCGAAGCGAAATATCATCTTAATGATCC‑‑‑‑‑‑‑‑‑‑‑‑GTATTTCGGCTTCCTGAAACAACTGGCGCACGGTGATGTCGGTGCATCGTTTAAATATAAAGAGTATGCTATCAATGATTTGGTTGCGTCCGGTTTTCCCGTTTCGGCCAAACGTGGTGCCGCATCATTGTGCCTAGCTGTGATAGTTGGGGTTAGTACTGGCGTTa                                                                      <  2:232792/251‑1 (MQ=255)                                                                                                                                                      cGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcgttatt                                                                                >  1:146857‑M1/1‑83 (MQ=255)                                                                                                                                                             cccTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcgttattgccgcat                                                                         <  2:152256‑M1/251‑176 (MQ=255)                                                                                                                                                                                                                  cGAAATATCATCTTAATGATCCGtatttcatctacctgaaacaactggcgcacggtgctttcggtccatcggttgaatataaagcttattcggtcaatgacttggatgcatccagttttcgcgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcgttattgccgcattaaaacaaaacaccaaatgggactataccgtgatggggctggcaatgaccggg                    >  1:391239‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                  cGAAATATCATCTTAATGATCCGtatttcaggtacctgaaacaactgacgcacggtggtttcggtccattgtttaaatataaagattattcgttcaatgacttgtttgcatccagtttttccgtttctgccatactgagagccgcagcatttttcctggcggttatactgggtgttagtgcttgcgttattgccgcattaaaacaaaacaccaaatgggactataccgtgattgtgttggcaatgatcggg                    >  1:231033‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                  cGAAATATCATCTTAATGATCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggagttagtgctggcgttattgccgcattaaaacaaaacaccaaatgggactataccgtgatggggctggcaatgaccg                      >  1:70494‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                           ggcttAATGATCCGtatttcagctaactgaaacaactgccgcacggttatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttccagtttctgccgaactaggagccgcagcatttttccttgcggtcatattgggtgttagtgatggcgttattgccgcattaaaacaaggcaccaagtgggactttaccttgattgggctggcaatgagcgggtttgtcatt           >  2:178225‑M1/3‑14 (MQ=255)                                                                                                                                                                                                                                    aTCCGtatttcagctacctgaaacaactggcgcacggtgatttcggtccatcgtttaaatataaagattattcggtcaatgacttggttgcatccagttttcccgtttctgccaaactgggagccgcagcatttttccttgcggtaatactgggtgttagtgctggcgttattgccgcattaaaacaaaacaccaaatgggactataccgtgatggggctggcaatgaccggggttgttatccccagtttt  <  2:99530‑M1/251‑247 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       GCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGATGCGCCTCGCGCCGGGAAGCCCTTTTACCGGCGAACGTACTTTACCGCCAGAAGTGATGGCCAATATCGAAGCGAAATATCATCTTAATGATCCGATCATGACACAGTATTTCAGCTACCTGAAACAACTGGCGCACGGTGATTTCGGTCCATCGTTTAAATATAAAGATTATTCGGTCAATGACTTGGTTGCATCCAGTTTTCCCGTTTCTGCCAAACTGGGAGCCGCAGCATTTTTCCTTGCGGTAATACTGGGTGTTAGTGCTGGCGTTATTGCCGCATTAAAACAAAACACCAAATGGGACTATACCGTGATGGGGCTGGCAATGACCGGGGTTGTT  >  NC_000913/1302841‑1303318

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