New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	= 4572954	15 (0.650)	4 (0.180)	4/256	NT	19.4%	intergenic (‑89/+112)	hsdM/hsdR	DNA methyltransferase M/pseudogene, premature stop codon, host restriction; endonuclease R
	?	CP009273	= 4573041	19 (0.860)					intergenic (‑176/+25)	hsdM/hsdR	DNA methyltransferase M/pseudogene, premature stop codon, host restriction; endonuclease R

CATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGA  >  CP009273/4572900‑4573184                                                                                                                                              |                                                                                                                                                cATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCCAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACCTGACCGGGGTGAACGAGCGCAGCCAACCCACAGCCAGccc                                                                                                                                              <  1:148567‑M2/145‑4 (MQ=255)            agaagaagaAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAg                                                                                                                                   <  2:136020/145‑1 (MQ=255)             gaagaagaaAGTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAgt                                                                                                                                  <  1:94876/145‑1 (MQ=255)                gaagaaATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGATGTGta                                                                                                                               >  2:235406/1‑145 (MQ=255)                gaagaaATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGAAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGta                                                                                                                               <  2:13271/145‑1 (MQ=255)                  agaaATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGtata                                                                                                                             <  2:28886/145‑1 (MQ=255)                   gaaGTTGCCGAGGGTAATATACACAAACTCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATAc                                                                                                                            <  1:33229/145‑1 (MQ=255)                       ttGCGGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACgtgt                                                                                                                        <  1:389379/145‑1 (MQ=255)                               ggTAATATACACAAAATCATTCAGGTTGCAGCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAg                                                                                                                <  2:114093/145‑1 (MQ=255)                                   atatACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCg                                                                                                            >  1:396448/1‑145 (MQ=255)                                       acacAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGcctgaatgattttgtgtatattaccctcggcaatttcttct                                                                                                        <  1:172136‑M2/145‑42 (MQ=255)                                                 aTTCAGGTTGCATCAAGGCGGCAGGTGAGTGAATCCCCGGGAGCGTACAGAAGAACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGcctgaatgattttgtgtatattaccctcggcaatttcttcttctgcggctc                                                                                              <  2:233382‑M2/145‑52 (MQ=255)                                                   tCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGcctgaatgattttgtgtatattaccctcggcaatttcttcttctgcggctcga                                                                                            <  2:369525‑M2/145‑54 (MQ=255)                                                      ggTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGa                                                                                         <  2:108794/145‑1 (MQ=255)                                                                    ggcAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGACGGGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATGGCTGAATTTGCCCCGCAgg                                                                           >  1:257135/1‑145 (MQ=255)                                                                                      cGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCCGGGTCTCGAGATTATCGTcc                                                         >  1:104741/1‑144 (MQ=255)                                                                                        ggAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCaaa                                                       <  2:62021/145‑1 (MQ=255)                                                                                            cGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAATGtt                                                   <  2:149718/145‑1 (MQ=255)                                                                                                 aGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGATGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTg                                              <  1:235406/145‑1 (MQ=255)                                                                                                               ggtgTGAACGAGCGCAGCCAACTCATATGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTcc                                <  2:305973/142‑1 (MQ=255)                                                                                                                     aaCGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCcgccgc                          <  2:391897/145‑1 (MQ=255)                                                                                                                      aCGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTGCCGCcgcg                         >  1:316339/1‑145 (MQ=255)                                                                                                                                ccAACGCAGAGGCAGcctgaatgattttgtgtatattaccctcggcaatttcttcttctgcggctcgatgaatttgggccgctgcttaatttacggaactcacaatgaacaataacgatctggtcgcgaagctgtggaagctgtg               >  1:369525‑M2/1‑15 (MQ=255)                                                                                                                                     gctgagGCAGCCTTAAGGGTGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAGGTGGcc          >  1:355242/4‑145 (MQ=255)                                                                                                                                             aGCCGGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGa  <  2:396448/145‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                CATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGAATTTGCCCAGCAGGGTATCGAGATTATCGTCAAAGGTTCTTTGCAGCATCGCCTTCCCGCCGCGACGGTGGAAGTTGCCGGTTTTGA  >  CP009273/4572900‑4573184

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