New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP009273 = 457295415 (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.