New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 158231 =64 (0.310)25 (0.120) 4/198 30.4 0.34 intergenic (+34/‑196) parC/tuf1 DNA topoisomerase IV subunit A/translation elongation factor Tu
?NC_014034 1253896 = 32 (0.150)intergenic (+54/+37) RCAP_rcc01181/hisH hypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Rejected: Position hash score below cutoff.

CGCCGGGGGTGCGGGGGCGGAG.CCCC...................................................................................................  <  NC_014034/158256‑158231
..........................CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCG  >  NC_014034/1253896‑1253994
                                                                                                                              
cctCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGC                           >  2:1828390/4‑101 (MQ=33)
 ctCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCA                          >  1:1035545/3‑101 (MQ=33)
 ctCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCA                          >  1:332042/3‑101 (MQ=33)
   CGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGC                        >  2:2552263/1‑101 (MQ=32)
   CGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGC                        >  2:1524725/1‑101 (MQ=32)
    GGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGac                       >  1:2006029/1‑99 (MQ=17)
     GGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCC                      >  2:2151746/1‑101 (MQ=32)
        GTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGT                   >  2:3871944/1‑101 (MQ=32)
         TGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTG                  >  2:4012485/1‑101 (MQ=31)
           CGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCG                >  1:679077/1‑101 (MQ=1)
             GGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCT              >  2:1053110/1‑101 (MQ=1)
              GGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTG             >  1:1688970/1‑101 (MQ=1)
               GGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGC            >  2:913407/1‑101 (MQ=1)
               GGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGC            >  2:2937174/1‑101 (MQ=1)
                gggccccCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCG           <  1:265479/94‑1 (MQ=0)
                  ggagcCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAT         >  2:3724872/6‑101 (MQ=1)
                  ggagcCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAT         >  1:3004295/6‑101 (MQ=1)
                   gagcCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATT        >  2:1990273/5‑101 (MQ=1)
                    cccCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTT       <  1:1951069/98‑1 (MQ=1)
                      cCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCT     >  2:1052656/2‑101 (MQ=1)
                      cCCCCCCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCT     <  2:1133907/100‑1 (MQ=1)
                       CCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTC    >  2:2124978/1‑101 (MQ=1)
                       CCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTC    <  1:580820/101‑1 (MQ=1)
                        CCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCC   <  2:2042140/101‑1 (MQ=1)
                         CCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCG  <  2:3905717/101‑1 (MQ=1)
                                                                                                                              
CGCCGGGGGTGCGGGGGCGGAG.CCCC...................................................................................................  <  NC_014034/158256‑158231
..........................CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCG  >  NC_014034/1253896‑1253994

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG