New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 33336559 (0.410)3 (0.150) 3/152 1.6 19.3% coding (22/978 nt) yrbG calcium/sodium antiporter
?NZ_CP009273 = 3333682 17 (0.870)coding (49/978 nt) yrbG calcium/sodium antiporter
Rejected: Frequency below/above cutoff threshold.

TTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCAT  >  NZ_CP009273/3333582‑3333737
                                                                         |                                                                                  
tttCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACt                                                                    >  2:577622/1‑90 (MQ=255)
   cTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACtttt                                                                 <  1:528012/90‑1 (MQ=255)
        aaTCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCg                                                            <  1:529793/90‑1 (MQ=255)
          tCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGtt                                                          <  1:162443/90‑1 (MQ=255)
                   gACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTacagtgcc                                                 <  1:582718‑M2/90‑9 (MQ=255)
                     ccGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTacagtgccgt                                               >  2:305824‑M2/1‑80 (MQ=255)
                               taCCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGtt                                     >  1:39968/1‑90 (MQ=255)
                                    gaagaaAGGACGTTAGATGCTTTTAGCTACGGCACTGTaaacgaccaaaagtaaaccaacaattaacagtgccgtagctaaaagcatcta                                >  1:382386‑M1/1‑38 (MQ=255)
                                       gaaAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTacagtgccgtagctaaaagcatctaacg                             >  1:343496‑M2/1‑62 (MQ=255)
                                           ggACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCt                         >  1:253125/1‑90 (MQ=255)
                                              cGTTAGATGCTTTTAGCTACGGCACTGTaaacgaccaaaagtaaaccaacaattaacagtgccgtagctaaaagcatctaacgtcctttc                      >  1:343496‑M1/1‑28 (MQ=255)
                                                 tAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTacagtgccgtagctaaaagcatctaacgtcctttcttc                   >  1:382386‑M2/1‑52 (MQ=255)
                                                     tGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGcc               >  1:204968/1‑90 (MQ=255)
                                                         tttAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAAc           >  2:355240/1‑90 (MQ=255)
                                                                aCGGCACTGTaaacgaccaaaagtaaaccaacaattaacagtgccgtagctaaaagcatctaacgtcctttcttcaggtatactcgccgg    >  2:305824‑M1/1‑10 (MQ=255)
                                                                 cGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCa   >  1:333672/1‑90 (MQ=255)
                                                                  ggCACTGTaaacgaccaaaagtaaaccaacaattaacagtgccgtagctaaaagcatctaacgtcctttcttcaggtatactcgccggtc  <  1:582718‑M1/90‑83 (MQ=255)
                                                                         |                                                                                  
TTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCAT  >  NZ_CP009273/3333582‑3333737

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