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.

GACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCATCCCGCCGCTGATCATCGGCATGACGGTGGTCA  >  NZ_CP009273/3333601‑3333769
                                                                                 |                                                                                       
gACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTacagtgcc                                                                                 <  1:582718‑M2/90‑9 (MQ=255)
  ccGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTacagtgccgt                                                                               >  2:305824‑M2/1‑80 (MQ=255)
            taCCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGtt                                                                     >  1:39968/1‑90 (MQ=255)
                    gaaAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTacagtgccgtagctaaaagcatctaacg                                                             >  1:343496‑M2/1‑62 (MQ=255)
                        ggACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCt                                                         >  1:253125/1‑90 (MQ=255)
                              tAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTacagtgccgtagctaaaagcatctaacgtcctttcttc                                                   >  1:382386‑M2/1‑52 (MQ=255)
                                  tGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGcc                                               >  1:204968/1‑90 (MQ=255)
                                      tttAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAAc                                           >  2:355240/1‑90 (MQ=255)
                                              cGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCa                                   >  1:333672/1‑90 (MQ=255)
                                                        aaTTGTTGGTTTACTTTTGGTCTTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCATCCCGCCGCt                         >  2:247144/1‑90 (MQ=255)
                                                               ggTTTCCTTTTGGTCGTTTACAGTGCCGCCCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCATCCCGCCGCTGatcatc                  <  2:405062/90‑1 (MQ=255)
                                                               ggTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCATCCCGCCGCTGatcatc                  >  2:381164/1‑90 (MQ=255)
                                                                     tttttGGTCGTTTACAGTCCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCATCCCGCCGCTGATCATCGGCATg            <  2:549673/89‑1 (MQ=255)
                                                                      ttttGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCATCCCGCCGCTGATCATCGGCATGa           >  2:611185/1‑90 (MQ=255)
                                                                         tGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCATCCCGCCGCTGATCATCGGCATGACgg        <  2:333672/90‑1 (MQ=255)
                                                                            tCGTTTACAGTGCCGCCCGCTTGGTTTTTGCCGCGTTTATTTTTTGCCGAACTTTTGGCATCCCGCCGCTGATCATCGGCATGACggtgg     <  1:319464/90‑1 (MQ=255)
                                                                            tCGTTTACAGTGCCGCCCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGACCCTTTGGCATCCCGCCGCTGATCATCGGCATGACggtgg     <  2:253125/90‑1 (MQ=255)
                                                                              gTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCATCCCGCCGCTGATCATCGGCATGACGGTGGTc   >  2:372871/1‑90 (MQ=255)
                                                                               tttACAGTGCCGCCCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCATCCCGCCGCTGATCATCGGCATGACGGTGGTCa  <  1:354257/90‑1 (MQ=255)
                                                                               tttACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTTTTCTTTGCCGAACTTTTGGCATCCCGCCGCTGATCATCGGCATGACGGTGGTCa  <  1:577622/90‑1 (MQ=255)
                                                                               tttAAAGTGCCGCCCGCCGGGTTTTTGCCGCGTCTATTCTTTGCCGACCCTTTGGCATCCCGCCGCTGATCATCGGCATGACGGTGGTCa  <  2:204968/90‑1 (MQ=255)
                                                                                 |                                                                                       
GACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGACCGCCTGGTTTTTGCCGCGTCTATTCTTTGCCGAACCTTTGGCATCCCGCCGCTGATCATCGGCATGACGGTGGTCA  >  NZ_CP009273/3333601‑3333769

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