New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	4368615 =	16 (0.920)	4 (0.250)	3/284	NT	20.5%	coding (147/1134 nt)	ampC	penicillin‑binding protein; beta‑lactamase, intrinsically weak
	?	CP009273	4368733 =	16 (0.990)					coding (29/1134 nt)	ampC	penicillin‑binding protein; beta‑lactamase, intrinsically weak

CCGGCACCTGCAATGGCAGGCCGCCAGCAGTGTAGGTTGCGAGATGTAATAGTGTGATCCCATTCCACTGTTTAGCGGTAAGTTCAGGCCAGTATTTTGTTGTGGGATCGCTTAACTTGATTTCCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGT  >  CP009273/4368347‑4368879                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                         ccGGCACCTGCAATGGCAGGCCGCCAGCAGTGTAGGTTGCGAGATGTAATAGTGTGATCCCATTCCACTGTTTAGCGGTAAGTTCAGGCCAGTATTTTGTTGTGGGATCGCTTAACTTGATTTCCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCANTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTAACCCTGATAAATTac                                                                                                                                                                                                                                                                     >  1:310252/1‑274 (MQ=255)      aCCTGCAATGGCAGGCCGCCAGCAGTGTAGGTTGCGAGATGTAATAGTGTGATCCCATTCCACTGTTTAGCGGTAAGTTCAGGCCAGTATTTTGTTGTGGGATCGCTTAACTTGATTTCCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTTCTTTTTGGCGATTTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgccac                                                                                                                                                                                                                                                               >  1:87565/1‑275 (MQ=255)                                                               atCCATTGTTTAGCGGTAAGTTCAGGCCAGTATTTTGTTGTGGGATCGCTTAACTTGATTTCCCCGCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACaa                                                                                                                                                                                                       <  1:56739/273‑1 (MQ=255)                                                                                                       tgGGATCGCTTAACTTGATTTCCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATaa                                                                                                                                                                                                                               >  2:64627/1‑210 (MQ=255)                                                                                                           aTCGCTTAACTTGATTTCCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgg                                                                                                                                                                                                                           <  1:64629/210‑1 (MQ=255)                                                                                                                                                                                                  cGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATATACAtt                                                                   >  1:5251/1‑275 (MQ=255)                                                                                                                                                                                                                      tctTTTTGGCAATGTCCGCACAGCCCCAGTTAAAGTTATAAGTTTTACCCTTATCCATTTCCGCCGCCGCCATACCCGTGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGTTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTaa                                               <  2:5251/274‑1 (MQ=255)                                                                                                                                                                                                                                    aaaccagaccctatgttcaaaacgacgctctgcgccttattAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGctc                                                                                                                                                                                                                                    <  1:203559‑M1/39‑1 (MQ=255)                                                                                                                                                                                                                                        cagaccctatgttcaaaacgacgctctgcgccttattAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                                                                                                                                >  2:203553‑M1/38‑80 (MQ=255)                                                                                                                                                                                                                                                        aacgacactctgcgccttattAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGctc                                                                                                                                                                                                                                    <  1:122766‑M1/39‑1 (MQ=255)                                                                                                                                                                                                                                                            taaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                                                                                                                                     >  2:164191/1‑71 (MQ=255)                                                                                                                                                                                                                                                            taaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                                                                                                                                               >  2:319704/1‑125 (MQ=255)                                                                                                                                                                                                                                                            acgctctgcgccttattAATTACCGCCACCGCCATACCCGTGATCTTTTTTTTCACtata                                                                                                                                                                                                                                >  2:122764‑M1/18‑60 (MQ=255)                                                                                                                                                                                                                                                                gTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgtg                                                                                                                                                                                                                 <  1:164195/71‑1 (MQ=255)                                                                                                                                                                                                                                                                gTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCTCAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                                                                                                                                           <  1:319715/125‑1 (MQ=255)                                                                                                                                                                                                                                                                   tACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcag                                                                                                                                                                               >  2:140758/1‑102 (MQ=255)                                                                                                                                                                                                                                                                   tACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAATTGTTGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAAAATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCTATGCGTTAGATCTTAACGACACCACTCAGCGAGACAACTGTCCGTATACCAGCCCGACCga  >  2:21717/1‑274 (MQ=255)                                                                                                                                                                                                                                                                   tACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                                                                                                                                               >  2:84676/1‑118 (MQ=255)                                                                                                                                                                                                                                                                       cTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaaa                                                                                                                                                                           <  1:140761/102‑1 (MQ=255)                                                                                                                                                                                                                                                                       cTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                                                                                                                                           <  1:84678/118‑1 (MQ=255)                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                         CCGGCACCTGCAATGGCAGGCCGCCAGCAGTGTAGGTTGCGAGATGTAATAGTGTGATCCCATTCCACTGTTTAGCGGTAAGTTCAGGCCAGTATTTTGTTGTGGGATCGCTTAACTTGATTTCCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGT  >  CP009273/4368347‑4368879

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