New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 28463000 (0.000)24 (0.940) 21/90 0.1 100% coding (172/1710 nt) hycE formate hydrogenlyase subunit HycE
?NC_000913 2867269 = 0 (0.000)coding (283/993 nt) rpoS RNA polymerase, sigma S (sigma 38) factor

GCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAA  >  NC_000913/2867222‑2867316
                                               |                                               
gcgtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatc                                                <  1:1338542‑M2/2‑1 (MQ=255)
 cgtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatca                                               >  1:1771169‑M2/47‑49 (MQ=255)
  gtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatcat                                              >  1:1441831‑M2/46‑49 (MQ=255)
          cattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                      >  1:294269‑M2/38‑49 (MQ=255)
            ttcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                      >  1:58076‑M2/36‑47 (MQ=255)
              cagtttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCga                                    >  1:1787802‑M2/34‑47 (MQ=255)
              cagtttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCga                                    >  1:2021656‑M2/34‑47 (MQ=255)
                 tttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGc                               <  1:253139‑M2/19‑1 (MQ=255)
                    gcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGa                             >  1:195113‑M2/28‑48 (MQ=255)
                      gttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAt                          >  1:870919‑M2/26‑49 (MQ=255)
                         cgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtctc                       >  1:1548915‑M2/23‑49 (MQ=255)
                         cgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtct                        <  1:250209‑M2/26‑1 (MQ=255)
                          gtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtctc                       >  1:2020788‑M2/22‑48 (MQ=255)
                          gtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtctc                       >  1:392546‑M2/22‑48 (MQ=255)
                            cgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCAc                    >  1:152286‑M2/20‑49 (MQ=255)
                             gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCAc                    <  1:785302‑M2/30‑1 (MQ=255)
                             gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                   >  1:1457886‑M2/19‑49 (MQ=255)
                             gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                   >  1:1378421‑M2/19‑49 (MQ=255)
                               taccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                   >  1:34960‑M2/17‑47 (MQ=255)
                                 ccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGt               <  1:773073‑M2/35‑1 (MQ=255)
                                     acagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgcgc           >  1:1636963‑M2/11‑49 (MQ=255)
                                       agcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGa         >  1:56729‑M2/9‑49 (MQ=255)
                                        gcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAc        >  1:1870676‑M2/8‑49 (MQ=255)
                                        gcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAc        <  1:438479‑M2/42‑1 (MQ=255)
                                              aTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaa  <  1:200518‑M2/48‑1 (MQ=255)
                                               tcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaa  >  1:1554893/1‑48 (MQ=255)
                                               |                                               
GCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAA  >  NC_000913/2867222‑2867316

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