New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 28463000 (0.000)23 (0.900) 18/90 0.2 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

CGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCA  >  NC_000913/2846255‑2846347
                                             |                                               
cGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtca                                              >  1:1632373‑M1/1‑46 (MQ=255)
cGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtca                                              >  1:668056‑M1/1‑46 (MQ=255)
 gTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcat                                             >  1:142263‑M1/1‑45 (MQ=255)
    aTGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccg                                          <  1:164465‑M1/49‑8 (MQ=255)
        ccATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcgg                                      <  1:2040853‑M1/49‑12 (MQ=255)
          aTTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggc                                     >  1:1336118‑M1/1‑36 (MQ=255)
              aGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagag                                <  1:393157‑M1/49‑18 (MQ=255)
              aGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgaga                                 <  1:1278745‑M1/48‑17 (MQ=255)
                tttGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggc                              <  1:1586795‑M1/49‑20 (MQ=255)
                    cGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgaca                          >  1:373341‑M1/1‑26 (MQ=255)
                    cGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgaca                          <  1:1782740‑M1/49‑24 (MQ=255)
                     gTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgacat                         <  1:692587‑M1/49‑25 (MQ=255)
                     gTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgacat                         <  1:1746917‑M1/49‑25 (MQ=255)
                      ttCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgacatc                        >  1:1941028‑M1/1‑24 (MQ=255)
                      ttCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgacatc                        >  1:244377‑M1/1‑24 (MQ=255)
                        cgtcgtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatct                       <  1:1555228‑M1/48‑27 (MQ=255)
                           cgtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccac                   >  1:1742014‑M1/1‑19 (MQ=255)
                           cgtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccac                   >  1:318230‑M1/1‑19 (MQ=255)
                            gtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacg                  <  1:793261‑M1/49‑32 (MQ=255)
                                 cAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgcagtg             >  1:2043546‑M1/1‑13 (MQ=255)
                                  aaaCAGCACCGAtcatccggcggcgagaggcgacatctccacgcagtgc            >  1:1939755‑M1/1‑12 (MQ=255)
                                     cAGCACCGAtcatccggcggcgagaggcgacatctccacgcagtgcgc          >  1:588149‑M1/1‑9 (MQ=255)
                                     cAGCACCGAtcatccggcggcgagaggcgacatctccacgcagtgcgc          >  1:1227811‑M1/1‑9 (MQ=255)
                                           cGAtcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcg    >  1:756937‑M1/1‑3 (MQ=255)
                                             atcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgca  >  1:1826292‑M1/1‑1 (MQ=255)
                                             |                                               
CGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCA  >  NC_000913/2846255‑2846347

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.