New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1340634 =16 (0.770)4 (0.210) 3/154 1.7 24.4% coding (366/1986 nt) pdeR cyclic di‑GMP phosphodiesterase
?NZ_CP009273 1340740 = 10 (0.520)coding (260/1986 nt) pdeR cyclic di‑GMP phosphodiesterase

TGACGTCGTGTTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCA  >  NZ_CP009273/1340546‑1340709
                                                                                        |                                                                           
tGACGTCGTGTTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCgg                                                                            >  1:58299/1‑90 (MQ=255)
   tcatcacctccagcctgatgatgagtctgaccgtcgatgaaagcgatctttctgtgcatctggtaggacgaaaaatcaataaacgGGAAt                                                                         >  1:37529‑M1/86‑90 (MQ=255)
    gtcgtGTTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAAtg                                                                        >  1:13659/1‑90 (MQ=255)
            ccagcctgatgatgagtctgaccgtcgatgaaagcgatctttctgtgcatctggtaggacgaaaaatcaataaacgGGAATGTGCTTCAg                                                                <  2:349969‑M1/14‑1 (MQ=255)
            ccagcctgatgatgagtctgaccgtcgatgaaagcgatctttctgtgcatctggtaggacgaaaaatcaataaacgGGAATGTGCTTCAg                                                                <  2:200243‑M1/14‑1 (MQ=255)
                     cTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCt                                                       >  2:596475/1‑90 (MQ=255)
                      tgtgTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTg                                                      >  2:21177/1‑90 (MQ=255)
                       gtgtAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGc                                                     >  1:572486/1‑90 (MQ=255)
                       gtgtAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGc                                                     <  1:79443/90‑1 (MQ=255)
                         gtAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTc                                                   <  1:21177/90‑1 (MQ=255)
                            tctgaccgtcgatgaaagcgatctttctgtgcatctggtaggacgaaaaatcaataaacgGGAATGTGCTTCAGAAACTACCTGCTCAGc                                                >  1:349969‑M1/61‑90 (MQ=255)
                             tCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCa                                               <  1:596475/90‑1 (MQ=255)
                             ctgaccgtcgatgaaagcgatctttctgtgcatctggtaggacgaaaaatcaataaacgGGAATGTGCTTCAGAAACTACCTGCTCAGCa                                               >  1:200243‑M1/60‑90 (MQ=255)
                                 cacaTAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAgg                                           >  2:232993/1‑90 (MQ=255)
                                      aaCCGATTGAAGCGTTGGATATTGCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGc                                      >  2:284392/1‑90 (MQ=255)
                                          gATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCgt                                  >  2:248799/1‑90 (MQ=255)
                                          gATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCgt                                  >  1:75055/1‑90 (MQ=255)
                                                             cccccGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAAc               >  1:98052/1‑90 (MQ=255)
                                                                ccGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAAccgc            >  1:389988/1‑90 (MQ=255)
                                                                      gTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGgtgt      >  2:421795/1‑90 (MQ=255)
                                                                          agaATCACTATTGCGAAATGTGCTTCAAAACCTCCCGGCTCAGCAAAGAAAAGCCCGTGTGATAAATCACGAGAAACCGCCGGTGTGTca  <  2:65024/90‑1 (MQ=255)
                                                                                        |                                                                           
TGACGTCGTGTTCTTTCAACCCTGTGTAATCTTCACATAACCGATTGAAGCGTTGGATATTCCCCCGGCTGTCGAGAATCACTATTGCGGAATGTGCTTCAGAAACTACCTGCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCA  >  NZ_CP009273/1340546‑1340709

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