New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 158228 =66 (0.300)53 (0.240) 20/194 21.9 0.52 intergenic (+31/‑199) parC/tuf1 DNA topoisomerase IV subunit A/translation elongation factor Tu
?NC_014034 = 2330172 30 (0.130)intergenic (+147/+165) RCAP_rcc02161/mrcA hypothetical protein/penicillin‑binding protein 1A
Rejected: Position hash score below cutoff.

GCTGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAGGCG  >  NC_014034/158128‑158327
                                                                                                    |                                                                                                   
gCTGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcg                                                                                                     <  2:2279687/101‑1 (MQ=255)
gCTGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcg                                                                                                     >  2:625203/1‑101 (MQ=255)
gCTGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcg                                                                                                     <  2:51739/101‑1 (MQ=255)
 cTGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcgc                                                                                                    >  2:1749080/1‑101 (MQ=255)
  tGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcgcg                                                                                                   <  2:689069/101‑1 (MQ=255)
  tGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcgcg                                                                                                   <  1:4111752/101‑1 (MQ=255)
  tGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcgcc                                                                                                   >  2:4377055/1‑100 (MQ=255)
     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:353666/101‑1 (MQ=255)
               gAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGcccc                                                                                      <  2:2184386/101‑1 (MQ=255)
                agggcCGGGCGGATGGCGACGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccccc                                                                                     <  2:2653834/98‑1 (MQ=255)
                aGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccccc                                                                                     <  1:2608649/101‑1 (MQ=255)
                aGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccccc                                                                                     <  1:720864/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:663020/101‑1 (MQ=255)
                   gcCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCa                                                                                  <  1:4313415/101‑1 (MQ=255)
                   gcCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCa                                                                                  <  2:1139184/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                                                     <  1:4377326/101‑1 (MQ=255)
                                                gcgcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCa                                                     <  2:172691/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:918412/101‑1 (MQ=255)
                                                                           tttCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCgg                          <  1:1790807/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)
                                                                                                    |                                                                                                   
GCTGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAGGCG  >  NC_014034/158128‑158327

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