New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5866803 =	17 (1.020)	3 (0.190)	3/272	1.8	14.6%	coding (1326/2253 nt)	cadA‑III	cadmium translocating P‑type ATPase
	?	NC_002947	5866845 =	19 (1.200)					coding (1284/2253 nt)	cadA‑III	cadmium translocating P‑type ATPase
Rejected: Frequency below/above cutoff threshold.

ACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATGGTCACCGGGGTCGAGATCACCAGCGCGCACGGGCAGGCCACCACCAGC  >  NC_002947/5866680‑5866935                                                                                                                            |                                                                                                                                      aCGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTa                                                                                                            <  1:152882/150‑1 (MQ=255)        actgggtttaccgcgccttggtgctgctggtggtggcctgcccgtgcgcgctggtgatctcgaccccggtgaccatcgtcagtggcctggccgccgccgcgcgcaaaggcatcctgATCAAGGGTGGCGGCCACCCTCCAGGTACACGcc                                                                                                      >  2:9377‑M1/117‑150 (MQ=255)                          gCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCt                                                                                   >  2:266104/1‑150 (MQ=255)                            ttGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCttt                                                                                 <  1:243048/150‑1 (MQ=255)                            ttGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCttt                                                                                 <  2:347393/150‑1 (MQ=255)                                   cGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGcgcgcg                                                                          <  1:207032/150‑1 (MQ=255)                                     gTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCgcggc                                                                        <  1:240809/150‑1 (MQ=255)                                                gcggcgCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCggccaggc                                                             <  2:104447/150‑1 (MQ=255)                                                           aCCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATg                                                  <  1:93649/150‑1 (MQ=255)                                                                   tcgaccccggtgaccatcgtcagtggcctggccgccgccgcgcgcaaaggcatcctgATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCtt                                                                                                >  2:113399‑M1/58‑96 (MQ=255)                                                                   tcgaccccggtgaccatcgtcagtggcctggccgccgccgcgcgcaaaggcatcctgATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCtt                                                                                                <  1:113399‑M1/39‑1 (MQ=255)                                                                      ggTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATGGTCACCGGGGt                                       <  2:173441/150‑1 (MQ=255)                                                                         cTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATGGTCACCGGGGTCga                                    <  1:79810/150‑1 (MQ=255)                                                                              ccGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATGGTCACCGGGGTCGAGATCa                               <  1:3999/150‑1 (MQ=255)                                                                                         tGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATGGTCACCGGGGTCGAGATCACCAGCGCGCAc                    >  2:309967/1‑150 (MQ=255)                                                                                         tGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATGGTCACCGGGGTCGAGATCACCAGCGCGCAc                    <  2:356491/150‑1 (MQ=255)                                                                                               aTGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATGGTCACCGGGGTCGAGATCACCAGCGCGCACGGGCAg              <  1:280838/150‑1 (MQ=255)                                                                                                       ggTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATGGTCACCGGGGTCGAGATCACCAGCGCGCACGGGCAGGccaccac      >  2:163270/1‑150 (MQ=255)                                                                                                       ggTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATGGTCACCGGGGTCGAGATCACCAGCGCGCACGGGCAGGccaccac      <  1:266104/150‑1 (MQ=255)                                                                                                           ttGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATGGTCACCGGGGTCGAGATCACCAGCGCGCACGGGCAGGCCACCACcagc  <  1:177994/150‑1 (MQ=255)                                                                                                                            |                                                                                                                                      ACGCTCGCCCAGGCTGGCCGCCAAGGCTTGCGCACGGTCTTCGAACAGCGGCGCCAGCACCTTGGCGTCGGTCTGCACCGGCTTGCCATGGGTGATGGTGCCGGTCTTGTCCAGGGCCAGGAAATCAAGG‑TGGCGGCCACCCTCCAGGTACACGCCACCCTTGATCAGGATGCCTTTGCGCGCGGCGGCGGCCAGGCCACTGACGATGGTCACCGGGGTCGAGATCACCAGCGCGCACGGGCAGGCCACCACCAGC  >  NC_002947/5866680‑5866935

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