Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913_3_pae_tpiA 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_3_pae_tpiA 1299499–1300697 1300697 1–1199 23 [0] [0] 23 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_3_pae_tpiA = 12994980 (0.000)22 (0.670) 20/274 0.7 100% intergenic (+253/‑1684) ychE/oppA UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein
?NC_000913_3_pae_tpiA 1300698 = 0 (0.000)intergenic (+1453/‑484) ychE/oppA UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

GGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAA  >  NC_000913_3_pae_tpiA/1300553‑1300843
                                                                                                                                                 |                                                                                                                                                 
tcaaattacagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaataat                                                                                                                                               >  1:404994‑M2/146‑150 (MQ=255)
        cagttagttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAAtt                                                                                                                                       >  1:569587‑M2/138‑150 (MQ=255)
             agttataaggatttccttaattgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTaa                                                                                                                                  <  1:58820‑M2/18‑1 (MQ=255)
               ttataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAat                                                                                                                                >  2:38241‑M2/131‑150 (MQ=255)
                tataaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAatt                                                                                                                               <  1:452077‑M2/21‑1 (MQ=255)
                   aaggatttccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAattatt                                                                                                                            <  1:38241‑M2/24‑1 (MQ=255)
                           ccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAttt                                                                                                                    <  2:505364‑M2/32‑1 (MQ=255)
                           ccttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCAttt                                                                                                                    <  1:580754‑M2/32‑1 (MQ=255)
                            cttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagcggaaaatgctaaATAATAATCAATTGTTAaattattctgactttt                                                                                                                   >  2:242989‑M2/118‑144 (MQ=255)
                            cttaactgcttctcctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATGGTTAaattattttgactttt                                                                                                                   >  2:422098‑M2/118‑144 (MQ=255)
                                        tcctcaccatcatggtattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACt                                                                                                       <  1:51413‑M2/45‑1 (MQ=255)
                                          ctcaccatcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCACTTGTTAAATTATTTCAAAGTTCACTACTGGAACTGt                                                                                                     >  2:593514‑M2/104‑150 (MQ=255)
                                               catcatgttattttcgccacatcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCa                                                                                                >  1:416567‑M2/99‑150 (MQ=255)
                                                        ttgttcggcacatcataagccttggattggtgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATa                                                                                       <  2:579161‑M2/61‑1 (MQ=255)
                                                                   atcataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTCCATTTCACTACTGGAACTGTACTCAGAAAAGATAGACATTCTtac                                                                            >  2:407764‑M2/79‑149 (MQ=255)
                                                                      ataatcctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGAATTTAACTACTGGAACTGTAATCAGAAAAGATAGAGATGAttaggta                                                                         >  1:65054‑M2/76‑147 (MQ=255)
                                                                        aaacctgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAAt                                                                       <  2:337297‑M2/77‑1 (MQ=255)
                                                                          tcctggacttgctgaagaataaatgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTATTAAATTATTGTGCATTTCACAACTGGAACTGTAATCAGAATAGATAGACATGCTTAGCCAAtct                                                                     <  2:567297‑M2/79‑1 (MQ=255)
                                                                             tgggcttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTa                                                                  <  2:247173‑M2/82‑1 (MQ=255)
                                                                                  ttgctggtgaagaatcgaaatgatattattaattgcaccgtctctggtagagtaaagtgctaaATAATAATCAATTGTTACATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGTCATGCTTACCCCATCTCTATtgga                                                             <  2:438506‑M2/87‑1 (MQ=255)
                                                                                  ttgctgaagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGCTCGACCTGCGTAGCAAATCTCtatttca                                                             >  2:566933‑M2/64‑148 (MQ=255)
                                                                                        aagaataattgaaatgatattattaattccactgcctttggtagaggaaagtgctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGattgaat                                                       <  1:24839‑M2/93‑1 (MQ=255)
                                                                                                                                             ctaaATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCaa  <  2:437841‑M2/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATAATAATCAATTGTTAAATTATTGTGCATTTCACTACTGGAACTGTAATCAGAAAAGATAGACATGCTTAGCCAATCTCTATTTGATTGAATTGAAAGATGTTTGTTAAGGCATGGATGCAAGCTATAGATTCTGATACGGTCAA  >  NC_000913_3_pae_tpiA/1300553‑1300843

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