New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1899984 =39 (0.540)5 (0.100) 3/46 3.7 10.9% coding (1558/1599 nt) pdeD putative c‑di‑GMP phosphodiesterase PdeD
?NC_000913 1900017 = 55 (1.120)coding (1591/1599 nt) pdeD putative c‑di‑GMP phosphodiesterase PdeD
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTT  >  NC_000913/1899950‑1900018
                                  |                                  
cAAATGGCTGGCGGGCAGCCAACCGCCgcccgccc                                    >  1:1298799/1‑35 (MQ=255)
    tggctggcGGGCAGCCAACCGCCGCCCGCCCGGCa                                <  1:231582/35‑1 (MQ=255)
    tggctggcGGGCAGCCAACCGCCGCCCGCCCGGCa                                <  1:4696054/35‑1 (MQ=255)
      gctggcGGGCAGCCAACCGCCGCCCGCCCGGCATa                              <  1:7799469/35‑1 (MQ=255)
       ctggcGGGCAGCCAGCCGCCGCCCGCCCGGCATaa                             <  1:1307865/35‑1 (MQ=255)
       ctggcGGGCAGCCAACCGCCGCCCGCCCGGCATaa                             <  1:1194730/35‑1 (MQ=255)
       ctggcGGGCAGCCAACCGCCGCCCGCCCGGCATaa                             <  1:6471720/35‑1 (MQ=255)
       ctggcGGGCAGCCAACCGCCGCCCGCCCGGCATaa                             <  1:4151338/35‑1 (MQ=255)
        tggcGGGCAGCCAACCGCCGCCCGCCCGGCATAAt                            <  1:86908/35‑1 (MQ=255)
        tggcGGGCAGCCAACCGCCGCCCGCCCGGCATAAt                            <  1:3088587/35‑1 (MQ=255)
          gcGGGCAGCCAACCGCCGCCCGCCCGGCATAATgg                          >  1:2457904/1‑35 (MQ=255)
           cGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGa                         <  1:247474/35‑1 (MQ=255)
           cGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGa                         <  1:4384508/35‑1 (MQ=255)
             ggCAGCCAACCGCCGCCCGCCCGGCATAAAGGGCa                       >  1:4663413/1‑35 (MQ=37)
               cAGCCAACCGCCGCCCGCCCGGCATAATGGACata                     <  1:2730325/35‑1 (MQ=255)
                aGCCAACCGCCGCCCGCCCGGCATAATGGACatat                    >  1:5302723/1‑35 (MQ=255)
                aGCCAACCGCCGCCCGCCCGGCATAATGGACatat                    <  1:3359255/35‑1 (MQ=255)
                 gCCAACCGCCGCCCGCCCGGCATAATGGACATATc                   <  1:336101/35‑1 (MQ=255)
                 gCCAACCGCCGCCCGCCCGGCATAATGGACATATc                   <  1:5562266/35‑1 (MQ=255)
                 gCCAACCGCCGCCCGCCCGGCATAATGGACATATc                   <  1:7283018/35‑1 (MQ=255)
                  ccAACCGCCGCCCGCCCGGCATAATGGACATATCa                  <  1:3133445/35‑1 (MQ=255)
                  ccAACCGCCGCCCGCCCGGCATAATGGACATATCa                  <  1:1628436/35‑1 (MQ=255)
                    aaCCGCCGCCCGCCCGGCATAATGGACATATCACg                <  1:3975250/35‑1 (MQ=255)
                     aCCGCCGCCCGCCCGGCATAATGGACATATCACGc               >  1:5178444/1‑35 (MQ=255)
                       cgccgcCCGCCCGGCATAATGGACATATCACGccc             <  1:1433733/35‑1 (MQ=255)
                         ccgcccgcCCGGCATAATGGACATATCACGCCCAt           <  1:6094150/35‑1 (MQ=255)
                         ccgcccgcCCGGCATAATGGACATATCACGCCCAt           <  1:3066341/35‑1 (MQ=255)
                          cgcccgccCGGCATAATGGACATATCACGCCCAtt          <  1:4136578/35‑1 (MQ=255)
                           gcccgcccGGCATAATGGACATATCCCGCCCATTa         <  1:478352/35‑1 (MQ=255)
                            cccgcccgGCATAATGGACATATCACGCCCATTAt        <  1:344111/35‑1 (MQ=255)
                             ccgcccgGCATAATGGACATATCACGCCCATTATg       <  1:205032/35‑1 (MQ=255)
                             ccgcccgGCATAATGGACATATCACGCCCATTATg       <  1:4636122/35‑1 (MQ=255)
                             ccgcccgGCATAATGGACATATCACGCCCATTATg       <  1:241710/35‑1 (MQ=255)
                              cgcccgGCATAATGGACATATCACGCCCATTATGc      <  1:5981279/35‑1 (MQ=255)
                               gcccgGCATAATGGACATATCACGCCCATTATGcc     <  1:383371/35‑1 (MQ=255)
                                cccgGCATAATGGACATATCACGCCCATTATGCcg    >  1:5103086/1‑35 (MQ=255)
                                cccgGCATAATGGACATATCACGCCCATTATGCcg    <  1:6905883/35‑1 (MQ=255)
                                 ccgGCATAATGGACATATCACGCCCATTATGCcgt   <  1:4181813/35‑1 (MQ=255)
                                 ccgGCATAATGGACATATCACGCCCATTATGCcgt   <  1:7754315/35‑1 (MQ=255)
                                 acgGCATAATGGACATATCACGCCCATTATGCcgt   >  1:715104/2‑35 (MQ=255)
                                 acgGCATAATGGACATATCACGCCCATTATGCcgt   >  1:2896141/2‑35 (MQ=255)
                                  cgGCATAATGGACATATCACGCCCATTATGCcgtt  >  1:4292626/1‑35 (MQ=255)
                                  cgGCATAATGGACATATCACGCCCATTATGCcgta  >  1:36305/1‑34 (MQ=255)
                                  |                                  
CAAATGGCTGGCGGGCAGCCAACCGCCGCCCGCCCGGCATAATGGACATATCACGCCCATTATGCCGTT  >  NC_000913/1899950‑1900018

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