New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 158228 =65 (0.320)37 (0.180) 17/194 20.9 0.43 intergenic (+31/‑199) parC/tuf1 DNA topoisomerase IV subunit A/translation elongation factor Tu
?NC_014034 = 2330172 35 (0.170)intergenic (+147/+165) RCAP_rcc02161/mrcA hypothetical protein/penicillin‑binding protein 1A
Rejected: Position hash score below cutoff.

CTGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGG  >  NC_014034/158129‑158321
                                                                                                   |                                                                                             
cTGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcgc                                                                                              >  2:2031088/1‑101 (MQ=255)
 tGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcgcg                                                                                             >  1:3511439/1‑101 (MQ=255)
 tGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcgcg                                                                                             <  2:625656/101‑1 (MQ=255)
  gggCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCgg                                                                                            <  1:2754677/101‑1 (MQ=255)
   gccAAGCGGGCGAGCGACGGTCGGGTGGCGGCGCGCGGTTTGCCGCGCTGCAACCGGTTCACCTGACGGGGTTTCGGGCGCTGCCGCACCCGACACGCggg                                                                                           <  2:2717809/99‑1 (MQ=255)
     cAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCCGCCGCGCCCGACCCGCggggg                                                                                         <  1:2400178/101‑1 (MQ=255)
       aGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCt                                                                                       <  1:1443798/101‑1 (MQ=255)
        gCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTc                                                                                      <  1:3720198/101‑1 (MQ=255)
        gCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTc                                                                                      <  1:4138733/101‑1 (MQ=255)
         cGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAGCCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTcc                                                                                     <  1:2718556/101‑1 (MQ=255)
         cGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTcc                                                                                     <  2:531793/101‑1 (MQ=255)
          gggCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCcgccgacgggccc                                                                                    >  2:219437/1‑91 (MQ=255)
           ggCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGc                                                                                   <  1:875704/101‑1 (MQ=255)
            gCGCGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCACGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGTCCTTCCCGCGTTCGCTCCGcc                                                                                  <  2:3511186/101‑1 (MQ=255)
             cGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGccc                                                                                 <  1:1590976/101‑1 (MQ=255)
              gAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGcccc                                                                                <  2:1454663/101‑1 (MQ=255)
                   ccGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAc                                                                           <  1:3735030/101‑1 (MQ=255)
                     gggCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAccc                                                                         <  2:516418/101‑1 (MQ=255)
                        cGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCAcccccc                                                                      <  1:2873971/101‑2 (MQ=255)
                         tgATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCAGCACCCCCgg                                                                     <  1:3970716/100‑1 (MQ=255)
                          gATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGc                                                                    <  1:437954/101‑1 (MQ=255)
                             ggCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTa                                                                 <  2:2261613/101‑1 (MQ=255)
                             ggCGCCGCGCGGTTCCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTa                                                                 <  2:1242703/101‑1 (MQ=255)
                             ggCCCCGCGCGGGTTCCCGCGCCACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCCCCCGACCCGCGGGGGCTCCCGCCCCGGACCCCCGGCgga                                                                 <  2:4075815/101‑3 (MQ=255)
                              gcgcCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAt                                                                <  1:2292356/101‑1 (MQ=255)
                              gcgcCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAt                                                                <  1:133170/101‑1 (MQ=255)
                               cgccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtt                                                               <  1:1132999/101‑1 (MQ=255)
                                gccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAttt                                                              <  2:2516726/101‑1 (MQ=255)
                                 ccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtttt                                                             <  2:228368/101‑1 (MQ=255)
                                 ccgcGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTAtttt                                                             <  2:3562982/101‑1 (MQ=255)
                                  cgcgcgGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTg                                                            <  1:2031226/101‑1 (MQ=255)
                                  cgcgcgGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTg                                                            <  1:3009398/101‑1 (MQ=255)
                                       ggTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAg                                                       <  2:702820/101‑1 (MQ=255)
                                           tCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaaa                                                   <  1:2208524/101‑1 (MQ=255)
                                            cccGCGCGACAACCGTTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaaaa                                                  <  2:1332295/101‑1 (MQ=255)
                                            cccGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaaaa                                                  <  1:2908053/101‑1 (MQ=255)
                                            cccGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGaaaa                                                  <  1:2619926/101‑1 (MQ=255)
                                                cgcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCaa                                              <  2:2034292/101‑1 (MQ=255)
                                                cgcgACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCaa                                              <  1:2708767/101‑1 (MQ=255)
                                                   gACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGt                                           <  1:3477323/101‑1 (MQ=255)
                                                     cacccGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGCCCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGc                                         <  2:3753061/98‑1 (MQ=255)
                                                      aaCCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGcc                                        <  2:3767569/101‑1 (MQ=255)
                                                       aggggTTCCCCTGACGGGGGTGCGGGCGCTGCCGCGCCCGCCCCGCGGGGGCTCCGCCCCCGCACGCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGccc                                       <  2:1569950/98‑1 (MQ=37)
                                                        ccGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGcccc                                      <  2:3233003/101‑1 (MQ=255)
                                                          gtttCCCCTGACGGGGTTTCGGGCGGTGCCGCCCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGcccccc                                    <  1:3964507/99‑1 (MQ=255)
                                                          ggTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGcccccc                                    <  2:3986989/101‑1 (MQ=255)
                                                           gTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCg                                   <  2:691190/101‑1 (MQ=255)
                                                            ttCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCgg                                  <  2:2535233/101‑1 (MQ=255)
                                                            ttCCCCGGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCgg                                  <  1:755045/101‑1 (MQ=255)
                                                              ccccTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGAc                                <  2:3954211/101‑1 (MQ=255)
                                                                ccTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgt                              <  1:2412699/101‑1 (MQ=255)
                                                                 cTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgtg                             <  2:3994755/101‑1 (MQ=255)
                                                                 cTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgtg                             <  2:918737/101‑1 (MQ=255)
                                                                 cTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgtg                             <  2:812729/101‑1 (MQ=255)
                                                                  tGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACgtgt                            <  2:2209994/101‑1 (MQ=255)
                                                                     cGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGt                         <  2:1883179/101‑1 (MQ=255)
                                                                            tCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCTCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCgggg                  <  1:3571874/101‑1 (MQ=255)
                                                                            tCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCgggg                  <  1:3666088/101‑1 (MQ=255)
                                                                               ggCGCTGCCGCGCCCGGCCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcg               <  2:3841580/101‑1 (MQ=255)
                                                                               ggCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcg               <  2:3231975/101‑1 (MQ=255)
                                                                                gcgcGGGCGCGCCCGGCCCGCGGGGGCCCCGCCCCCGGACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcgc              <  2:1837796/101‑1 (MQ=38)
                                                                                  gcTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcgcgc            <  2:3410770/101‑1 (MQ=255)
                                                                                  gcTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcgcgc            <  1:823773/101‑1 (MQ=255)
                                                                                   cTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGAcgcgcg           <  2:3027202/101‑1 (MQ=255)
                                                                                      ccGCGGCTGCCCCGCGGGGGCTCCGCCCCCGCTCCCCCGGCGTATTTTGCCAAGGAGAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTg        <  2:2512590/101‑1 (MQ=38)
                                                                                      ccGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTg        <  2:2177597/101‑1 (MQ=255)
                                                                                       cgcgCCCGACCCGCGGGGGCTCCGCCCCCGCGCCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGc       <  2:2117500/101‑1 (MQ=255)
                                                                                            ccGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTgg  <  2:2242647/101‑1 (MQ=255)
                                                                                            ccGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTgg  <  2:535078/101‑1 (MQ=255)
                                                                                                   |                                                                                             
CTGGGCAAGCGGGCGAGCGCCGGGCGGATGGCGCCGCGCGGTTTCCCGCGCGACAACCGGTTCCCCTGACGGGGTTTCGGGCGCTGCCGCGCCCGACCCGCGGGGGCTCCGCCCCCGCACCCCCGGCGTATTTTGCCAAGGAAAAGCAAGGTGCCCCCCGGACGTGTCGTTTCGGGGACGCGCGCTGCCCTGG  >  NC_014034/158129‑158321

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 38 ≤  ATCG