New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 270718820 (0.920)3 (0.160) 3/148 1.6 13.2% coding (934/1335 nt) srmB ATP‑dependent RNA helicase SrmB
?NZ_CP009273 = 2707267 22 (1.160)coding (1013/1335 nt) srmB ATP‑dependent RNA helicase SrmB
Rejected: Frequency below/above cutoff threshold.

GACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGC  >  NZ_CP009273/2707191‑2707349
                                                                            |                                                                                  
gACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGcc                                                                       <  2:142285/90‑1 (MQ=255)
     tCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCg                                                                  >  2:402238/1‑90 (MQ=255)
              cAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCAc                                                         >  1:8198/1‑90 (MQ=255)
              cAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCAc                                                         >  1:470859/1‑90 (MQ=255)
                             cTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGtaccgcgcgcggcaacatcggttgcgaccagtacgtttacgc                                          >  1:47693‑M2/1‑48 (MQ=255)
                              ttCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCtggtg                                         >  2:242/1‑90 (MQ=255)
                               tCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCtggtgg                                        >  2:144161/1‑90 (MQ=255)
                                cGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGa                                       >  2:316834/1‑90 (MQ=255)
                                cGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGa                                       <  1:312330/90‑1 (MQ=255)
                                  atatGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAg                                     <  2:470859/90‑1 (MQ=255)
                                  atatGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAg                                     <  1:181125/90‑1 (MQ=255)
                                   tatGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGtaccgcgcgcggcaacatcggttgcgaccagtacgtttacgcgacctt                                    >  1:83687‑M2/1‑42 (MQ=255)
                                       ccGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCAt                                <  1:616076/90‑1 (MQ=255)
                                         gcgcAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGa                              >  2:596627/1‑90 (MQ=255)
                                          cgcAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGAc                             <  2:246052/90‑1 (MQ=255)
                                                ggCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGGGGAAGCCCATGACCATctg                       >  1:509054/1‑90 (MQ=255)
                                                    aTACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATctgctgc                   <  2:96834/90‑1 (MQ=255)
                                                    aTACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATctgctgc                   <  1:208546/90‑1 (MQ=255)
                                                          aTTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTa             >  1:73940/1‑90 (MQ=255)
                                                                aCCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAg       >  2:199768/1‑90 (MQ=255)
                                                                  cGTATCGGACGtaccgcgcgcggcaacatcggttgcgaccagtacgtttacgcgaccttcggtcaaacgcttgatcgcttcgttacgctt     >  2:356570‑M2/1‑11 (MQ=255)
                                                                   gTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGcc    <  1:391204/90‑1 (MQ=255)
                                                                   gTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGcc    <  1:376815/90‑1 (MQ=255)
                                                                     aTCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGc  <  1:596627/90‑1 (MQ=255)
                                                                     aTCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGc  <  2:73940/90‑1 (MQ=255)
                                                                            |                                                                                  
GACATTCCTGACGTCAGCCACGTCTTTAACTTCGATATGCCGCGCAGTGGCGATACTTATTTGCACCGTATCGGACGTACCGCGCGCGCCGGTCGTAAAGGCACCGCAATTTCGCTGGTGGAAGCCCATGACCATCTGCTGCTGGGTAAAGTAGGCCGC  >  NZ_CP009273/2707191‑2707349

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.