New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 284630033 (0.570)28 (0.520) 24/92 0.6 48.7% coding (172/1710 nt) hycE formate hydrogenlyase subunit HycE
?NC_000913 2867269 = 28 (0.520)coding (283/993 nt) rpoS RNA polymerase, sigma S (sigma 38) factor

GCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAA  >  NC_000913/2867222‑2867318
                                               |                                                 
gcgtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatca                                                 >  1:3569986‑M2/48‑50 (MQ=255)
gcgtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatca                                                 >  1:3685947‑M2/48‑50 (MQ=255)
gCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGatca                                                 >  1:582729/1‑50 (MQ=255)
 ccATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGatcat                                                <  1:2897259/50‑1 (MQ=255)
 ccATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGatca                                                 <  1:658455/49‑1 (MQ=255)
  gtaatgaccattcagtttgcgttcgtcgttaccaaacagcaccgatcatc                                               >  1:4296574‑M2/46‑50 (MQ=255)
  cATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGatcatc                                               <  1:357433/50‑1 (MQ=255)
  cATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGatcatc                                               >  1:412378/1‑50 (MQ=255)
   taatgaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATcc                                              <  1:2249968‑M2/6‑1 (MQ=255)
   aTAACGGCGGGCAATTTTTACCACCAGACGCAAATTACTCTCGATCATcc                                              <  1:2424315/50‑1 (MQ=255)
     atgaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCcgg                                            <  1:2222873‑M2/8‑1 (MQ=255)
      aCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCcggc                                           <  1:1569049/50‑1 (MQ=255)
       gaccattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCcggcg                                          >  1:2184922‑M2/41‑50 (MQ=255)
        accattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCcggcgg                                         <  1:1113096‑M2/11‑1 (MQ=255)
          cggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCggcggc                                        >  1:591117/1‑49 (MQ=255)
          cattcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                        >  1:3387138‑M2/38‑49 (MQ=255)
           gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCggcggc                                        >  1:3152236/1‑48 (MQ=255)
           gggCAATTTTTAACACCAGACGCAAGTTAATCTCGATCATCCggcggc                                        >  1:311874/1‑48 (MQ=39)
            ttcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                        >  1:2306482‑M2/36‑47 (MQ=255)
            ttcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCggcggc                                        >  1:3499965‑M2/36‑47 (MQ=255)
            ttcagtttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCgag                                     <  1:2520831‑M2/15‑1 (MQ=255)
             gCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCgaga                                    <  1:3135841/50‑1 (MQ=255)
              cagtttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCgaga                                    <  1:126653‑M2/16‑1 (MQ=255)
              cagtttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGagag                                   >  1:3171632‑M2/34‑50 (MQ=255)
               aaTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAgg                                  >  1:3007303/1‑50 (MQ=255)
               aaTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAgg                                  >  1:1358707/1‑50 (MQ=255)
               aaTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAgg                                  >  1:1884875/1‑50 (MQ=255)
                 tttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGc                                 >  1:1234501‑M2/31‑49 (MQ=255)
                  ttgcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGa                               >  1:2660193‑M2/30‑50 (MQ=255)
                    ttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGAc                              <  1:3066697/49‑1 (MQ=255)
                    gcgttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACa                             >  1:3038532‑M2/28‑50 (MQ=255)
                      gttcgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACAtc                           >  1:4286559‑M2/26‑50 (MQ=255)
                        caccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACAtctc                         >  1:2317416/1‑50 (MQ=255)
                         cgtcgttaccaaacagcaccgaTCATCCGGCGGCGAGAGTCGACAtctc                         >  1:3594496‑M2/23‑49 (MQ=255)
                            aGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACg                     <  1:2316634/50‑1 (MQ=255)
                            aGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACg                     <  1:3001065/50‑1 (MQ=255)
                             gttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                     >  1:50154‑M2/19‑49 (MQ=255)
                              ttaccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACg                     >  1:2650602‑M2/18‑48 (MQ=255)
                              aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACg                     >  1:449879/1‑48 (MQ=255)
                               taccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAg                  >  1:2688416‑M2/17‑50 (MQ=255)
                               taccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAg                  >  1:1560148‑M2/17‑50 (MQ=255)
                                gCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGt                 >  1:1603836/1‑50 (MQ=255)
                                gCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGt                 >  1:1503699/1‑50 (MQ=255)
                                 ccaaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGt                 <  1:2851489‑M2/35‑1 (MQ=255)
                                  aaGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgc               >  1:512264/1‑50 (MQ=255)
                                   aaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgcg              >  1:2818394‑M2/13‑50 (MQ=255)
                                   aaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgcg              <  1:2125655‑M2/38‑1 (MQ=255)
                                   aaacagcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgcg              >  1:1744800‑M2/13‑50 (MQ=255)
                                   aGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTgcg              >  1:1617699/1‑50 (MQ=255)
                                       agcaccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAc          >  1:573178‑M2/9‑50 (MQ=255)
                                       aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAc          >  1:3116953/1‑50 (MQ=255)
                                       aCTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAc          >  1:128321/1‑50 (MQ=255)
                                        ctctCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAcg         <  1:4002455/50‑1 (MQ=255)
                                        ctctCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAc          <  1:2043063/49‑1 (MQ=255)
                                         tctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGAcgc        >  1:1638457/1‑50 (MQ=255)
                                           ccgaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACgcgc      <  1:1279705‑M2/46‑1 (MQ=255)
                                             gaTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGCGCGCGACGCGCa     >  1:1549978‑M2/3‑49 (MQ=255)
                                              aTCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaa    >  1:1717430‑M2/2‑49 (MQ=255)
                                               tcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaaaa  <  1:3214280/50‑1 (MQ=255)
                                               tcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCaaaa  >  1:1227203/1‑50 (MQ=255)
                                               |                                                 
GCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAA  >  NC_000913/2867222‑2867318

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