New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 731453NA (NA)7 (0.050) 7/142 12.9 NA coding (2648/4194 nt) rhsC RhsC
?NC_000913 = 731461 NA (NA)coding (2656/4194 nt) rhsC RhsC
Rejected: Coverage evenness skew score above cutoff.

TACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAA  >  NC_000913/731434‑731528
                   |                                                                           
gccagctccagcACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGtata                     <  2:2923555/74‑1 (MQ=9)
  cagctacagcACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGCATAcc                   >  2:4645452/1‑76 (MQ=21)
  cagcgacagcACCCCCGGCTGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATAcc                   <  1:964682/76‑1 (MQ=11)
   agctacagcACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGCATCCGGATATCCGCATCCCGcataacc                  >  1:1641732/1‑73 (MQ=11)
    gctacagcACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACcac                 >  1:3039070/1‑76 (MQ=34)
     ctacagcACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACcaca                >  2:2873903/1‑76 (MQ=34)
     ctacagcACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACcaca                <  1:5603222/76‑1 (MQ=35)
     cgacagcACCCCCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACcaca                <  2:3764317/74‑1 (MQ=255)
       acagcACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGa              >  2:1856429/1‑76 (MQ=34)
        cagcACCCCCGGCGGGCGGACCGGCGTTCACACCCCCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGAc             <  2:5602381/76‑1 (MQ=11)
        cagcACCACCGGCAGGCTGACCGGCGTTCACCCCACCGCAGCGAATCTGGATATCCGCATCCCGCATACCACAGAc             >  1:1539694/1‑76 (MQ=255)
         cgcACCACCTGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATAACACAGAcc            <  1:717780/75‑1 (MQ=11)
         agcCCCCCCGGCGGGCTGACCGGCGTTCACACCCCCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGAcc            <  2:3506992/76‑1 (MQ=11)
         agcACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGAcc            >  1:4955542/1‑76 (MQ=35)
         agcACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGAcc            <  1:3779856/76‑1 (MQ=35)
         agcACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGAcc            <  1:920757/76‑1 (MQ=35)
          gcACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGAccc           <  2:3466713/76‑1 (MQ=35)
           ccccaccGGCGGGCGGCCCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg          <  2:5119990/74‑1 (MQ=11)
           caccgcCGGCAGGCTGAGCGGCGTGCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg          <  2:876768/76‑1 (MQ=11)
           caccacCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg          >  2:2531236/1‑76 (MQ=34)
           caccacCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCg          <  1:995995/76‑1 (MQ=35)
            cccaccGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCgg         <  1:2242605/75‑1 (MQ=25)
            accaccGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCGGGATATCCGCATCCCGTATACCACAGACCCgg         >  2:1261806/1‑76 (MQ=21)
             ccaccGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGc        >  2:3209010/1‑76 (MQ=35)
             ccaccGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGc        >  2:4352600/1‑76 (MQ=34)
              caccGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCa       <  1:5594755/76‑1 (MQ=35)
               cccGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAg      <  1:5728742/75‑1 (MQ=25)
               accGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAg      >  1:1324209/1‑76 (MQ=34)
                ccGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAgg     >  1:5620905/1‑76 (MQ=35)
                  ggCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTa   >  2:2154791/1‑76 (MQ=34)
                  ggCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTa   >  1:2503616/1‑76 (MQ=34)
                  ggCAGGCGGCCCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTa   <  2:2665586/76‑1 (MQ=255)
                   gCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTaa  >  2:121437/1‑76 (MQ=34)
                   |                                                                           
TACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAA  >  NC_000913/731434‑731528

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