New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 158231 =65 (0.290)27 (0.120) 2/198 36.0 0.36 intergenic (+34/‑196) parC/tuf1 DNA topoisomerase IV subunit A/translation elongation factor Tu
?NC_014034 1253896 = 30 (0.130)intergenic (+54/+37) RCAP_rcc01181/hisH hypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Rejected: Position hash score below cutoff.

GCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAGGCG  >  NC_014034/158133‑158327
                                                                                                  |                                                                                                
gCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCgggg                                                                                                <  2:1589798/101‑1 (MQ=255)
 cAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCggggg                                                                                               <  1:4117725/101‑1 (MQ=255)
     cGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTcc                                                                                           <  1:1441535/101‑1 (MQ=255)
      gggCGAGCTCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCGGCCGCGCCCGACCCGCGGGGGCTCCg                                                                                          <  2:1512441/101‑1 (MQ=255)
      gggCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTcct                                                                                          <  1:2834904/101‑2 (MQ=255)
      gggCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCg                                                                                          <  1:3376142/101‑1 (MQ=255)
      gggCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCCCCCGACCCCCggcgccctgc                                                                                          >  2:3819392/1‑93 (MQ=255)
          gggcgcCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGcccc                                                                                      <  1:2033084/99‑1 (MQ=255)
          gggcgcCGGGCGCATGGCGCCGCGCGGTTTCCCGCGCGACATCCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGcccc                                                                                      <  2:4181207/99‑1 (MQ=255)
          gAGCGCCGGGCGGATGGCTCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGcccc                                                                                      <  2:484575/101‑1 (MQ=255)
          gAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGcccc                                                                                      <  2:2184386/101‑1 (MQ=255)
          gAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGcccc                                                                                      <  2:353666/101‑1 (MQ=255)
           agggcCGGGCGGATGGCGACGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccccc                                                                                     <  2:2653834/98‑1 (MQ=255)
           aGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccccc                                                                                     <  1:720864/101‑1 (MQ=255)
           aGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccccc                                                                                     <  1:2608649/101‑1 (MQ=255)
            gcgcCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCg                                                                                    <  2:3805921/101‑1 (MQ=255)
            gcgcCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCg                                                                                    <  2:3025028/101‑1 (MQ=255)
            gcgcCGGGCGGATGGCGCCCCGCGATTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCg                                                                                    <  2:495442/101‑1 (MQ=255)
             cgcCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGc                                                                                   <  2:2249550/101‑1 (MQ=255)
              gcCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCa                                                                                  <  2:1139184/101‑1 (MQ=255)
              gcCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCa                                                                                  <  1:4313415/101‑1 (MQ=255)
              gcCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCa                                                                                  <  2:663020/101‑1 (MQ=255)
                cGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAcc                                                                                <  1:955762/101‑1 (MQ=255)
                   gCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAccccc                                                                             <  1:3691672/101‑1 (MQ=255)
                   gCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAccccc                                                                             <  2:175079/101‑1 (MQ=255)
                   gCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAccccc                                                                             <  1:2741757/101‑1 (MQ=255)
                   gCGGAGGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAccccc                                                                             <  1:576453/101‑1 (MQ=255)
                    cggcgggCGTCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCg                                                                            <  2:1679148/96‑1 (MQ=255)
                        tGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGt                                                                        <  1:3362715/101‑1 (MQ=255)
                         ggCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTa                                                                       <  2:3853540/101‑1 (MQ=255)
                          gcgcCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAt                                                                      <  1:4085890/101‑1 (MQ=255)
                           cgccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtt                                                                     <  2:2897878/101‑1 (MQ=255)
                           cgccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtt                                                                     <  1:1976542/101‑1 (MQ=255)
                             cgaggcgGTTTCGCGCGAAACAACCGGTTCCGCTGACGGGGTTGCGGGCGCTCCCGCCCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtttt                                                                   <  1:1057569/97‑1 (MQ=25)
                              cgcgcgGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTg                                                                  <  1:2991045/101‑1 (MQ=255)
                                      ttCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaa                                                          <  1:2491918/101‑1 (MQ=255)
                                      ttCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaa                                                          <  1:3942385/101‑1 (MQ=255)
                                        cccGCCCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaaaa                                                        <  2:125873/101‑1 (MQ=255)
                                         ccGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAg                                                       <  1:2764256/101‑1 (MQ=255)
                                           gcgcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCa                                                     <  2:172691/101‑1 (MQ=255)
                                           gcgcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCa                                                     <  1:4377326/101‑1 (MQ=255)
                                            cgcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCaa                                                    <  2:151094/101‑1 (MQ=255)
                                              cgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAgg                                                  <  1:2000441/101‑1 (MQ=255)
                                                   aCCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGccc                                             <  2:1665133/101‑1 (MQ=255)
                                                   aCCGGTCCACCTGAGGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGCGCCCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGccc                                             <  1:4423768/101‑1 (MQ=37)
                                                     cGGTCCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGccccc                                           <  1:625271/101‑1 (MQ=255)
                                                       gTTCCCCTGACGGGGTTTCGGGCGCTGCCGCTCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCg                                         <  1:1046618/101‑1 (MQ=255)
                                                       gTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCg                                         <  2:1914900/101‑1 (MQ=255)
                                                        ttCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCgg                                        <  1:3835057/101‑1 (MQ=255)
                                                           cccTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACg                                     <  1:1545111/101‑1 (MQ=255)
                                                              tGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgtgt                                  <  2:2169496/101‑1 (MQ=255)
                                                                aCGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCg                                <  1:2692984/101‑1 (MQ=255)
                                                                  ggggTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGtt                              <  2:1710604/101‑1 (MQ=255)
                                                                    ggTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTc                            <  1:4347960/101‑1 (MQ=255)
                                                                    ggTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTc                            <  2:3917462/101‑1 (MQ=255)
                                                                      tttCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCgg                          <  1:1790807/101‑1 (MQ=255)
                                                                      tttCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCgg                          <  1:918412/101‑1 (MQ=255)
                                                                      tttCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCgg                          <  1:1283717/101‑1 (MQ=255)
                                                                         cGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGa                       <  1:3929033/101‑1 (MQ=255)
                                                                             cgcTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcgcg                   <  2:3913439/101‑1 (MQ=255)
                                                                                      ccccGGCCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCt          <  2:4143109/100‑1 (MQ=255)
                                                                                         cgccccGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGc       <  1:872604/98‑1 (MQ=255)
                                                                                            cccGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAgg    <  1:3639902/101‑1 (MQ=255)
                                                                                             ccGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAggc   <  1:1600302/101‑1 (MQ=255)
                                                                                              cgcgGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAggcg  <  1:3695260/101‑1 (MQ=255)
                                                                                              cgcgGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAggcg  <  1:1232443/101‑1 (MQ=255)
                                                                                                  |                                                                                                
GCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAGGCG  >  NC_014034/158133‑158327

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 38 ≤  ATCG