Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913 1,299,499 Δ1,199 bp 100% intergenic (+254/‑485) ychE → / → oppA UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1299499–1300693 1300697 5–1199 18 [0] [0] 18 ychE/oppA UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12994980 (0.000)16 (0.330) 13/406 NT 100% intergenic (+253/‑1684) ychE/oppA UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein
?NC_000913 1300698 = 0 (0.000)intergenic (+1453/‑484) ychE/oppA UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

TGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAAC  >  NC_000913/1300552‑1300917
                                                                                                                                                  |                                                                                                                                                                                                                           
ctcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTa                                                                                                                                            >  2:114975‑M2/147‑228 (MQ=255)
                         tttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttgt                                                                                                                   <  1:114975‑M2/107‑1 (MQ=255)
                          ttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaagaaagtgctaaATAATAATCAATTGTTAAATTATtgt                                                                                                                                                                                                    >  2:348169‑M2/121‑146 (MQ=255)
                                     cttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGAt                                                                                                       >  1:66104‑M2/110‑228 (MQ=255)
                                                                    atcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATaaa                                                                        >  1:112799‑M2/79‑228 (MQ=255)
                                                                                ggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGa                                                                                                        >  1:298750‑M2/67‑184 (MQ=255)
                                                                                ggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTTTGATACGGTCAATAAAAGAGAATTGCtt                                                            >  2:504488‑M2/67‑228 (MQ=255)
                                                                                 gctttctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGAt                                                                                                       <  2:298750‑M2/119‑1 (MQ=255)
                                                                                                attgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCATTAAAAGAGAATTGCTTAACAATTTTCAAAAAt                                            >  2:98187‑M2/51‑228 (MQ=255)
                                                                                                              attaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGaaa                              >  2:381713‑M2/37‑228 (MQ=255)
                                                                                                                    tccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCg                        >  1:290889‑M2/31‑228 (MQ=255)
                                                                                                                       actgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAACCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTACGGCATGGATGCAAGCTATAGATTCTGTTACGGTCACTAAAAAAGAATTGCTTAACAATTTTGCAAAAAGAATTGGCAAGTAAAAAACGCAt                     >  2:283943‑M2/28‑228 (MQ=255)
                                                                                                                        ctgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCAtt                    >  2:384969‑M2/27‑228 (MQ=255)
                                                                                                                        ctgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGACTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATCGATTCTGATATGGTCAATAAAAGAGATATGCTTACCAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCCtt                    >  1:365404‑M2/27‑228 (MQ=255)
                                                                                                                                     aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTATCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCACTAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAAAACCATTTGGTACTTTCCg        >  1:268108‑M2/14‑227 (MQ=255)
                                                                                                                                     aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCgg       >  1:612527‑M2/14‑228 (MQ=255)
                                                                                                                                          agtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGGACTTTCCGGGCAAc  >  1:307038‑M2/9‑228 (MQ=255)
                                                                                                                                            tgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGTCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAAc  >  1:191381‑M2/7‑226 (MQ=255)
                                                                                                                                                  |                                                                                                                                                                                                                           
TGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAAAAGAGAATTGCTTAACAATTTTGCAAAATGTATTGGCGAGTAAGAACCGCATTTGGTACTTTCCGGGCAAC  >  NC_000913/1300552‑1300917

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