New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 39246233 (0.140)26 (0.110) 2/198 38.6 0.43 intergenic (‑110/+33) RCAP_rcc00363/trmD hypothetical protein/tRNA (guanine‑N1)‑methyltransferase
?NC_014034 1253896 = 37 (0.160)intergenic (+54/+37) RCAP_rcc01181/hisH hypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Rejected: Position hash score below cutoff.

CTCCTTGTCACCGCGCAAGACTATGCGGCCCAAACCGGAGGGACAATGGCAGGCCATTTCATCTTGCCTCAAATATCCCGGGGGTCCGGGGGCGGAG.CCCC..................................................................................................  >  NC_014034/392362‑392462
.....................................................................................................CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCC  >  NC_014034/1253896‑1253993
                                                                                                                                                                                                        
                                                                            ctCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCA                         >  1:2045009/3‑101 (MQ=33)
                                                                            ctCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCA                         >  1:586145/3‑101 (MQ=32)
                                                                              CGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGC                       >  2:2014047/1‑101 (MQ=32)
                                                                              CGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGC                       >  1:2731371/1‑101 (MQ=32)
                                                                              CGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGC                       >  1:4502706/1‑101 (MQ=32)
                                                                               GGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCC                      >  2:2977579/1‑101 (MQ=32)
                                                                               GGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCC                      >  1:675259/1‑101 (MQ=32)
                                                                                  GGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGG                   >  1:4245193/1‑101 (MQ=32)
                                                                                       GGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGC              >  1:1073280/1‑101 (MQ=1)
                                                                                        gggggggcgcCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCT             <  2:4091353/91‑1 (MQ=11)
                                                                                        GGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCT             >  2:3606556/1‑101 (MQ=1)
                                                                                         ggccgccccCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTG            <  2:183050/92‑1 (MQ=0)
                                                                                          ggcgccccCCCGGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGC           <  2:2234407/93‑1 (MQ=0)
                                                                                             ggagcCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAT        >  2:3873057/6‑101 (MQ=1)
                                                                                             ggagcCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAT        >  2:2257756/6‑101 (MQ=1)
                                                                                              gcgcCCCGGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATT       <  1:2919127/97‑1 (MQ=1)
                                                                                               agcCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTT      >  1:523554/4‑101 (MQ=1)
                                                                                               agcCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTT      >  2:1102081/4‑101 (MQ=1)
                                                                                                gcCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTC     >  1:778416/3‑101 (MQ=1)
                                                                                                 cCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCT    >  2:4127384/2‑101 (MQ=1)
                                                                                                 cCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCT    >  2:4169726/2‑101 (MQ=1)
                                                                                                 cCCCCCCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCT    <  1:4069719/100‑1 (MQ=1)
                                                                                                  CCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTC   >  1:3548841/1‑101 (MQ=1)
                                                                                                   CCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCC  <  1:693072/101‑1 (MQ=1)
                                                                                                   CCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCC  <  1:3177504/101‑1 (MQ=1)
                                                                                                   CCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCC  <  1:2262720/101‑1 (MQ=1)
                                                                                                                                                                                                        
CTCCTTGTCACCGCGCAAGACTATGCGGCCCAAACCGGAGGGACAATGGCAGGCCATTTCATCTTGCCTCAAATATCCCGGGGGTCCGGGGGCGGAG.CCCC..................................................................................................  >  NC_014034/392362‑392462
.....................................................................................................CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCC  >  NC_014034/1253896‑1253993

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