New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 39246232 (0.150)24 (0.110) 4/198 30.4 0.43 intergenic (‑110/+33) RCAP_rcc00363/trmD hypothetical protein/tRNA (guanine‑N1)‑methyltransferase
?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.

CTCCTTGTCACCGCGCAAGACTATGCGGCCCAAACCGGAGGGACAATGGCAGGCCATTTCATCTTGCCTCAAATATCCCGGGGGTCCGGGGGCGGAG.CCCC...................................................................................................  >  NC_014034/392362‑392462
.....................................................................................................CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCG  >  NC_014034/1253896‑1253994
                                                                                                                                                                                                         
                                                                            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/4‑101 (MQ=32)
                                                                                    tgCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTG                  >  2:4012485/3‑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)
                                                                                                                                                                                                         
CTCCTTGTCACCGCGCAAGACTATGCGGCCCAAACCGGAGGGACAATGGCAGGCCATTTCATCTTGCCTCAAATATCCCGGGGGTCCGGGGGCGGAG.CCCC...................................................................................................  >  NC_014034/392362‑392462
.....................................................................................................CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCG  >  NC_014034/1253896‑1253994

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