New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 521169044 (0.830)3 (0.060) 3/238 NT 6.3% coding (260/1638 nt) PP_4592 metal dependent hydrolase
?NC_002947 = 5211711 48 (0.970)coding (239/1638 nt) PP_4592 metal dependent hydrolase

GCACACCCAGCAGCAGGATATACGCGGCAAACTTGAACGACATCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGACTCAGCCACAACCCCGAAAGGCCGGCAGGAACCACCGCTTGCAGG  >  NC_002947/5211565‑5211857
                                                                                                                                                   |                                                                                                                                                  
gCACACCCAGCAGCAGGATATACGCGGCAAACTTGAACGACATCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACt                                                                                                                                                   <  4:157941‑M2/148‑2 (MQ=255)
      ccAGCAGCAGGATATACGCGGCAAACTTGAACGACATCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACtcatgaa                                                                                                                                             <  5:72803‑M2/148‑8 (MQ=255)
         gcagcagGATATACGCGGCAAACTTGAACGACATCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAACGCCACGCCCGCGCCACTCATGaacaa                                                                                                                                          <  1:3066/149‑1 (MQ=255)
           agcagGATATACGCGGCAAACTTGAACGACATCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAacac                                                                                                                                       <  4:125599/149‑1 (MQ=255)
                 atatACGCGGCAAACTTGAACGACATCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAg                                                                                                                                 <  7:200895/149‑1 (MQ=255)
                    taCGCGGCAAACTTGAACGACATCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCAACGCCCGCGCCACTCATGAACAACACGATAAGTAc                                                                                                                              >  6:218153/1‑149 (MQ=255)
                        cgGCAAACTTGAACGACATCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACtc                                                                                                                                                  <  2:8063‑M2/125‑3 (MQ=255)
                          gCAAACTTGAACGACATCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCt                                                                                                                         >  6:220892/1‑148 (MQ=255)
                             aaCTTGAACGACATCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCAc                                                                                                                     <  6:72803/149‑1 (MQ=255)
                              aattgaACGACATCAAGTCGCGAACCTACGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGGCCGCGCCACtcatg                                                                                                                                               <  5:304146‑M2/120‑6 (MQ=255)
                                         aTCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCgg                                                                                                           >  3:143718/1‑147 (MQ=255)
                                         aTCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCgg                                                                                                           >  3:204241/1‑147 (MQ=255)
                                          tCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAg                                                                                                        >  8:165697/1‑149 (MQ=255)
                                          tCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAg                                                                                                        >  7:267310/1‑149 (MQ=255)
                                            aaGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAgg                                                                                                       >  3:20881/1‑148 (MQ=255)
                                                     aCCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTGACGt                                                                                             >  4:78916/1‑149 (MQ=255)
                                                     aCCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAAc                                                                                               <  5:220892/147‑1 (MQ=255)
                                                       cTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGt                                                                                             >  5:238220/1‑147 (MQ=255)
                                                       cTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGa                                                                                           <  4:143718/149‑1 (MQ=255)
                                                       cTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGa                                                                                           <  4:204241/149‑1 (MQ=255)
                                                           gCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAgg                                                                                         <  5:262525/147‑1 (MQ=255)
                                                           gCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAgg                                                                                         >  6:262525/1‑147 (MQ=255)
                                                                    gTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGaa                                                                              >  5:261234/1‑149 (MQ=255)
                                                                     tctcGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAAc                                                                             <  4:20881/149‑1 (MQ=255)
                                                                         ggCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAAcgc                                                                           <  3:157941/147‑1 (MQ=255)
                                                                            cATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAaga                                                                      <  8:267310/149‑1 (MQ=255)
                                                                              tatTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCaa                                                                                                    >  5:130665/1‑117 (MQ=255)
                                                                              tatTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCaa                                                                                                    <  6:130665/117‑1 (MQ=255)
                                                                              tatTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGa                                                                                   <  7:105853/134‑1 (MQ=255)
                                                                              tatTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGa                                                                                   >  8:105853/1‑134 (MQ=255)
                                                                                ttACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCAt                                                                  <  6:238220/149‑1 (MQ=255)
                                                                                        cATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATcagc                                                              <  4:192087/145‑1 (MQ=255)
                                                                                        cATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATcagc                                                              >  3:192087/1‑145 (MQ=255)
                                                                                           gCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCa                                                                                                     <  4:234708/103‑1 (MQ=255)
                                                                                           gCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCa                                                                                                     >  3:234708/1‑103 (MQ=255)
                                                                                            ccGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGa                                                                               <  7:97712/124‑1 (MQ=255)
                                                                                            ccGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGa                                                                               >  8:97712/1‑124 (MQ=255)
                                                                                                  tCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGc                                                   >  7:333512/1‑146 (MQ=255)
                                                                                                    aaTAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGAc                                              >  1:22545/1‑149 (MQ=255)
                                                                                                    aaTAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGAc                                              <  8:333512/149‑1 (MQ=255)
                                                                                                         gTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTGACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGACTCAGc                                         <  3:78916/149‑1 (MQ=255)
                                                                                                            cGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACTATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAACATCAGCACAGCAAACGCCAGACTCAGCcac                                      <  7:20180/149‑1 (MQ=255)
                                                                                                                 tACTGGTTCATGAAGTAAGCCACGCCCGCGCCACtcatgaaccagtacggcgtacttattgacgccggcatgttccgtaatatggccgagaccaatattgcggaggtgcgcgacttgatgtcgttcaagtttgccgcgtatatcctgct                                 <  6:108113‑M2/149‑116 (MQ=255)
                                                                                                                    tGGTTCATGAAGTAAGCCACGCCCGCGCCACtcatgaaccagtacggcgtacttattgacgccggcatgttccgtaatatggccgagaccaatattgcggaggtgcgcgacttgatgtcgttcaagtttgccgcgtatatcctg                                   >  5:296757‑M2/1‑31 (MQ=255)
                                                                                                                    tGGTTCATGAAGTAAGCCACGCCCGCGCCACtcatgaaccagtacggcgtacttattgacgccggcatgttccgtaatatggccgagaccaatattgcggaggtgcgcgacttgatgtcgttcaagtttgccgcgtatatcctg                                   <  6:296757‑M2/144‑114 (MQ=255)
                                                                                                                      gTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATc                                                                 >  7:210275/1‑112 (MQ=255)
                                                                                                                      gTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATc                                                                 <  8:210275/112‑1 (MQ=255)
                                                                                                                      gTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTCAACGCGAAGAGCATCAGCACAGCAAACGCCAGACTCAGCCACAACCCCGa                              >  5:212279/1‑147 (MQ=255)
                                                                                                                       ttCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGACTCAGCCACAACCCCGAAAg                           <  7:165697/149‑1 (MQ=255)
                                                                                                                             aaGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGACTCAGCCACAAcc                                  >  5:29220/1‑136 (MQ=255)
                                                                                                                                 aaGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGACTCAGCCACAACCCCGAAAGGc                         <  8:10061/141‑1 (MQ=255)
                                                                                                                                 aaGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGACTCAGCCACAACCCCGAAAGGc                         >  7:10061/1‑141 (MQ=255)
                                                                                                                                    gCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGACTCAGCCACAACCCCGAAAGGCCGGCAGGAac               <  2:22545/149‑1 (MQ=255)
                                                                                                                                     ccACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGACTCAGCCACAACCCCGAAAGGCCGGCAGGAacc              >  2:9130/1‑149 (MQ=255)
                                                                                                                                       aCGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGACTCAGCCACAACCCCGAAAGGCCGGCAGGAaccac            <  5:218153/149‑1 (MQ=255)
                                                                                                                                       aCGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGACTCAGCCACAACCCCGAAAGGCCGGCAGGAaccac            >  5:39283/1‑149 (MQ=255)
                                                                                                                                           ccGCGGCACTGATGAACAACACGATATGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCAt                                                                  <  2:91350/91‑1 (MQ=255)
                                                                                                                                                 cACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGACTCAGCCACAACCCCGAAAGGCCGGCAGGCACCCCCGCTAGCAgg  >  4:256482/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCACACCCAGCAGCAGGATATACGCGGCAAACTTGAACGACATCAAGTCGCGCACCTCCGCAATATTGGTCTCGGCCATATTACGGAACATGCCGGCGTCAATAAGTACGCCGTACTGGTTCATGAAGTAAGCCACGCCCGCGCCACTCATGAACAACACGATAAGTACCGGCTTCACTACATACCGGAAGGCAATTAACGTGAGGATCAGATTGAACGCGAAGAGCATCAGCACAGCAAACGCCAGACTCAGCCACAACCCCGAAAGGCCGGCAGGAACCACCGCTTGCAGG  >  NC_002947/5211565‑5211857

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.