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.

TACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTG  >  NC_000913/731434‑731535
                           |                                                                          
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                       <  1:5603222/76‑1 (MQ=35)
     ctacagcACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACcaca                       >  2:2873903/1‑76 (MQ=34)
     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:3779856/76‑1 (MQ=35)
         agcACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGAcc                   >  1:4955542/1‑76 (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)
                    cAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACcac                        <  2:704950/60‑1 (MQ=35)
                    cAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACcac                        >  1:704950/1‑60 (MQ=38)
                     cggCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAAcc       <  1:121437/75‑1 (MQ=25)
                     aGGCTGACCGGCGTTCACCCCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAAcc       <  2:2088930/76‑1 (MQ=25)
                     aGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAAcc       <  2:4604689/76‑1 (MQ=35)
                      ggCTGACCGGCGTTCCCACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCg      <  2:2487476/76‑1 (MQ=25)
                       gCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGc     <  1:4373801/76‑1 (MQ=35)
                       gCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGc     >  2:2764296/1‑76 (MQ=34)
                       gCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGc     >  1:2106818/1‑76 (MQ=34)
                       gCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGc     <  1:4761440/76‑1 (MQ=35)
                         tGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGcc  >  2:3733247/1‑75 (MQ=21)
                         tGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGcc  >  1:924815/1‑75 (MQ=25)
                          gACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCCTCCCGTATACACCAGACCCGGCGGGGGACAGCCTg  >  1:5650814/1‑76 (MQ=11)
                          gACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTg  <  1:3531167/76‑1 (MQ=35)
                          gACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTg  >  1:3386074/1‑76 (MQ=35)
                          gACCGCCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTg  <  1:1533876/76‑1 (MQ=25)
                           |                                                                          
TACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATACCACAGACCCGGCAGGTAACCGCCTG  >  NC_000913/731434‑731535

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