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

ACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGG  >  NC_000913/2846251‑2846349
                                                 |                                                 
aCGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGA                                                   <  1:1362640‑M1/50‑1 (MQ=255)
  ggCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAca                                                 <  1:712924/50‑1 (MQ=255)
   gCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtca                                                >  1:3569986‑M1/1‑47 (MQ=255)
   gCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtca                                                >  1:3685947‑M1/1‑47 (MQ=255)
   gCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAcag                                                >  1:2819614/1‑50 (MQ=255)
     gTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatc                                              >  1:4296574‑M1/1‑45 (MQ=255)
      tAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatcc                                             <  1:2249968‑M1/50‑7 (MQ=255)
      tAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAcagcca                                             <  1:1716254/50‑1 (MQ=255)
        aTGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccgg                                           <  1:2222873‑M1/50‑9 (MQ=255)
        aTGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAcagccagc                                           <  1:3530052/50‑1 (MQ=255)
         tGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAcagccagc                                           >  1:236377/1‑49 (MQ=255)
          gACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcg                                         >  1:2184922‑M1/1‑40 (MQ=255)
           aCCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcgg                                        <  1:1113096‑M1/50‑12 (MQ=255)
             cATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggc                                       >  1:3387138‑M1/1‑37 (MQ=255)
              aTTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCt                                     <  1:1827100/50‑1 (MQ=255)
              aTTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCt                                     <  1:1368227/50‑1 (MQ=255)
               ttCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgag                                    <  1:2520831‑M1/50‑16 (MQ=255)
               ttCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggc                                       >  1:3499965‑M1/1‑35 (MQ=255)
               ttCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggc                                       >  1:2306482‑M1/1‑35 (MQ=255)
               ttCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTg                                    >  1:54989/1‑50 (MQ=255)
                 cAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagag                                  >  1:3171632‑M1/1‑33 (MQ=255)
                 cAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgaga                                   <  1:126653‑M1/49‑17 (MQ=255)
                  aGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGttt                                 <  1:3467499/50‑1 (MQ=255)
                   gTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTg                                >  1:1854837/1‑50 (MQ=255)
                   gTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTg                                <  1:3156179/50‑1 (MQ=255)
                    tttGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggc                                >  1:1234501‑M1/1‑30 (MQ=255)
                     ttGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcga                              >  1:2660193‑M1/1‑29 (MQ=255)
                      tGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTgta                              <  1:2312371/49‑1 (MQ=255)
                       gCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgaca                            >  1:3038532‑M1/1‑27 (MQ=255)
                       gCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTgtag                             >  1:3683884/1‑49 (MQ=255)
                        cGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTgtagta                           <  1:1725002/50‑1 (MQ=255)
                         gTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgacatc                          >  1:4286559‑M1/1‑25 (MQ=255)
                         gTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTgtagta                           >  1:314813/1‑49 (MQ=255)
                          ttCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTgtagta                           >  1:2047538/1‑48 (MQ=255)
                          ttCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTaaa                         <  1:3465779/50‑1 (MQ=255)
                            cgtcgtTACCAAACAGCACCGAtcatccggcggcgagagtcgacatctc                        >  1:3594496‑M1/1‑22 (MQ=255)
                            cgtcgtTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAg                        <  1:798602/49‑1 (MQ=255)
                                gtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacg                    >  1:50154‑M1/1‑18 (MQ=255)
                                gtTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACt                   >  1:4038683/1‑50 (MQ=255)
                                 ttACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacg                    >  1:2650602‑M1/1‑17 (MQ=255)
                                 ttACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTc                  <  1:3608449/50‑1 (MQ=255)
                                  tACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgcag                 >  1:2688416‑M1/1‑16 (MQ=255)
                                  tACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgcag                 >  1:1560148‑M1/1‑16 (MQ=255)
                                  tACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTc                  >  1:3461726/1‑49 (MQ=255)
                                   aCCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTcca                >  1:3324879/1‑50 (MQ=255)
                                    ccAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgcagt                <  1:2851489‑M1/49‑36 (MQ=255)
                                    ccAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTccac               >  1:967032/1‑50 (MQ=255)
                                    ccAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTcca                <  1:216964/49‑1 (MQ=255)
                                     cAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTccacc              <  1:243503/50‑1 (MQ=255)
                                      aaaCAGCACCGAtcatccggcggcgagaggcgacatctccacgcagtgcg             >  1:2818394‑M1/1‑12 (MQ=255)
                                      aaaCAGCACCGAtcatccggcggcgagaggcgacatctccacgcagtgcg             <  1:2125655‑M1/50‑39 (MQ=255)
                                      aaaCAGCACCGAtcatccggcggcgagaggcgacatctccacgcagtgcg             >  1:1744800‑M1/1‑12 (MQ=255)
                                       aaCAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCcaccac            <  1:105525/50‑1 (MQ=255)
                                        aCAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACt           <  1:2587745/50‑1 (MQ=255)
                                          aGCACCGAtcatccggcggcgagaggcgacatctccacgcagtgcgcgac         >  1:573178‑M1/1‑8 (MQ=255)
                                           gCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTcc        <  1:3363501/50‑1 (MQ=255)
                                            cACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCg       <  1:1942458/50‑1 (MQ=255)
                                              ccGAtcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgc     <  1:1279705‑M1/50‑47 (MQ=255)
                                              ccGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGc     <  1:1994040/50‑1 (MQ=255)
                                               cGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCa    <  1:2107245/50‑1 (MQ=255)
                                                gAtcatccggcggcgagaggcgacatctccacgcagcgcgcgacgcgca    >  1:1549978‑M1/1‑2 (MQ=255)
                                                gACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAg   >  1:4235774/1‑50 (MQ=255)
                                                 atcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaa   >  1:1717430‑M1/1‑1 (MQ=255)
                                                 aCAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAgg  >  1:338781/1‑50 (MQ=255)
                                                 |                                                 
ACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGG  >  NC_000913/2846251‑2846349

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