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

TGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGA  >  CP009273/4572820‑4573097
                                                                                                                                      |                                                                                                                                               
tGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTGGTTGGTCAGTGTGAGTGCCGTAGAGTAAGCAGCGGCCCAAATTGATCGAGGCGCAGAAGAAGAAATGGCCGAGGGGAAGAGACACAAACTCATTCAGGTGGCagccag                                                                                                                                       >  1:374465/1‑140 (MQ=255)
     cacagcttccacagcttcGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGctgcctctgcgttgg                                                                                                                                  >  1:369525‑M1/1‑130 (MQ=255)
                               aTCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTAc                                                                                                        >  2:50297/1‑145 (MQ=255)
                                                        gTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGacaga                                                                               >  1:246728/1‑141 (MQ=255)
                                                          aaaTTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAgc                                                                             <  2:284790/145‑1 (MQ=255)
                                                                       gCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCag                                                                <  2:272552/145‑1 (MQ=255)
                                                                        cccAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAACCATTCAGGTTCCATCAAGGCGGCAAGTGAGTGAATCCCCAGGAGCGTACAGAAGTACGTGACCGGGGAGAACCAGCGCAGCCAACGCaga                                                               >  2:363620/1‑145 (MQ=255)
                                                                          cAACTTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAgg                                                             <  1:309708/145‑1 (MQ=255)
                                                                                cATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCCAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACCTGACCGGGGTGAACGAGCGCAGCCAACCCACAGCCAGccc                                                       <  1:148567‑M2/145‑4 (MQ=255)
                                                                                  tCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGctgcctctgcgttggctgcgctcgttcaccccggtcacgtacttctgtacgctcccggggattcactcacttgccgccttgatgcaacctga                                                     <  2:369525‑M1/145‑93 (MQ=255)
                                                                                    gAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGctgcctctgcgttggctgcgctcgttcaccccggtcacgttcttctgtacgctcccggggattcactcacctgccgccttgatgcaacctgaat                                                   <  2:233382‑M1/145‑95 (MQ=255)
                                                                                           agaagaagaAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAg                                            <  2:136020/145‑1 (MQ=255)
                                                                                            gaagaagaaAGTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAgt                                           <  1:94876/145‑1 (MQ=255)
                                                                                              agaagaaATTGCCGAGGGTAATATACACAAAATCATTCAGGctgcctctgcgttggctgcgctcgttcaccccggtcacgtacttctgtacgctcccggggattcactcacttgccgccttgatgcaacctgaatgattttgtgt                                         <  1:172136‑M1/145‑105 (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)
                                                                                                                                    gggctggctgtgggttggctgcgctcgttcaccccggtcaggtacttctgtacgctcccggggattcactcacttgccgcctggatgcaacctgaatgattttgtgtatattaccctcggcaatttcttcttctgcggctcgatg   <  1:148567‑M1/144‑143 (MQ=255)
                                                                                                                                     ggTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGa  <  2:108794/145‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                               
TGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCATTCAGGTTGCATCAAGGCGGCAAGTGAGTGAATCCCCGGGAGCGTACAGAAGTACGTGACCGGGGTGAACGAGCGCAGCCAACGCAGAGGCAGCCTGAAGGATGAAGTGTATACGTGTCAGGCCAGCTCGTCCCAGATATAATCGCTGA  >  CP009273/4572820‑4573097

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.