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

TGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAG  >  NC_014034/158130‑158324
                                                                                                  |                                                                                                
tGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcgcg                                                                                                <  1:1900087/101‑1 (MQ=255)
   gCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCgggg                                                                                             <  1:3247032/101‑1 (MQ=255)
   gCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCgggg                                                                                             <  1:2230775/101‑1 (MQ=255)
    cAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCcgggc                                                                                            >  2:3033388/1‑100 (MQ=255)
     aaGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcgccgcggg                                                                                           >  1:2258199/1‑95 (MQ=255)
     aaGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGc                                                                                           <  2:1362372/101‑1 (MQ=255)
     aaGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCCCgggggg                                                                                           >  1:2304141/1‑100 (MQ=255)
      aGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCcggccgg                                                                                          >  2:672012/1‑97 (MQ=255)
      aGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCt                                                                                          <  2:833329/101‑1 (MQ=255)
        cGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTcc                                                                                        <  1:3154870/101‑1 (MQ=255)
          ggCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGc                                                                                      <  1:4309948/101‑1 (MQ=255)
            cGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccc                                                                                    <  2:4362751/101‑1 (MQ=255)
             gAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGcccc                                                                                   <  1:2450127/101‑1 (MQ=255)
              aGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccccc                                                                                  <  2:4664806/101‑1 (MQ=255)
              aGCGCCGGGCGGAGGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccccc                                                                                  <  2:3274741/101‑1 (MQ=255)
               gcgcCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCg                                                                                 <  2:2485604/101‑1 (MQ=255)
                  ccGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAc                                                                              <  2:4197749/101‑1 (MQ=255)
                     ggCGGATGGCGCCGCGGGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCCCGGGGGCTCCGCCCCCGCAcccc                                                                           <  1:2851723/101‑1 (MQ=255)
                      gCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAccccc                                                                          <  2:3247259/101‑1 (MQ=255)
                      gCGGATGGCGCCGCGCGGTTTCCCGCGCAAGACCCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAccccc                                                                          <  1:120491/101‑1 (MQ=255)
                        ggATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCgg                                                                        <  1:1700025/101‑1 (MQ=255)
                          aTGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCg                                                                      <  1:893014/101‑1 (MQ=255)
                            ggCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTa                                                                    <  2:660108/101‑1 (MQ=255)
                               gccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAttt                                                                 <  1:2949860/101‑1 (MQ=255)
                               gccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAttt                                                                 <  1:1281937/101‑1 (MQ=255)
                                     cgGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCaa                                                           <  1:1128108/101‑1 (MQ=255)
                                     cgGTTGCCCGCGCGACAACCGGGTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCaa                                                           <  2:799993/101‑1 (MQ=255)
                                      ggTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAg                                                          <  1:293362/101‑1 (MQ=255)
                                          tCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaaa                                                      <  1:3903266/101‑1 (MQ=255)
                                          tCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaaa                                                      <  1:1619007/101‑1 (MQ=255)
                                           cccGCGCGACAACCGGGTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaaaa                                                     <  2:1169000/101‑1 (MQ=255)
                                                 cgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAgg                                               <  2:3027095/101‑1 (MQ=255)
                                                    cAACCGGTTCCCCTGACGGGGTNTCGNGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGc                                            <  1:4163224/101‑1 (MQ=255)
                                                     aaCCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGca                                           <  2:4657794/101‑2 (MQ=255)
                                                      aCCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGccc                                          <  1:550975/101‑1 (MQ=255)
                                                      aCCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGccc                                          <  1:1840385/101‑1 (MQ=255)
                                                          gTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCg                                      <  2:2094862/101‑1 (MQ=255)
                                                           ttCCCCTGACGGGGTTTCGGGCTGTGCCGCGACCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCgg                                     <  1:4146215/101‑1 (MQ=255)
                                                             ccccTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGAc                                   <  2:983427/101‑1 (MQ=255)
                                                                cTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgtg                                <  1:3306497/101‑1 (MQ=255)
                                                                cTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgtg                                <  2:4212519/101‑1 (MQ=255)
                                                                  gACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTc                              <  2:4564513/101‑1 (MQ=255)
                                                                  gACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTc                              <  1:590811/101‑1 (MQ=255)
                                                                   gcGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCg                             <  1:2048120/100‑1 (MQ=255)
                                                                   gcGGGGGTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCg                             <  2:3687782/100‑1 (MQ=255)
                                                                      gggTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGttt                          <  1:3607299/101‑1 (MQ=255)
                                                                        gTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCg                        <  2:2021784/101‑1 (MQ=255)
                                                                            cGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGa                    <  2:1597340/101‑1 (MQ=255)
                                                                             gggCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAc                   <  2:4400186/101‑1 (MQ=255)
                                                                                    gCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCt            <  2:1146514/101‑1 (MQ=255)
                                                                                       gcgcCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGcc         <  1:2057366/101‑1 (MQ=255)
                                                                                       gcgcCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGcc         <  2:181649/101‑1 (MQ=255)
                                                                                         gcCCGGCCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCt       <  1:630731/101‑1 (MQ=255)
                                                                                         gcCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCt       <  1:1983597/101‑1 (MQ=255)
                                                                                          cccGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTg      <  1:2666837/101‑1 (MQ=255)
                                                                                           ccGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGATGTGTCGTTTCGGGGACGCGCGCTGCCCTgg     <  1:4148353/101‑1 (MQ=255)
                                                                                              gcccGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAg  <  2:103294/100‑1 (MQ=255)
                                                                                                  |                                                                                                
TGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGGCAG  >  NC_014034/158130‑158324

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