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–1299500 2–1199 20 [0] [0] 21 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)18 (0.730) 18/248 0.5 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

GATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAA  >  NC_000913/1300556‑1300846
                                                                                                                                              |                                                                                                                                                    
aattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaataataa                                                                                                                                                <  4:53762‑M2/7‑1 (MQ=255)
                      ttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCa                                                                                                                          >  5:27654‑M2/121‑149 (MQ=255)
                                                atgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgtg                                                                                                                            <  8:59375‑M2/27‑1 (MQ=255)
                                                atgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATtgtg                                                                                                                            >  7:59375‑M2/95‑121 (MQ=255)
                                                      ttttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATa                                                                                          <  1:72641‑M2/61‑1 (MQ=255)
                                                          cgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGAc                                                                                       <  6:27654‑M2/64‑1 (MQ=255)
                                                             cacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAAc                                                                                                           <  3:77523‑M2/44‑1 (MQ=255)
                                                             cacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAAc                                                                                                           >  4:77523‑M2/82‑125 (MQ=255)
                                                                          tgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAAtctc                                                                       >  1:28240‑M2/69‑148 (MQ=255)
                                                                                       gaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGattgaatt                                                         <  2:28240‑M2/94‑1 (MQ=255)
                                                                                                                                 aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGc                                                                              >  5:118022‑M2/14‑86 (MQ=255)
                                                                                                                                 aggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGc                                                                              <  6:118022‑M2/73‑1 (MQ=255)
                                                                                                                                     aagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATAc           >  7:56663‑M2/10‑149 (MQ=255)
                                                                                                                                       gtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGattgaattga                                                       >  5:16766‑M2/8‑103 (MQ=255)
                                                                                                                                       gtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGattgaattga                                                       <  6:16766‑M2/96‑1 (MQ=255)
                                                                                                                                        tgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGt        >  3:68042‑M2/7‑149 (MQ=255)
                                                                                                                                         gctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATgg                                  >  5:97665‑M2/6‑121 (MQ=255)
                                                                                                                                         gctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATgg                                  <  6:97665‑M2/116‑1 (MQ=255)
                                                                                                                                             aATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttgtt                                           >  2:103902‑M2/2‑109 (MQ=255)
                                                                                                                                             aATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGAtgtttgtt                                           <  1:103902‑M2/108‑1 (MQ=255)
                                                                                                                                              ataataatCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATaa  <  4:66481/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAATAA  >  NC_000913/1300556‑1300846

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: