New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 564995144 (0.780)5 (0.090) 5/254 NT 10.3% coding (1242/2127 nt) PP_4959 diguanylate cyclase
?NC_002947 = 5649985 44 (0.790)coding (1208/2127 nt) PP_4959 diguanylate cyclase

ATCAATCGCGAGCGCGCCCGTCCGTTGAGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTGCTCCGGGTAGGGCAAGCTGTCAAACAGTGCAGTGAACTCATCACCGCCCATACGCGC  >  NC_002947/5649838‑5650133
                                                                                                                                                   |                                                                                                                                                    
aTCAATCGCGAGCGCGCCCGTCCGTTGAGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCc                                                                                                                                                     <  2:54698‑M2/149‑2 (MQ=255)
   aaTCGCGAGCGCGCCCGTCCGTTGAGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGAt                                                                                                                                                  >  3:316484/1‑149 (MQ=255)
    aTCGCGAGCGCGCCCGTCCGTTGAGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGAt                                                                                                                                                  >  6:210256/1‑148 (MQ=255)
         gAGCGCGCCCGTCCGTTGAGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGAc                                                                                                                                            <  4:82923/149‑1 (MQ=255)
             gcgcCCGTCCGTTGAGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCa                                                                                                                                        <  8:78235/149‑1 (MQ=255)
             gcgcCCGTCCGTTGAGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAATCCGGGTAGGTGGCAATCCCGATGCTGTcccc                                                                                                                                         >  3:30392/1‑148 (MQ=255)
                 ccgtccgtTGAGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAgg                                                                                                                                     >  4:158703/1‑148 (MQ=255)
                    tccgtTGAGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGcaca                                                                                                                                 <  4:171636/149‑1 (MQ=255)
                        ttGAGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATc                                                                                                                             <  5:198676/149‑1 (MQ=255)
                        ttGAGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATc                                                                                                                             <  4:316484/149‑1 (MQ=255)
                        ttGAGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATc                                                                                                                             >  3:40298/1‑149 (MQ=255)
                           aGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGa                                                                                                                           <  3:158703/148‑1 (MQ=255)
                           aGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAg                                                                                                                          <  1:54698/149‑1 (MQ=255)
                                 tGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTc                                                                                                                    >  8:302463/1‑149 (MQ=255)
                                   gTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTa                                                                                                                  >  5:350770/1‑149 (MQ=255)
                                          aaGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATg                                                                                                           >  4:116386/1‑149 (MQ=255)
                                                tATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGcc                                                                                                                        >  7:286235/1‑130 (MQ=255)
                                                tATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGcc                                                                                                                        <  8:286235/130‑1 (MQ=255)
                                                  ttGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATg                                                                                                    <  6:350770/148‑1 (MQ=255)
                                                            cccGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTc                                                                                         >  8:99869/1‑149 (MQ=255)
                                                                 ttgcttgGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGc                                                                                     >  1:166040/1‑148 (MQ=255)
                                                                         cTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCg                                                                            <  7:99869/149‑1 (MQ=255)
                                                                            cGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCTAGCCAGTCTACCTCATGGCTGATGGTAACACGTTCGAGCAACCGTTCGGc                                                                          <  3:115336/148‑1 (MQ=255)
                                                                            cGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCgg                                                                           <  2:166040/147‑1 (MQ=255)
                                                                            cGTACATTGCTGCATCGGCCGAACCCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACAGCTAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACAGGTTCGAGCAACCGTTCgg                                                                           <  4:110437/147‑1 (MQ=255)
                                                                               aCATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCaaaa                                                                      >  2:38794/1‑149 (MQ=255)
                                                                                 aTTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGa                                                                                       <  5:265797/130‑1 (MQ=255)
                                                                                 aTTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGa                                                                                       >  6:265797/1‑130 (MQ=255)
                                                                                 aTTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACg                                                                    >  1:303534/1‑149 (MQ=255)
                                                                                     ctgcATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGccc                                                                 <  2:303534/148‑1 (MQ=255)
                                                                                        cATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGcc                                                             <  7:302463/149‑1 (MQ=255)
                                                                                          tCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTg                                                           >  7:84568/1‑149 (MQ=255)
                                                                                          tCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTg                                                           >  5:125421/1‑149 (MQ=255)
                                                                                            ggCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTGCt                                                         <  1:38794/149‑1 (MQ=255)
                                                                                              ccGAACGCAGCATGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCACCCGTTCGGCAAAACGCCCCGCCTGCTcc                                                       >  8:145911/1‑149 (MQ=255)
                                                                                                    gcagcaGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTGCTCCGGGTAg                                                 <  4:40298/149‑1 (MQ=255)
                                                                                                         agGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCcgaaggcctgctgcgttcggccgatgcagcaatgtacgcagccaagcaagcggggcgtcagcaataccgcttctacgaccaggaactcaacggacggg                                                    <  8:288484‑M2/141‑99 (MQ=255)
                                                                                                         agGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCcgaaggcctgctgcgttcggccgatgcagcaatgtacgcagccaagcaagcggggcgtcagcaataccgcttctacgaccaggaactcaacggacggg                                                    >  7:288484‑M2/1‑43 (MQ=255)
                                                                                                            ccTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCcgaaggcctgctgcgttcggccgatgcagcaatgtacgcagccaagcaagcggggcgtcagcaataccgcttctacgaccaggaactcaacggacgggcgcgctcgcga                                         <  5:3579‑M2/149‑110 (MQ=255)
                                                                                                             cTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCcgaaggcctgctgcgttcggccgatgcagcaatgtacgcagccaagcaagcggggcgtcagcaataccgcttctacgaccaggaactcaacggacgggcgcgctcgcgat                                        <  2:102834‑M2/149‑111 (MQ=255)
                                                                                                             cTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTGCTCCGGGTAgg                                                >  6:288732/1‑141 (MQ=255)
                                                                                                             cTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTGCTCCGGGTAgg                                                <  5:288732/141‑1 (MQ=255)
                                                                                                                cGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCcgaaggcctgctgcgttcggccgatgcagcaatgtacgcagccaagcaagcggggcgtcagcaataccgcttctacgaccaggaactcaacggacgggcgcgctcgcgattg                                      >  6:3579‑M2/1‑36 (MQ=255)
                                                                                                                    aTTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACCTCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTGCTCCGGGTAGGGCAAGCTGTCAAAca                                 <  6:291742/149‑1 (MQ=255)
                                                                                                                          gcCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTGCTCCGGGTAGGGCAAGCTGTCAAAcagtgcag                           >  6:40256/1‑149 (MQ=255)
                                                                                                                          gcCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCAGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTGCTCCGGGTAGGGCAAGCTGTCAAAcagtgcag                           >  2:94772/1‑149 (MQ=255)
                                                                                                                            cacaGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCg                                                                >  2:274952/1‑110 (MQ=255)
                                                                                                                            cacaGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCg                                                                <  1:274952/110‑1 (MQ=255)
                                                                                                                               aGTCCGGGTAGGTTGCATTTCCGATGCTGACCCCAAGGCAGACATCGATCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCAGGCCTGCTCCGGGTAGGGCAAGCTGTCAAATAGTGCAGTGAAc                      <  7:145911/149‑1 (MQ=255)
                                                                                                                               aGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTGCTCCGGGTAGGGCAAGCTGTCAAACAGTGCAGTGAAc                      <  5:210256/149‑1 (MQ=255)
                                                                                                                                             gCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTGCTCCGGGTAGGGCAAGCTGTCAAACAGTGCAGTGAACTCATCACCGCCCa         >  8:18827/1‑148 (MQ=255)
                                                                                                                                              cAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTGCTCCGGGTAGGGCAAGCTGTCAAACAGTGCAGTGAACTCATCACCGCCCATa       >  1:42/1‑149 (MQ=255)
                                                                                                                                                   ccGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTGCTCCGGGTAGGGCAAGCTGTCAAACAGTGCAGTGAACTCATCACCGCCCATACgcgc  <  6:24445/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCAATCGCGAGCGCGCCCGTCCGTTGAGTTCCTGGTCGTAGAAGCGGTATTGCTGACGCCCCGCTTGCTTGGCTGCGTACATTGCTGCATCGGCCGAACGCAGCAGGCCTTCGACATTGGCGCCACAGTCCGGGTAGGTGGCAATGCCGATGCTGACCCCAAGGCACACATCGAGCCCGTCTACCTCATGGCTGATGGAAACACGTTCGAGCAACCGTTCGGCAAAACGCCCGGCCTGCTCCGGGTAGGGCAAGCTGTCAAACAGTGCAGTGAACTCATCACCGCCCATACGCGC  >  NC_002947/5649838‑5650133

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