New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_012967 = 54333165 (0.230)7 (0.030) 7/126 17.1 15.5% coding (549/552 nt) ybcL Raf kinase inhibitor‑like protein YbcL
?NC_012967 = 543565 19 (0.070)coding (222/798 nt) ybcM helix‑turn‑helix transcriptional regulator
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCA  >  NC_012967/543261‑543402
                                                                      |                                                                       
cGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGt                                                                        <  1:1469927/72‑1 (MQ=255)
cGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGt                                                                        <  2:7778787/72‑1 (MQ=255)
 gTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTa                                                                       >  2:9501891/1‑72 (MQ=255)
   tAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAgg                                                                     <  2:5174571/72‑1 (MQ=255)
   tAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAgg                                                                     <  2:8871388/72‑1 (MQ=255)
    aGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAggg                                                                    <  2:8477095/72‑1 (MQ=255)
     gTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGt                                                                   <  1:8509907/72‑1 (MQ=255)
     gTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGt                                                                   <  2:5534175/72‑1 (MQ=255)
     gTTGGTTATATGCTTAANGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGt                                                                   >  1:883643/1‑72 (MQ=255)
      ttGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTg                                                                  <  1:5360602/72‑1 (MQ=255)
      tgggTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTg                                                                  <  2:418910/70‑1 (MQ=255)
       tGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTga                                                                 <  2:6244307/72‑1 (MQ=255)
          ttATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGagag                                                              <  1:1289251/72‑1 (MQ=255)
           tataTGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGt                                                             >  1:3849452/1‑72 (MQ=255)
           tataTGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGt                                                             >  1:3880591/1‑72 (MQ=255)
            atatGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTa                                                            >  2:7615417/1‑72 (MQ=255)
            atatGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTa                                                            >  1:4189158/1‑72 (MQ=255)
            atatGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTa                                                            <  1:4862476/72‑1 (MQ=255)
             tatGCTTAATGCTAATGAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTAt                                                           <  2:371993/72‑1 (MQ=255)
              atGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATg                                                          <  2:6239806/72‑1 (MQ=255)
                 cTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTg                                                       <  2:5879751/72‑1 (MQ=255)
                   tAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGc                                                     <  1:5179535/72‑1 (MQ=255)
                    aaTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCa                                                    <  1:8105606/72‑1 (MQ=255)
                     aTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCaa                                                   <  1:9719155/72‑1 (MQ=255)
                     aTGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCaa                                                   <  1:1654399/72‑1 (MQ=255)
                      tGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAag                                                  <  2:9520342/72‑1 (MQ=255)
                      tGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAag                                                  <  2:4134519/72‑1 (MQ=255)
                       gCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAaga                                                 >  2:1932790/1‑72 (MQ=255)
                       acTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAaga                                                 <  2:6011072/71‑1 (MQ=255)
                        cTAATAAAATCGCAACCGCTGAGATANCACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAagag                                                <  2:1056881/72‑1 (MQ=255)
                        cTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAagag                                                <  1:2624830/72‑1 (MQ=255)
                        cTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAagag                                                <  2:4986559/72‑1 (MQ=255)
                           ataaAAGCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTa                                             <  1:2086537/72‑1 (MQ=255)
                            taaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTaa                                            <  1:6981827/72‑1 (MQ=255)
                            taaAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTaa                                            <  2:9825662/72‑1 (MQ=255)
                            aaaaaNCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTaa                                            <  2:876226/71‑1 (MQ=255)
                              aaaTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGa                                          <  1:1449071/72‑1 (MQ=255)
                              aaaTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGa                                          <  1:1947862/72‑1 (MQ=255)
                               aaTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGAc                                         >  1:7677149/1‑72 (MQ=255)
                               aaTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGAc                                         >  2:5868151/1‑72 (MQ=255)
                               aaTCCCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGAc                                         <  2:5633213/72‑1 (MQ=255)
                                ctCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                        <  2:3481978/71‑1 (MQ=255)
                                aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                        <  2:9064495/72‑1 (MQ=255)
                                aTCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACt                                        <  1:8187501/72‑1 (MQ=255)
                                 tCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTa                                       <  2:8996319/72‑1 (MQ=255)
                                  cGCAACCGCTGAGATAACNCCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTaa                                      <  1:3819314/72‑1 (MQ=255)
                                    cAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACt                                    <  1:3575758/72‑1 (MQ=255)
                                     aaCCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACtt                                   <  1:3337066/72‑1 (MQ=255)
                                     aaCCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACtt                                   <  1:6190851/72‑1 (MQ=255)
                                     aaCCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACtt                                   <  2:154178/72‑1 (MQ=255)
                                           tGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCANGAGGTAAGACTAACTTAAAGAt                             >  1:4357629/1‑72 (MQ=255)
                                           tGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGAt                             >  1:6820734/1‑72 (MQ=255)
                                           tGAGATAACACCAGTTTATAGGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGAt                             <  1:4024857/72‑1 (MQ=255)
                                               aTAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCgaa                         <  2:1386513/72‑1 (MQ=255)
                                                  ncacCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCgaagaa                      >  1:5394969/2‑72 (MQ=255)
                                                   cacCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCCAGAGGTAAGACTAACTTAAAGATCGAAGaaa                     >  1:328938/1‑72 (MQ=255)
                                                    acCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAAt                    >  1:8953937/1‑72 (MQ=255)
                                                     ccAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATa                   >  1:2818626/1‑72 (MQ=255)
                                                     ccAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATa                   >  1:1744088/1‑72 (MQ=255)
                                                            aTGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGc            >  1:4603867/1‑72 (MQ=255)
                                                            aTGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGc            >  1:3363370/1‑72 (MQ=255)
                                                            aTGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGc            <  2:4677182/72‑1 (MQ=255)
                                                             tGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCa           <  2:3521969/72‑1 (MQ=255)
                                                                    aaGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAcat    <  2:9797674/72‑1 (MQ=255)
                                                                      gTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAAcatca  <  2:8189502/72‑1 (MQ=255)
                                                                      |                                                                       
CGTTAGTTGGTTATATGCTTAATGCTAATAAAATCGCAACCGCTGAGATAACACCAGTTTATGAGATAAAGTAGGGTGAGAGTATGCTGGCAAGAGGTAAGACTAACTTAAAGATCGAAGAAATACGGATGCATAAACATCA  >  NC_012967/543261‑543402

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