New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4507779 =52 (0.880)4 (0.070) 4/246 NT 7.9% coding (1002/1392 nt) cysG multifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase
?NC_002947 4507834 = 45 (0.810)coding (947/1392 nt) cysG multifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase

CATGTAGAACACCAAGGTCTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAA  >  NC_002947/4507632‑4507925
                                                                                                                                                    |                                                                                                                                                  
atgtaggcaagcgccgcgctgatcattcagtgccgcaggaccagatcaaccgcctgttggtcgatctggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagag                                                                                                                                                    <  7:2582‑M1/1‑1 (MQ=255)
  gtaggcaagcgccgcgctgatcattcagtgccgcaggaccagatcaaccgcctgttggtcgatctggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCt                                                                                                                                                  <  1:187716‑M1/3‑1 (MQ=255)
   taggcaagcgccgcgctgatcattcagtgccgcaggaccagatcaaccgcctgttggtcgatctggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCTg                                                                                                                                                 <  2:159013‑M1/4‑1 (MQ=255)
     ggcaagcgccgcgctgatcattcagtgccgcaggaccagatcaaccgcctgttggtcgatctggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCTGgc                                                                                                                                               <  2:21672‑M1/6‑1 (MQ=255)
      gcaagcgccgcgctgatcattcagtgccgcaggaccagatcaaccgcctgttggtcgatctggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCTGgcc                                                                                                                                              <  6:61362‑M1/7‑1 (MQ=255)
           ccAAGGTCTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGt                                                                                                                                         >  5:323808/1‑148 (MQ=255)
            cAAGGTCTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTaa                                                                                                                                       >  6:242371/1‑149 (MQ=255)
            cAAGGTCTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTaa                                                                                                                                       >  7:9908/1‑149 (MQ=255)
            cAAGGTCTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTaa                                                                                                                                       >  7:334195/1‑149 (MQ=255)
                  cTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTg                                                                                                                                 <  5:39119/149‑1 (MQ=255)
                     cgcCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTgg                                                                                                                                <  7:269370/147‑1 (MQ=255)
                        cGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTCGGCATAGTCGCGGTGAGACAGCCCAATGCCGCCATACGCGGAGCAGCCGCAGGCCCGCGTAATGCCTGGcaaca                                                                                                                           >  8:160803/1‑149 (MQ=255)
                            gcCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTg                                                                                                                       <  8:334195/149‑1 (MQ=255)
                              caccaGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTgg                                                                                                                      >  2:311234/1‑148 (MQ=255)
                              caccaGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTgg                                                                                                                      >  3:336197/1‑148 (MQ=255)
                                 caGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACg                                                                                                                  >  3:358784/1‑149 (MQ=255)
                                     tCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAAt                                                                                                              <  6:323808/149‑1 (MQ=255)
                                       aTTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGc                                                                                                            >  1:41611/1‑149 (MQ=255)
                                         tCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGcc                                                                                                           >  1:152008/1‑148 (MQ=255)
                                         tCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCAGGAACGGAATGcc                                                                                                           >  4:132477/1‑148 (MQ=255)
                                           gatcaaccgcctgttggtcgatctggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCAt                                                                                                         >  8:2582‑M1/106‑149 (MQ=255)
                                            aaGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGc                                                                                                       >  4:78096/1‑149 (MQ=255)
                                             aGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCt                                                                                                      <  1:311234/149‑1 (MQ=255)
                                               gCAGGTTGCTGGTGCCGTCCTTGAGGTGGCCCGGTGACGAAGCGCACCGACTGGGAATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGATGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGATc                                                                                                     >  1:21672/1‑149 (MQ=255)
                                                cAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGcc                                                                                                           <  5:108624/141‑1 (MQ=255)
                                                cAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGcc                                                                                                           >  6:108624/1‑141 (MQ=255)
                                                 aGGTTGCTGGTGCCGTCCTTTAGGTGCCCTGTGACGAAGCGCAGCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGc                                                                                                  <  7:160803/149‑1 (MQ=255)
                                                   gTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGc                                                                                                  >  4:337704/1‑147 (MQ=255)
                                                        ttggtcgatctggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGc                                                                                              >  1:58162‑M1/93‑147 (MQ=255)
                                                          gTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTc                                                                                         <  4:24155/149‑1 (MQ=255)
                                                          gTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTc                                                                                         <  3:337704/149‑1 (MQ=255)
                                                               gTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGAtct                                                                                    >  7:284389/1‑149 (MQ=255)
                                                               gTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGAtct                                                                                    >  3:236623/1‑149 (MQ=255)
                                                                  tggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTc                                                                                         <  5:379276‑M1/60‑1 (MQ=255)
                                                                  tggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTc                                                                                         >  6:379276‑M1/83‑142 (MQ=255)
                                                                   cccgcgTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTc                                                                                <  1:132097/144‑1 (MQ=255)
                                                                        gTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTcgccgc                                                                           >  5:192393/1‑149 (MQ=255)
                                                                            cccGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGcccc                                                                        >  3:42249/1‑149 (MQ=255)
                                                                              cGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCgg                                                                      >  6:114477/1‑149 (MQ=255)
                                                                              cGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGAGCTCTTCGCCGCCCCgg                                                                      <  3:132477/149‑1 (MQ=255)
                                                                                   aCGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCTATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCg                                                                   >  7:73246/1‑147 (MQ=255)
                                                                                     gAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCg                                                                   >  1:175765/1‑145 (MQ=255)
                                                                                          gcACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAAc                                                          <  5:114477/149‑1 (MQ=255)
                                                                                               gACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTc                                                     <  3:113405/149‑1 (MQ=255)
                                                                                                 cTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCTATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGc                                                   <  8:73246/149‑1 (MQ=255)
                                                                                                     gCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCAcc                                               >  4:33973/1‑149 (MQ=255)
                                                                                                       aTAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCt                                             >  5:336357/1‑149 (MQ=255)
                                                                                                         aGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCtt                                            >  5:61362/1‑148 (MQ=255)
                                                                                                          gTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACTGGTCGCCACCCTTCa                                          <  4:311933/149‑1 (MQ=255)
                                                                                                              cgGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGgcg                                      <  8:284389/149‑1 (MQ=255)
                                                                                                                gTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCgca                                    <  8:9908/149‑1 (MQ=255)
                                                                                                                     tCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCAcc                               >  3:153417/1‑149 (MQ=255)
                                                                                                                       aGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCg                             >  1:242011/1‑149 (MQ=255)
                                                                                                                        gCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGc                            <  2:152008/149‑1 (MQ=255)
                                                                                                                            aaTGCCGCCATACGCGGAGCAGCCGCTGACCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGc                        >  1:210673/1‑149 (MQ=255)
                                                                                                                                 cgccATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGc                   <  5:242371/149‑1 (MQ=255)
                                                                                                                                      tACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGc              <  3:33973/149‑1 (MQ=255)
                                                                                                                                        cgcgGAGCAGCCGCTGGCCGCCGTAATGCATGGCACCACCTAGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGcc             >  2:187716/1‑148 (MQ=255)
                                                                                                                                         gcgGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGc              >  3:21634/1‑146 (MQ=255)
                                                                                                                                           ggAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGc              >  7:33509/1‑144 (MQ=255)
                                                                                                                                            gAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATc        <  3:78096/149‑1 (MQ=255)
                                                                                                                                                  ccGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACAGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGCCCaa  >  4:225992/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATGTAGAACACCAAGGTCTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAA  >  NC_002947/4507632‑4507925

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