Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,299,499 Δ1,199 bp 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–1300697 1300697 1–1199 25 [0] [0] 25 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)25 (0.830) 20/226 0.4 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

CTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGAT  >  NC_000913/1300585‑1300829
                                                                                                                 |                                                                                                                                   
actgcttctcctcaccataatgttattttcgccacatcataatcctgggcttgatgaagaataatagaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgt                                                                                                            >  2:519319‑M2/114‑139 (MQ=255)
 ctgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgtg                                                                                                           <  2:372178‑M2/27‑1 (MQ=255)
 ctgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgtg                                                                                                           <  1:376821‑M2/27‑1 (MQ=255)
     ttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAtt                                                                                                       >  1:309735‑M2/109‑139 (MQ=255)
                   atgttattttcgccacatcataatccttggcttgctgaagaataattgaaatgatattaataattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACt                                                                                         <  2:309735‑M2/45‑1 (MQ=255)
                         ttttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATc                                                                                   >  2:279381‑M2/89‑139 (MQ=255)
                            tcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggacagtgctaaATAATAATCAATTGTTAAATTATTCTGCATTTCAATACTGGAACTGTACTCAGa                                                                                >  2:561158‑M2/86‑139 (MQ=255)
                             cgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGaa                                                                               >  2:472626‑M2/85‑139 (MQ=255)
                                          tcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGc                                                                  <  1:279381‑M2/68‑1 (MQ=255)
                                               ggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGc                                                             >  1:478235‑M2/67‑139 (MQ=255)
                                               ggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGc                                                             >  2:331747‑M2/67‑139 (MQ=255)
                                               ggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGc                                                             >  2:399366‑M2/67‑139 (MQ=255)
                                                 cttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCa                                                           <  1:472626‑M2/75‑1 (MQ=255)
                                                              aattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgtg                                                                                                           <  1:555869‑M2/27‑1 (MQ=255)
                                                              aattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgtg                                                                                                           >  2:555869‑M2/52‑78 (MQ=255)
                                                                                 attccactacctttgttagaggaaagtgctagATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttgt                           <  1:296776‑M2/107‑1 (MQ=255)
                                                                                               ggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGaaa                                                                              <  2:338499‑M2/56‑1 (MQ=255)
                                                                                               ggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGaaa                                                                              >  1:338499‑M2/19‑74 (MQ=255)
                                                                                               ggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGc             <  2:280055‑M2/121‑1 (MQ=255)
                                                                                                gtagaggaaagtgctaaATAATAATCAATTGTTAAACTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCa            <  1:399366‑M2/122‑1 (MQ=255)
                                                                                                    aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTa                                                    >  2:346459‑M2/14‑95 (MQ=255)
                                                                                                    aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTa                                                    <  1:346459‑M2/82‑1 (MQ=255)
                                                                                                    aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCt        >  2:635757‑M2/14‑139 (MQ=255)
                                                                                                      gaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCtat      <  1:635757‑M2/128‑1 (MQ=255)
                                                                                                          gtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGAt  >  1:246788‑M2/8‑139 (MQ=255)
                                                                                                                 |                                                                                                                                   
CTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGAT  >  NC_000913/1300585‑1300829

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