New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 284630031 (0.430)23 (0.340) 21/92 1.2 40.7% coding (172/1710 nt) hycE formate hydrogenlyase subunit HycE
?NC_000913 2867269 = 38 (0.560)coding (283/993 nt) rpoS RNA polymerase, sigma S (sigma 38) factor

ACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGG  >  NC_000913/2846251‑2846349
                                                 |                                                 
aCGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGA                                                   >  1:3548705‑M1/1‑50 (MQ=255)
aCGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGA                                                   >  1:4390599‑M1/1‑50 (MQ=255)
  ggCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAt                                                  >  1:2164690‑M1/1‑48 (MQ=255)
  ggCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAca                                                 <  1:3488393/50‑1 (MQ=255)
  ggCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGA                                                   >  1:1240201‑M1/1‑48 (MQ=255)
    cGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcat                                               >  1:4209574‑M1/1‑46 (MQ=255)
    cGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcat                                               <  1:3844335‑M1/50‑5 (MQ=255)
      tAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatcc                                             >  1:3216006‑M1/1‑44 (MQ=255)
      tAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatcc                                             >  1:4148015‑M1/1‑44 (MQ=255)
      tAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAcagcca                                             <  1:2725373/50‑1 (MQ=255)
        aTGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccgg                                           <  1:1128710‑M1/50‑9 (MQ=255)
        aTGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAcagccagc                                           >  1:2170929/1‑50 (MQ=255)
        aTGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAcagccagc                                           >  1:2920784/1‑50 (MQ=255)
         tGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggc                                          <  1:2224265‑M1/50‑10 (MQ=255)
          gACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcg                                         >  1:167760‑M1/1‑40 (MQ=255)
            ccATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggc                                       >  1:4082540‑M1/1‑38 (MQ=255)
             cATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCAccc                                      >  1:2503811/1‑50 (MQ=255)
              aTTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCt                                     >  1:3399/1‑50 (MQ=255)
              aTTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCt                                     >  1:5032512/1‑50 (MQ=255)
                tCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGt                                   >  1:1035767/1‑50 (MQ=255)
                tCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGt                                   >  1:2690847/1‑50 (MQ=255)
                 cAGTTTGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagag                                  >  1:3873008‑M1/1‑33 (MQ=255)
                 cAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGt                                   <  1:2552070/49‑1 (MQ=255)
                  aGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGttt                                 >  1:3596120/1‑50 (MQ=255)
                  aGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGttt                                 >  1:3187581/1‑50 (MQ=255)
                   gTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTg                                <  1:3025770/50‑1 (MQ=255)
                      tGCGTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcga                              >  1:1742197‑M1/1‑28 (MQ=255)
                      tGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTgtag                             >  1:4590766/1‑50 (MQ=255)
                      tGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTgtag                             >  1:3227300/1‑50 (MQ=255)
                      tGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTgta                              <  1:741449/49‑1 (MQ=255)
                        cGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTgtagt                            >  1:1242653/1‑49 (MQ=255)
                        cGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTgtagt                            <  1:1660968/49‑1 (MQ=255)
                         gTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgacatc                          >  1:1261417‑M1/1‑25 (MQ=255)
                         gTTCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgacatc                          >  1:1740016‑M1/1‑25 (MQ=255)
                          ttCGTCGTTACCAAACAGCACCGAtcatccggcggcgagaggcgacatct                         >  1:211586‑M1/1‑24 (MQ=255)
                            cgtcgtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatctc                        >  1:3204020‑M1/1‑22 (MQ=255)
                             gtcgtTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGaa                      <  1:4613901/50‑1 (MQ=255)
                             gtcgtTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGa                       <  1:1883517/49‑1 (MQ=255)
                             gtcgtTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGa                       >  1:4054291/1‑49 (MQ=255)
                              tcgtTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGaaa                     >  1:3573860/1‑50 (MQ=255)
                               cgtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacg                    >  1:4064180‑M1/1‑19 (MQ=255)
                                gtTACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacg                    >  1:3743780‑M1/1‑18 (MQ=255)
                                gtTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAAc                    <  1:1170475/49‑1 (MQ=255)
                                 ttACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTc                  >  1:2402089/1‑50 (MQ=255)
                                  tACCAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgcag                 >  1:3204152‑M1/1‑16 (MQ=255)
                                    ccAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTccac               >  1:2075334/1‑50 (MQ=255)
                                     cAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgcagtgc              >  1:4982883‑M1/1‑13 (MQ=255)
                                     cAAACAGCACCGAtcatccggcggcgagaggcgacatctccacgcagtgc              <  1:675099‑M1/50‑38 (MQ=255)
                                     cAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTccac               >  1:5564329/1‑49 (MQ=255)
                                       aaCAGCACCGAtcatccggcggcgagaggcgacatctccacgcagtgcgc            <  1:1008855‑M1/50‑40 (MQ=255)
                                        aCAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACt           <  1:2540477/50‑1 (MQ=255)
                                          aGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTc         <  1:5074961/50‑1 (MQ=255)
                                           gCACCGAtcatccggcggcgagaggcgacatctccacgcagtgcgcgacg        >  1:1587528‑M1/1‑7 (MQ=255)
                                             aCCGAtcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcg      <  1:49869‑M1/50‑46 (MQ=255)
                                              ccGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCgg      >  1:123637/1‑49 (MQ=255)
                                               cGAtcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgc     >  1:2415825‑M1/1‑3 (MQ=255)
                                               cGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCgg      >  1:3778754/1‑48 (MQ=255)
                                                gAtcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaa   <  1:497930‑M1/50‑49 (MQ=255)
                                                 atcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaaa  <  1:43522‑M1/50‑50 (MQ=255)
                                                 atcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaa   >  1:1096114‑M1/1‑1 (MQ=255)
                                                 atcatccggcggcgagaggcgacatctccacgcagtgcgcgacgcgcaa   >  1:2355170‑M1/1‑1 (MQ=255)
                                                 aCAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAgg  <  1:1611072/50‑1 (MQ=255)
                                                 |                                                 
ACGGCGTAATGACCATTCAGTTTGCGTTCGTCGTTACCAAACAGCACCGACAGCCAGCCACCCTGTTTGTAGTAAAGAAACTCCACCACTTCCGGCAGG  >  NC_000913/2846251‑2846349

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