New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 1699965 =30 (0.170)21 (0.120) 10/168 22.8 0.26 intergenic (+38/+5) RCAP_rcc01568/purF RES domain family protein/amidophosphoribosyltransferase
?NC_014034 = 2590221 90 (0.420)intergenic (‑120/‑136) glnE/RCAP_rcc02415 glutamate‑ammonia‑ligase adenylyltransferase/hypothetical protein
Rejected: Position hash score below cutoff.

GGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTT  >  NC_014034/1699866‑1700065
                                                                                                     |                                                                                                    
ggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTAGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATAcc                                                                                                     >  2:1102100/1‑101 (MQ=255)
  agggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAAGAcacc                                                                                                   >  1:4171026/1‑101 (MQ=255)
    ggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCt                                                                                                 <  2:3230491/101‑1 (MQ=255)
     gagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTc                                                                                                <  1:1576207/101‑1 (MQ=255)
     gagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTc                                                                                                <  2:2194860/101‑1 (MQ=255)
        ccgccgATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGc                                                                                             <  2:1440701/101‑1 (MQ=255)
                                       gCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATc                                                              <  2:4170764/101‑1 (MQ=255)
                                                              ccccgccgggtCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCg                                        >  2:3043121/11‑101 (MQ=25)
                                                                 cgccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGccc                                     >  2:2800946/8‑101 (MQ=37)
                                                                        tcgtgtgggC.TGC.GCCGAGTGTAAATGCCCCCTCAGCCGTTCATCTTGACGCCCTGGCCGATGTGATCGGAGGTCTTCACCGGATAATCGCCCGAGAAACa                             <  1:1507165/94‑1 (MQ=255)
                                                                        tCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGGTCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAAc                              <  2:56039/101‑1 (MQ=38)
                                                                             ctcggcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGc                          <  2:15004/91‑1 (MQ=255)
                                                                             ctcggcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGc                          <  1:2283838/91‑1 (MQ=255)
                                                                              tcggcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCg                         >  2:1695697/10‑101 (MQ=255)
                                                                               cggcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCNTCACCGGATAATCGCCCGAGAACCAGGCGt                        >  1:1697574/9‑101 (MQ=255)
                                                                                ggcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTc                       >  1:582912/8‑101 (MQ=255)
                                                                                ggcgtcagcCGAGGCGAAATAACCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTc                       >  1:2406995/8‑101 (MQ=255)
                                                                                 gcgtcagcCGAGGCGAAATACACCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCg                      >  1:3909811/7‑101 (MQ=255)
                                                                                  cgtcagcCGAGGCGAAATCCCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGc                     <  2:4229939/96‑1 (MQ=255)
                                                                                   gtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCa                    >  2:2355603/5‑101 (MQ=255)
                                                                                   gtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCa                    >  2:2532474/5‑101 (MQ=255)
                                                                                    tcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCaa                   >  1:307762/4‑101 (MQ=255)
                                                                                     cagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAat                  >  2:132512/3‑101 (MQ=255)
                                                                                     cagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAat                  >  1:3035313/3‑101 (MQ=255)
                                                                                       gcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAatat                >  2:2657110/1‑101 (MQ=255)
                                                                                       gcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAatat                >  2:849994/1‑101 (MQ=255)
                                                                                       gcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCACTTGCCGATCTGATCGGAGGGATTCACCGGATAATAGCCCGAGAAACAGGCGTCGCAatat                >  2:1907400/1‑101 (MQ=255)
                                                                                        ccGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATAtt               >  1:89252/1‑101 (MQ=255)
                                                                                        ccGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATAtt               >  1:3592610/1‑101 (MQ=255)
                                                                                        ccGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATAtt               >  1:218289/1‑101 (MQ=255)
                                                                                         cGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTg              <  2:1281045/101‑1 (MQ=255)
                                                                                         cGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTg              <  1:1234641/101‑1 (MQ=255)
                                                                                           aGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCg            >  1:2898819/1‑101 (MQ=255)
                                                                                            ggCGAAATACCCCGCCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTGACCGGATAAGCGCCCGAGAAAGAGGCGTCGCAAGATTggcg           >  2:4217425/1‑98 (MQ=255)
                                                                                             gCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGa          <  2:2911840/101‑1 (MQ=255)
                                                                                              cGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGAc         <  2:2031949/101‑1 (MQ=255)
                                                                                              cGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATAGCCCGAGAAACAGGCGTCGCAATATTACGGAc         >  2:2599244/1‑101 (MQ=255)
                                                                                               gAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACa        >  2:906090/1‑101 (MQ=255)
                                                                                               gAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACa        >  1:3849869/1‑101 (MQ=255)
                                                                                                aaaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAg       <  2:3352092/101‑1 (MQ=255)
                                                                                                aaaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAg       >  2:3941422/1‑101 (MQ=255)
                                                                                                aaaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAg       >  1:4095947/1‑101 (MQ=255)
                                                                                                 aaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAgg      <  2:3081929/101‑1 (MQ=255)
                                                                                                 aaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAgg      >  2:864458/1‑101 (MQ=255)
                                                                                                  aTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGc     >  1:1148214/1‑101 (MQ=255)
                                                                                                   tACCCCCTCAGCCCTTCATCTTGAAGCCCTCGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCg    >  2:2138061/1‑101 (MQ=255)
                                                                                                    aCCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGt   >  1:3335842/1‑101 (MQ=255)
                                                                                                     cccccTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGtt  <  1:2661107/101‑1 (MQ=255)
                                                                                                     |                                                                                                    
GGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTT  >  NC_014034/1699866‑1700065

Base quality scores:  ATCG  < 3 ≤  ATCG  < 21 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG