New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 3415693 =19 (0.660)7 (0.260) 7/248 2.4 27.5% coding (1591/1830 nt) PP_3029 DNA cytosine methyltransferase family protein
?NC_002947 3415734 = 19 (0.700)coding (1550/1830 nt) PP_3029 DNA cytosine methyltransferase family protein

GATCAGGCGGATCTGGTCGGTCTTATTGATGGGCCGCCATTCCAGATGGCCGGTTTCCTGGTTCTCGAACAGGCCGCGATCAAGGATGTAGTCGGGGCTGAAGCCTTGAGCGAGCTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAGGTACTTGTGCAGGAAGTGAGCGCACTTGCGGGCCTTCACCAACAACTCCGGCGGCAGGATCGCTGCGGGCACCTGAACAACGGTGA  >  NC_002947/3415544‑3415842
                                                                                                                                                     |                                                                                                                                                     
gATCAGGCGGATCTGGTCGGTCTTATTGATGGGCCGCCATTCCAGATGGCCGGTTTCCTGGTTCTCGAACAGGCCGCGATCAAGGATGTAGTCGGGGCTGAAGCCTTGAGCGAGCTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGaa                                                                                                                                                       <  2:404154/150‑1 (MQ=255)
   atgcacaccctgcagaccaagaagcgcattgcgctggtcaccgttgttcatgtgaccgcagcgatcctgccgccggagtggttggtgaaggcccgcaagtgcgctcacttcctgcacaagtacctgccggagcacttccccgatctAGTc                                                                                                                                                    >  2:347105‑M1/147‑150 (MQ=255)
      gCGGATCTGGTCGGTCTTATTGATGGGCCGCCATTCCAGATGGCCGGTTTCCTGGTTCTCGAACAGGCCGCGATCAAGGATGTAGTCGGGGCTGAAGCCTTGAGCGAGCTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGAc                                                                                                                                                 >  2:130976/1‑150 (MQ=255)
                           gATGGGCCGCCATTCCAGATGGCCGGTTTCCTGGTTCTCGAACAGGCCGCGATCAAGGATGTAGTCGGGGCTGAAGCCTTGAGCGAGCTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGagcag                                                                                                                            <  1:130976/150‑1 (MQ=255)
                           cgcatggcgctggtcaccgttgttcaggtgcccgcagcgatcctgccgccggagttgttggtgaaggcccgcaagtgcgctcacttcctgcacaagtacctgccggagcacttccccgatctAGTCGACCAGTGCGTGCTCGCCGagcag                                                                                                                            <  2:677190‑M1/28‑1 (MQ=255)
                                 gcgctggtcaccgttgttcaggtgcccgcagcgatcctgccgccggagttgttggtgaaggcccgcaagtgcgctcacttcctgcacaagtacctgccggagcacttccccgatctAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAg                                                                                                                      >  1:677190‑M1/117‑150 (MQ=255)
                                      ggtcaccgttgttcaggtgcccgcagcgatcctgccgccggagttgttggtgaaggcccgcaagtgcgctcacttcctgcacaagtacctgccggagcacttccccgatctAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGa                                                                                                                 >  1:719421‑M1/112‑150 (MQ=255)
                                       ttCCAGATGGCCGGTTTCCTGGTTCTCGAACAGGCCGCGATCAAGGATGTAGTCGGGGCTGAAGCCTTGAGCGAGCTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGAc                                                                                                                <  1:239016/150‑1 (MQ=255)
                                               ggCCGGTTTCCTGGTTCTCGAACAGGCCGCGATCAAGGATGTAGTCGGGGCTGAAGCCTTGAGCGAGCTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCgg                                                                                                        <  1:466058/150‑1 (MQ=255)
                                                  cGGTTTCCTGGTTCTCGAACAGGCCGCGATCAAGGATGTAGTCGGGGCTGAAGCCTTGAGCGAGCTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGa                                                                                                     >  1:147722/1‑150 (MQ=255)
                                                               ctcGAACAGGCCGCGATCAAGGATGTAGTCGGGGCTGAAGCCTTGAGCGAGCTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAg                                                                                        <  1:298139/150‑1 (MQ=255)
                                                                   aaCAGGCCGCGATCAAGGATGTAGTCGGGGCTGAAGCCTTGAGCGAGCTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGAc                                                                                                                                                 >  1:192009/1‑89 (MQ=255)
                                                                   aaCAGGCCGCGATCAAGGATGTAGTCGGGGCTGAAGCCTTGAGCGAGCTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGAc                                                                                                                                                 <  2:192009/89‑1 (MQ=255)
                                                                                             ggggCTGAAGCCTTGAGCGAGCTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAGGTACTTGTGCAGGAAGTGAGCGCACTTGCg                                                          >  1:370757/1‑150 (MQ=255)
                                                                                             gcccgcaagtgcgctcacttcctgcacaagtacctgccggagcacttccccgatctAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTc                                                                                                                   <  1:391625‑M1/37‑1 (MQ=255)
                                                                                             gcccgcaagtgcgctcacttcctgcacaagtacctgccggagcacttccccgatctAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTc                                                                                                                   >  2:391625‑M1/57‑93 (MQ=255)
                                                                                                    aaGCCTTGAGCGAGCTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAGGTACTTGTGCAGGAAGTGAGCGCACTTGCGGGCCTTc                                                   >  1:202469/1‑150 (MQ=255)
                                                                                                             gcgagcTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAGGTACTTGTGCAGGAAGTGAGCGCACTTGCGGGCCTTCACCaacaac                                          >  2:4434/1‑150 (MQ=255)
                                                                                                                       gTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAGGTACTTGTGCAGGAAGTGAGCGCACTTGCGGGCCTTCACCAACAACTCCGGCGGCa                                <  2:147722/150‑1 (MQ=255)
                                                                                                                          gtacctgccggagcacttccccgatctAGTCGACCAGTGCGTGCTCGCCGagcagc                                                                                                                           >  2:379409‑M1/28‑56 (MQ=255)
                                                                                                                          gtacctgccggagcacttccccgatctAGTCGACCAGTGCGTGCTCGCCGagcagc                                                                                                                           <  1:379409‑M1/29‑1 (MQ=255)
                                                                                                                            ggggCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAGGTACTTGTGCAGGAAGTGAGCGCACTTGCGGGCCTTCACCAACAACTCCGGCGGCAGGATc                           <  1:4434/150‑1 (MQ=255)
                                                                                                                              ggCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAGGTACTTGTGCAGGAAGTGAGCGCACTTGCgg                                                         <  1:453561/118‑1 (MQ=255)
                                                                                                                              ggCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAGGTACTTGTGCAGGAAGTGAGCGCACTTGCgg                                                         >  2:453561/1‑118 (MQ=255)
                                                                                                                                          gCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAg                                                                                        <  1:670948/75‑1 (MQ=255)
                                                                                                                                          gCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAg                                                                                        >  2:670948/1‑75 (MQ=255)
                                                                                                                                                  tGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAGGTACTTGTGCAGGAAGTGAGCGCACTTGCGGGCCTTCACCAACAACTCCGGCGGCAGGATCGCTGCGGGCACCTGAACAACgg     >  2:194608/1‑150 (MQ=255)
                                                                                                                                                     AGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAGGTACTTGTGCAGGAAGTGAGCGCACTTGCGGGCCTTCACCAACAACTCCGGCGGCAGGATCGCTGCGGGCACCTGAACAACGGTGa  >  1:521363‑M1/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATCAGGCGGATCTGGTCGGTCTTATTGATGGGCCGCCATTCCAGATGGCCGGTTTCCTGGTTCTCGAACAGGCCGCGATCAAGGATGTAGTCGGGGCTGAAGCCTTGAGCGAGCTTCAGTTCCGGGGCCTTGAGCATGCGGAGGGTGAAGTCGACCAGTGCGTGCTCGCCGAGCAGCACCAGGTCGACTAGATCGGGGAAGTGCTCCGGCAGGTACTTGTGCAGGAAGTGAGCGCACTTGCGGGCCTTCACCAACAACTCCGGCGGCAGGATCGCTGCGGGCACCTGAACAACGGTGA  >  NC_002947/3415544‑3415842

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