New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 1146416 =45 (0.860)3 (0.060) 3/236 NT 7.1% coding (1653/2589 nt) PP_1006 heme receptor
?NC_002947 1146533 = 35 (0.690)coding (1536/2589 nt) PP_1006 heme receptor

TCGTTGAAGGCGTAGGTCAGGCCGAACTTGGGCGAGACGCGGTGCCACTTCTTGTCCGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCGCTGACCAAGGCCTGCCCGTCACGCGCGCTGTCCTGTCCGATGGCTGTACAG  >  NC_002947/1146267‑1146557
                                                                                                                                                     |                                                                                                                                               
aattggacaaggccttcagcatttgccaaaccgaccacctgtttacctacggcaccacactcacacaggagaaagtcacgggttcacgtagcggcaccgggccttgcctgcacatcggcggtacctgtaccgccatcggacaggacggcg                                                                                                                                                 <  1:911787‑M1/1‑1 (MQ=255)
 cGTTGAAGGCGTAGGTCAGGCCGAACTTGGGCGAGACGCGGTGCCACTTCTTGTCCGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAgc                                                                                                                                                <  1:71706/150‑1 (MQ=255)
  ttggacaaggccttcagcattggccaaaccgaccacctgtttacctacggcaccacactcaaacaggagaaagtcacgggttcacgtagcggcaccggcacttgcctgaacatcggtggtacctgtacagccatcggacaggacagcgcg                                                                                                                                               <  1:204190‑M1/3‑1 (MQ=255)
   tggacaaggccttcagcattggccaaaccgaccacctgtttacctacggcaccacactcaaacaggagaaagtcacgggttcacgtagcggcaccggcacttgcctgaacatcggtggtacctgtacagccatcggacaggacagcgcgc                                                                                                                                              >  2:463815‑M1/147‑150 (MQ=255)
        aaggccttcagcattggccaaaccgaccacctgtttacctacggcaccacactcaaacaggagacagtcacgggttcacgtagcggcacaggcaattgcctgaacatcggtggtacctgtacagccatcggacaggacagcGCGCGCGcc                                                                                                                                         <  2:690652‑M1/9‑1 (MQ=255)
           gTAGGTCAGGCCGAACTTGGGCGAGACGCGGTGCCACTTCTTGTCCGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCggg                                                                                                                                      >  1:618253/1‑150 (MQ=255)
                tCAGGCCGAACTTGGGCGAGACGCGGTGCCACTTCTTGTCCGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGa                                                                                                                                 <  2:1180504/150‑1 (MQ=255)
                       gAACTTGGGCGAGACGCGGTGCCACTTCTTGTCCGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCa                                                                                                                          <  1:1153468/150‑1 (MQ=255)
                       gAACTTGGGCGAGACGCGGTGCCACTTCTTGTCCGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCa                                                                                                                          >  2:142501/1‑150 (MQ=255)
                           ttGGGCGAGACGCGGTGCCACTTCTTGTCCGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAgttg                                                                                                                      >  1:1125296/1‑150 (MQ=255)
                           ttGGGCGAGACGCGGTGCCACTTCTTGTCCGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGAACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGAGCGCGCCGGGCATGAAGGTCCAgttg                                                                                                                      >  1:660882/1‑150 (MQ=255)
                               gCGAGACGCGGTGCCACTTCTTGTCCGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTcc                                                                                                                  >  2:1175348/1‑150 (MQ=255)
                                     cgcgGTGCCACTTCTTGTCCGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGa                                                                                                            <  2:1126241/150‑1 (MQ=255)
                                                ttcttGTCCGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTgc                                                                                                 <  2:1125296/150‑1 (MQ=255)
                                                      tCCGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAAc                                                                                           >  1:131322/1‑150 (MQ=255)
                                                       ccGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACa                                                                                          <  2:131322/150‑1 (MQ=255)
                                                         gAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAgg                                                                                        >  1:97187/1‑150 (MQ=255)
                                                         gAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAgg                                                                                        >  1:662035/1‑150 (MQ=255)
                                                          aaTCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGc                                                                                       <  2:682486/150‑1 (MQ=255)
                                                           aTCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCt                                                                                      >  2:738122/1‑150 (MQ=255)
                                                           aTCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCt                                                                                      <  1:1121051/150‑1 (MQ=255)
                                                               tCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAg                                                                                  >  1:399428/1‑150 (MQ=255)
                                                                       ggTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CAc                                                                        <  1:142501/150‑1 (MQ=255)
                                                                         tAGGCGCCGTACCGCTGGACTCCCGACCGCGCAGGAATTGCTCCGTGAATTTGTGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑TACgg                                                                      >  2:911787/1‑150 (MQ=255)
                                                                             cgcCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCgg                                                                  <  2:618253/150‑1 (MQ=255)
                                                                             cgcCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCgg                                                                  <  1:1175348/150‑1 (MQ=255)
                                                                                     acgtagcggcaccggcacttgcctgaacatcggtggtacctgtacagccatcggacaggacagcGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGa                                                          <  1:89534‑M1/86‑1 (MQ=255)
                                                                                      gCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTgcgc                                                                                               <  1:836524/114‑1 (MQ=255)
                                                                                      gCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTgcgc                                                                                               >  2:836524/1‑114 (MQ=255)
                                                                                      gCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGaa                                                         <  2:97187/150‑1 (MQ=255)
                                                                                        tGGACTCACGACAGCGCAGGAATTCCGCCGTTAATTTGGTCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGt                                                       <  2:662035/150‑1 (MQ=255)
                                                                                          gACTCCAGACCGCGCAGGAATTCATACGTGAATTTGCGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCg                                                     >  2:521357/1‑150 (MQ=255)
                                                                                             tCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCGCTg                                                  >  1:1104432/1‑150 (MQ=255)
                                                                                               cAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAg                                                                                                                         <  1:799620/79‑1 (MQ=255)
                                                                                               cAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAg                                                                                                                         >  2:799620/1‑79 (MQ=255)
                                                                                                    cgcgCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAg                                                                                                                               >  2:536593/1‑68 (MQ=255)
                                                                                                    cgcgCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAg                                                                                                                               <  1:536593/68‑1 (MQ=255)
                                                                                                     gcgcAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCGCTGACCAAGGc                                          <  2:399428/150‑1 (MQ=255)
                                                                                                         aGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAg                                                                                  >  2:954755/1‑108 (MQ=255)
                                                                                                         aGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAg                                                                                  <  1:954755/108‑1 (MQ=255)
                                                                                                         aGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCGCTGACCAAGGCCTGc                                      <  1:443615/150‑1 (MQ=255)
                                                                                                                catcggtggtacctgtacagccatcggacaggacagcGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCGCTGACCAAGGCCTGCCCGTCAc                               >  2:89534‑M1/38‑150 (MQ=255)
                                                                                                                  tcggtggtacctgtacagccatcggacaggacagcGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCGCTGACCAAGGCCTGCCCGTCAcgc                             >  1:690652‑M1/36‑150 (MQ=255)
                                                                                                                    cGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCGCt                                                   >  1:472256/1‑126 (MQ=255)
                                                                                                                    cGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCGCt                                                   <  2:472256/126‑1 (MQ=255)
                                                                                                                     gTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACCGGTGGTGAGTAGT‑‑GACGGTCGGGTCCGGGAAGTCGCTGACCAAGGCCTGCCCGTCAcgcgcg                          >  2:1261963/1‑150 (MQ=255)
                                                                                                                     gTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTc                                                      >  2:894472/1‑122 (MQ=255)
                                                                                                                     gTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTc                                                      <  1:894472/122‑1 (MQ=255)
                                                                                                                            tGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCGCTGACGCTGAAGTGCCCGTCACGCGCGctgtcct                   <  1:1261963/150‑1 (MQ=255)
                                                                                                                                   cATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACA‑‑CCGTTGCCATTACCACTGTTGGGTCCGGGAAGTCGCTGAAAAAGGCCTACCCGTCACGCGCGCTGTCCTGTCCGAt            >  1:738821/1‑150 (MQ=255)
                                                                                                                                       cGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCGCTGACCAAGGCCTGCCCGTCACGCGCGCTGTCCTGTCCGATGGCt        >  2:1207480/1‑150 (MQ=255)
                                                                                                                                           tgtAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCGCTGACCAAGGCCTGCCCGTCACGCGCGCTGTCCTGTCCGATGGCTGTAc    <  2:68351/150‑1 (MQ=255)
                                                                                                                                             tAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCGCTGACCAAGGCCTGCCCGTCACGCGCGCTGTCCTGTCCGATGGCTGTACAg  <  2:1104432/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TCGTTGAAGGCGTAGGTCAGGCCGAACTTGGGCGAGACGCGGTGCCACTTCTTGTCCGAATCGTCCTGGGCGGTAGGCGCCGTACCGCTGGACTCCAGACCGCGCAGGAATTCCTCCGTGAATTTGGGCTCCATGCGGGTGTAGTCGTAGCGCGCGCCGGGCATGAAGGTCCAGTTGTTCCAGCGGATTTCATCCTGCGCGAACAGGCTGTAGGTATT‑‑CACGGTCGGGTCCGGGAAGTCGCTGACCAAGGCCTGCCCGTCACGCGCGCTGTCCTGTCCGATGGCTGTACAG  >  NC_002947/1146267‑1146557

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