New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4049856 =33 (1.040)3 (0.110) 3/148 2.4 10.0% intergenic (‑157/+129) glnL/glnA nitrogen regulation protein NR(II)/glutamate‑‑ammonia ligase
?NZ_CP009273 4049890 = 25 (0.890)intergenic (‑191/+95) glnL/glnA nitrogen regulation protein NR(II)/glutamate‑‑ammonia ligase
Rejected: Frequency below/above cutoff threshold.

AGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACT  >  NZ_CP009273/4049804‑4049977
                                                                                      |                                                                                       
cctgggatcggtggcaagcacatcacgccggaaacgacgccaatgcgacttaTCCGGCCTACCCCATCATCACCGTGTAGGCCGGATAAg                                                                                      <  2:800917‑M1/38‑1 (MQ=255)
  tgggatcggtggcaagcacatcacgccggatgcgacgcaaatgcgtcttatccggcctacacggtgatgatgtggtaggccggaTAAGAc                                                                                    <  2:906187‑M2/6‑1 (MQ=255)
    gggCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGc                                                                                  >  2:210771/1‑90 (MQ=255)
     ggCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCa                                                                                 >  1:725168/1‑90 (MQ=255)
      gCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCAt                                                                                >  1:199364/1‑90 (MQ=255)
      gCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCAt                                                                                >  2:577271/1‑90 (MQ=255)
              tAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCg                                                                        >  2:483848/1‑90 (MQ=255)
              tAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCg                                                                        >  2:15132/1‑90 (MQ=255)
                 cTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCAt                                                                     >  1:613692/1‑90 (MQ=255)
                   gACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATcc                                                                   >  2:401309/1‑90 (MQ=255)
                     catcacgccggatgcgacgcaaatgcgtcttttcggcctacacggtgatgatgtggtaggccggaTAAGACGCATTTGCGTCGCATCCgg                                                                 <  2:515184‑M2/25‑1 (MQ=255)
                          ttCACGTTGGAGAGCGACTCACCTGCTCCGGACTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGa                                                            >  1:529377/1‑90 (MQ=255)
                           ccggatgcgacgcaaatgcgtcttaTCCGGCCTACCACATCATCACCGTGTAGGCCGAAAAGACGCATTTGCGTCGCATCCGGCGTGAtg                                                          <  2:515184‑M1/65‑1 (MQ=255)
                                 tGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCtt                                                     >  2:229202/1‑90 (MQ=255)
                                       gCGACTCACCTGCTCCGGCCTACCACATCATCCCCGTGTAGGCCGGTTAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCAcc                                               <  2:624434/90‑1 (MQ=255)
                                       gCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCAcc                                               <  1:53076/90‑1 (MQ=255)
                                        cGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCg                                              <  2:128859/90‑1 (MQ=255)
                                        cGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCg                                              <  1:577271/90‑1 (MQ=255)
                                        cGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCg                                              <  2:455156/90‑1 (MQ=255)
                                              gtcttaTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCa                                        <  2:906187‑M1/84‑1 (MQ=255)
                                              aCCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCa                                        <  2:199364/90‑1 (MQ=255)
                                                cttatccggcctacacggtgatgatggggtaggccggaTAAGTCGCATTGGCGTCGTTTCCGGCGTGATGTGCTTGCCACCGATCCCAgg                                      <  2:800917‑M2/52‑1 (MQ=255)
                                                  gCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGcc                                    <  2:516623/90‑1 (MQ=255)
                                                  gCTCCGGCCTACCACATCATCACCGTGTAGCCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGcc                                    <  1:854064/90‑1 (MQ=255)
                                                       GGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCa                               <  1:151166‑M1/90‑1 (MQ=255)
                                                       GGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCa                               <  2:792945‑M1/90‑1 (MQ=255)
                                                       GGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCa                               <  2:655484‑M1/90‑1 (MQ=255)
                                                       GGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCa                               <  1:282661‑M1/90‑1 (MQ=255)
                                                       GGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCa                               <  1:846130‑M1/90‑1 (MQ=255)
                                                            ACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGAc                          <  2:354779‑M1/90‑1 (MQ=255)
                                                            ACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGAc                          <  1:912893‑M1/90‑1 (MQ=255)
                                                                    atcaCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCAcc                                               >  2:701489/1‑61 (MQ=255)
                                                                    atcaCCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCAcc                                               <  1:701489/61‑1 (MQ=255)
                                                                          gtgtAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTc            >  2:570635/1‑90 (MQ=255)
                                                                          gtgtAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTc            >  1:362169/1‑90 (MQ=255)
                                                                          gtgtAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTc            >  1:251704/1‑90 (MQ=255)
                                                                          gtgtAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTc            >  1:836275/1‑90 (MQ=255)
                                                                          gtgtAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTc            >  2:894271/1‑90 (MQ=255)
                                                                              aggccggaTAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACg        >  1:442953‑M2/9‑90 (MQ=255)
                                                                               ggccggaTAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACgg       >  2:272320‑M2/8‑90 (MQ=255)
                                                                               ggccggaTAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACgg       <  1:570635‑M2/83‑1 (MQ=255)
                                                                                 ccggaTAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCa     <  2:818753‑M2/85‑1 (MQ=255)
                                                                                 ccggaTAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCa     >  2:86686‑M2/6‑90 (MQ=255)
                                                                                    gaTAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACt  >  2:107978‑M2/3‑90 (MQ=255)
                                                                                    gaTAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACt  >  1:443427‑M2/3‑90 (MQ=255)
                                                                                    gaTAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACt  >  2:825228‑M2/3‑90 (MQ=255)
                                                                                      |                                                                                       
AGATGGGCTACAGATAGCTGACAAACTTCACGTTGGAGAGCGACTCACCTGCTCCGGCCTACCACATCATCACCGTGTAGGCCGGATAAGACGCATTTGCGTCGCATCCGGCGTGATGTGCTTGCCACCGATCCCAGGCCTGCCAGAGACAGGCGAAAAGTTTCCACGGCAACT  >  NZ_CP009273/4049804‑4049977

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