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.

ATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTG  >  NZ_CP009273/3816531‑3816707
                                                                                        |                                                                                         
tgctctttataagcccagtgcgcggcaacacccatctccgccatctggtccatatcttcggtacggatctggacctcaaccggcacacCg                                                                                          <  1:307389‑M1/2‑1 (MQ=255)
aTGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGt                                                                                         >  1:370398/1‑90 (MQ=255)
      aTCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGc                                                                                   <  2:458194/90‑1 (MQ=255)
      aTCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGc                                                                                   <  1:426027/90‑1 (MQ=255)
       tCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGcc                                                                                  >  1:136679/1‑90 (MQ=255)
       tCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGcc                                                                                  >  2:359473/1‑90 (MQ=255)
       tCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGcc                                                                                  >  2:283094/1‑90 (MQ=255)
       tCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGcc                                                                                  >  2:202186/1‑90 (MQ=255)
       tCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGcc                                                                                  >  2:40174/1‑90 (MQ=255)
       tCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGcc                                                                                  >  1:541514/1‑90 (MQ=255)
       tCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGcc                                                                                  >  1:539338/1‑90 (MQ=255)
        cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCcg                                                                                 <  2:600891/90‑1 (MQ=255)
               ttCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaaa                                                                          >  2:153812/1‑90 (MQ=255)
               ttCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaaa                                                                          >  2:553304/1‑90 (MQ=255)
               ttCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaaa                                                                          >  2:22153/1‑90 (MQ=255)
               ttCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaaa                                                                          >  2:606848/1‑90 (MQ=255)
                tCCGCGTGACCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAg                                                                         <  1:395900/90‑1 (MQ=255)
                 ccGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGa                                                                        <  2:541506/90‑1 (MQ=255)
                       gATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTatat                                                                  >  2:381703/1‑90 (MQ=255)
                                      tGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGc                                                   >  2:363361/1‑90 (MQ=255)
                                       gACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCg                                                  >  2:469028/1‑90 (MQ=255)
                                           ccTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACg                                              >  1:377739/1‑90 (MQ=255)
                                           ccTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACg                                              <  2:458131/90‑1 (MQ=255)
                                            cTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACgg                                             >  1:165724/1‑90 (MQ=255)
                                            cTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACgg                                             >  2:456504/1‑90 (MQ=255)
                                             tGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGc                                            >  1:429697/1‑90 (MQ=255)
                                                tATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTAt                                         <  2:1182/90‑1 (MQ=255)
                                                  tCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAGAGACTATATCGCCATTCCAAAAGCGGACGGCTATCa                                       <  1:368093/90‑1 (MQ=255)
                                                   cgcgTGCTGGGCCAGATGCACAGCCTGTACAAGCCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCa                                       <  1:349489/90‑1 (MQ=255)
                                                    tatcttcggtacggatctggacctcaaccggcacacCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAg                                      <  1:514687‑M1/54‑1 (MQ=255)
                                                    gcgTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGt                                     <  1:443248/90‑1 (MQ=255)
                                                     cgTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTc                                    >  1:452490/1‑90 (MQ=255)
                                                      gTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGa                                                                        >  1:254797/1‑53 (MQ=255)
                                                      gTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGa                                                                        <  2:254797/53‑1 (MQ=255)
                                                      gTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCt                                   >  1:188368/1‑90 (MQ=255)
                                                        gCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCttt                                 >  1:238230/1‑90 (MQ=255)
                                                           ggtacggatctggacctcaaccggcacacCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGc                               <  2:398739‑M1/61‑1 (MQ=255)
                                                            ggCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGcac                             <  2:165724/90‑1 (MQ=255)
                                                            ggCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGcac                             <  1:40174/90‑1 (MQ=255)
                                                            ggCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGcac                             <  1:469028/90‑1 (MQ=255)
                                                                    ctggacctcaaccggcacacCGTGCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCga                      <  2:494371‑M1/70‑1 (MQ=255)
                                                                         gCCTGTACAAGCCGCGTCCGGGCCGCGCGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCg                >  2:341361/1‑90 (MQ=255)
                                                                             gTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGc                               >  2:640218/1‑71 (MQ=255)
                                                                             gTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGc                               <  1:640218/71‑1 (MQ=255)
                                                                             gTACAAGCCGCATCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGccc            >  1:151725/1‑90 (MQ=255)
                                                                               aCAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGc          <  2:635192/90‑1 (MQ=255)
                                                                                 aaGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCAc        <  2:539338/90‑1 (MQ=255)
                                                                                 aaGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCAc        <  2:541514/90‑1 (MQ=255)
                                                                                 aaGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCAc        <  1:310245/90‑1 (MQ=255)
                                                                                 aaGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCAc        <  1:283094/90‑1 (MQ=255)
                                                                                 aaGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCAc        >  2:343960/1‑90 (MQ=255)
                                                                                  aGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGGACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACg       <  1:550939/90‑1 (MQ=255)
                                                                                   gCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGaa                                                >  1:592230/1‑48 (MQ=255)
                                                                                   gCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGaa                                                <  2:592230/48‑1 (MQ=255)
                                                                                   gCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACgg      >  1:101183/1‑90 (MQ=255)
                                                                                        CGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg  <  2:218476‑M1/90‑1 (MQ=255)
                                                                                        CGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg  <  2:320980‑M1/90‑1 (MQ=255)
                                                                                        CGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg  <  2:554005‑M1/90‑1 (MQ=255)
                                                                                        |                                                                                         
ATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTG  >  NZ_CP009273/3816531‑3816707

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