New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3816618 =51 (2.200)3 (0.130) 3/164 1.9 6.7% coding (859/2109 nt) spoT bifunctional GTP diphosphokinase/guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
?NZ_CP009273 3816703 = 34 (1.520)coding (944/2109 nt) spoT bifunctional GTP diphosphokinase/guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Rejected: Frequency below/above cutoff threshold.

GCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCA  >  NZ_CP009273/3816614‑3816790
                                                                                         |                                                                                       
gCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACgg                                                                                         >  1:101183/1‑90 (MQ=255)
    CGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg                                                                                     <  2:554005‑M1/90‑1 (MQ=255)
    CGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg                                                                                     <  2:320980‑M1/90‑1 (MQ=255)
    CGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg                                                                                     <  2:218476‑M1/90‑1 (MQ=255)
      TCCGGGCCGCGTGAAAGACTATATCGCCACTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGcc                                                                                   >  2:386027‑M1/1‑90 (MQ=255)
         gggCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGt                                                                                >  1:359639/1‑90 (MQ=255)
          ggCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGtt                                                                               <  2:541307/90‑1 (MQ=255)
           gCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTg                                                                              <  2:489226/90‑1 (MQ=255)
           gCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTg                                                                              <  1:405049/90‑1 (MQ=255)
            ccGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGa                                                                             <  2:43541/90‑1 (MQ=255)
            ccGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGa                                                                             <  1:606848/90‑1 (MQ=255)
                 tGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTcc                                                                        <  1:359473/90‑1 (MQ=255)
                   tcgaggtgtgcaaagactgatagccgttcgcttttggaatggcgatatagtctttcacgcggcccgcacgGTGTGCCGGTTGAGGTCCAg                                                                      <  2:494371‑M2/20‑1 (MQ=255)
                        cTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCg                                                                 >  2:85352/1‑90 (MQ=255)
                            gcaaagactgatagccgttcgcttttggaatggcgatatagtctttcacgcggcccggacgGTGTGCCGGTTGAGGTCCAGATCCGTAcc                                                             <  2:398739‑M2/29‑1 (MQ=255)
                              cGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGa                                                           <  1:343960/90‑1 (MQ=255)
                                ccATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAg                                                         >  2:457610/1‑90 (MQ=255)
                                ccATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAg                                                         >  2:224059/1‑90 (MQ=255)
                                  aTTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGat                                                       >  1:423135/1‑90 (MQ=255)
                                   ctgatagccgttcgcttttggaatggcgatatagtctttcacgcggcccggacgGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGata                                                      <  1:514687‑M2/36‑1 (MQ=255)
                                             gAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGAt                                            <  2:101183/90‑1 (MQ=255)
                                             gAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGAt                                            >  1:531207/1‑90 (MQ=255)
                                                          tCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGa                                                           <  1:119004/62‑1 (MQ=255)
                                                          tCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGa                                                           >  2:119004/1‑62 (MQ=255)
                                                            tttGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGa                                                                             <  2:490714/42‑1 (MQ=255)
                                                            tttGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGa                                                                             >  1:490714/1‑42 (MQ=255)
                                                            tttGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGgtgt                             >  2:319636/1‑90 (MQ=255)
                                                            tttGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGgtgt                             <  2:359639/90‑1 (MQ=255)
                                                               gCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATgggc                                                 >  2:358919/1‑65 (MQ=255)
                                                               gCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAc                                                 <  1:358919/67‑1 (MQ=255)
                                                                   acCTCGATGATCGGCCCGCACGGTGTGCCGGTTGa                                                                             >  2:124866/1‑35 (MQ=255)
                                                                   acCTCGATGATCGGCCCGCACGGTGTGCCGGTTGa                                                                             <  1:124866/35‑1 (MQ=255)
                                                                       cGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACt                  <  2:429697/90‑1 (MQ=255)
                                                                       cGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACt                  <  2:435563/90‑1 (MQ=255)
                                                                       cGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACt                  <  1:22153/90‑1 (MQ=255)
                                                                       cGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACt                  <  1:153307/90‑1 (MQ=255)
                                                                       cGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACt                  <  1:569777/90‑1 (MQ=255)
                                                                         atgatCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACTgg                <  1:589009/90‑1 (MQ=255)
                                                                              cGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTta           >  2:441020/1‑90 (MQ=255)
                                                                                gCCCGCACGGTGTGCCGGTTGATGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTtata         <  2:577849/90‑1 (MQ=255)
                                                                                   cGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAag      >  2:307389/1‑90 (MQ=255)
                                                                                       cgGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCa  <  1:307389‑M2/88‑1 (MQ=255)
                                                                                         |                                                                                       
GCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCA  >  NZ_CP009273/3816614‑3816790

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