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

CTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGG  >  NC_000913/1299353‑1299718
                                                                                                                                                 |                                                                                                                                                                                                                            
ctcAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctcta                                                                                                                                            >  2:114975‑M1/1‑146 (MQ=255)
                         tttCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgt                                                                                                                   <  1:114975‑M1/228‑108 (MQ=255)
                          ttCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAAGAAAGTGCTAAataataatcaattgttaaattattgt                                                                                                                                                                                                    >  2:348169‑M1/1‑120 (MQ=255)
                                     cTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggat                                                                                                       >  1:66104‑M1/1‑109 (MQ=255)
                                                                    atcatAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaa                                                                        >  1:112799‑M1/1‑78 (MQ=255)
                                                                                ggCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattttgatacggtcaataaaagagaattgctt                                                            >  2:504488‑M1/1‑66 (MQ=255)
                                                                                ggCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatgga                                                                                                        >  1:298750‑M1/1‑66 (MQ=255)
                                                                                 gCTTTCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggat                                                                                                       <  2:298750‑M1/184‑120 (MQ=255)
                                                                                                atTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcattaaaagagaattgcttaacaattttcaaaaat                                            >  2:98187‑M1/1‑50 (MQ=255)
                                                                                                              attaATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagaaa                              >  2:381713‑M1/1‑36 (MQ=255)
                                                                                                                    tCCACTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccg                        >  1:290889‑M1/1‑30 (MQ=255)
                                                                                                                       aCTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttaaccaatctctatttgattgaattgaaagatgtttgttacggcatggatgcaagctatagattctgttacggtcactaaaaaagaattgcttaacaattttgcaaaaagaattggcaagtaaaaaacgcat                     >  2:283943‑M1/1‑27 (MQ=255)
                                                                                                                        cTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatt                    >  2:384969‑M1/1‑26 (MQ=255)
                                                                                                                        cTGCCTTTGGTAGAGGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgactagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatcgattctgatatggtcaataaaagagatatgcttaccaattttgcaaaatgtattggcgagtaagaaccgcctt                    >  1:365404‑M1/1‑26 (MQ=255)
                                                                                                                                     agGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttatccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcactaaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaaaccatttggtactttccg        >  1:268108‑M1/1‑13 (MQ=255)
                                                                                                                                     agGAAAGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatttggtactttccgg       >  1:612527‑M1/1‑13 (MQ=255)
                                                                                                                                          aGTGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattggcgagtaagaaccgcatttgggactttccgggcaac  >  1:307038‑M1/1‑8 (MQ=255)
                                                                                                                                            tGCTAAataataatcaattgttaaattattgtgcatttcactactggaactgtaatcagaaaagatagacatgcttagccaatctctatttgattgaattgaaagatgtttgttaaggcatggatgcaagctatagattctgatacggtcaataaaagagaattgcttaacaattttgcaaaatgtattgtcgagtaagaaccgcatttggtactttccgggcaac  >  1:191381‑M1/1‑6 (MQ=255)
                                                                                                                                                 |                                                                                                                                                                                                                            
CTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTTTCGCCACATCATAATCCTGGGCTTGCTGAAGAATAATTGAAATGATATTATTAATTCCACTGCCTTTGGTAGAGGAAAGTGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGG  >  NC_000913/1299353‑1299718

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.