New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 1699965 =52 (0.250)26 (0.130) 15/168 17.9 0.21 intergenic (+38/+5) RCAP_rcc01568/purF RES domain family protein/amidophosphoribosyltransferase
?NC_014034 = 2590221 141 (0.580)intergenic (‑120/‑136) glnE/RCAP_rcc02415 glutamate‑ammonia‑ligase adenylyltransferase/hypothetical protein
Rejected: Position hash score below cutoff.

GGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTT  >  NC_014034/1699870‑1700065
                                                                                                 |                                                                                                    
ggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCt                                                                                                 >  2:374325/1‑101 (MQ=255)
         gATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAAAACCCCCTCAGCCCTTc                                                                                        >  2:33599/1‑101 (MQ=255)
            ctcGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGCGGCGAAATAACCCCCCAGCCCTTCATc                                                                                     >  1:4406060/1‑101 (MQ=255)
                       tGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCtt                                                                          <  1:192407/101‑1 (MQ=255)
                            ggCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGa                                                                     <  1:374378/101‑1 (MQ=255)
                                 tctccAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGa                                                                <  2:319271/98‑1 (MQ=255)
                                    ccAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGaaga                                                             >  1:319324/1‑98 (MQ=255)
                                           gCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGc                                                       >  1:846498/1‑101 (MQ=38)
                                                       tGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAACCGCCCCCCCCGCCCTTCCTCTTGCCGCCCCTGCCGCTCTGCTCGGAGGGCTCCACAGGATaac                                          >  1:4103892/1‑100 (MQ=19)
                                                          gcggcccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCg                                        <  2:1831146/91‑1 (MQ=25)
                                                          ccccgccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCg                                        >  2:2831891/11‑101 (MQ=25)
                                                             ccccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGccc                                     <  2:1131170/94‑1 (MQ=37)
                                                              cccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCg                                    <  2:3552265/95‑1 (MQ=37)
                                                               ccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCga                                   <  1:1441376/96‑1 (MQ=37)
                                                                 gggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgcga                                 >  2:4291865/4‑101 (MQ=37)
                                                                 gggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgaga                                 >  1:198195/4‑101 (MQ=38)
                                                                    tCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAAc                              >  2:2734351/1‑101 (MQ=38)
                                                                    tCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAAc                              <  2:3480562/101‑1 (MQ=38)
                                                                    tCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAAc                              >  1:2897567/1‑101 (MQ=38)
                                                                     cGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACa                             >  2:681178/1‑101 (MQ=38)
                                                                         tCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGt                        <  1:2237897/101‑1 (MQ=255)
                                                                         ctcggcgtcagcCGAGGCGAAATACCCCCTCATCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGc                          >  1:3109504/11‑101 (MQ=39)
                                                                         ctcggcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGc                          >  2:3590649/11‑101 (MQ=255)
                                                                           cggcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGt                        <  1:4787737/93‑1 (MQ=255)
                                                                            gC.TGC.GCCGAGGCGAAATCCCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGc                     <  2:114850/101‑1 (MQ=255)
                                                                             gggtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCg                      <  1:4929857/95‑1 (MQ=255)
                                                                             gcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCg                      >  2:4983471/7‑101 (MQ=255)
                                                                             gcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCg                      >  2:4313377/7‑101 (MQ=255)
                                                                               gtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCa                    >  2:1980214/5‑101 (MQ=255)
                                                                               gtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCa                    >  2:3266497/5‑101 (MQ=255)
                                                                               gtcagcCGAGGCGAAATACCACCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCa                    >  2:4127522/5‑101 (MQ=255)
                                                                                tcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCaa                   >  1:639210/4‑101 (MQ=255)
                                                                                tcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCaa                   >  2:4072515/4‑101 (MQ=255)
                                                                                tcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCaa                   >  2:4769688/4‑101 (MQ=255)
                                                                                 cagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAat                  >  1:4598535/3‑101 (MQ=255)
                                                                                 cagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAat                  <  1:234444/99‑1 (MQ=255)
                                                                                 cagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAat                  >  2:4886977/3‑101 (MQ=255)
                                                                                 cagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGAGAATCGCCCGAGAAACAGGCGTCGCAat                  >  1:4460779/3‑101 (MQ=255)
                                                                                  agcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAata                 >  1:1517310/2‑101 (MQ=255)
                                                                                   gcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAatat                >  2:1176786/1‑101 (MQ=255)
                                                                                   gcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAatat                >  1:2255890/1‑101 (MQ=255)
                                                                                    ccGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATAtt               <  2:2431847/101‑1 (MQ=255)
                                                                                    ccGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATAtt               >  2:2961500/1‑101 (MQ=255)
                                                                                     cGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCGCCGCATAATCGCCCGAGAAACAGGCGTCGCAATATTg              >  1:4323530/1‑101 (MQ=255)
                                                                                      gAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGc             >  1:2159964/1‑101 (MQ=255)
                                                                                       aGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCg            >  2:2370086/1‑101 (MQ=255)
                                                                                          cGAAATCCCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGAc         <  2:4743563/101‑1 (MQ=255)
                                                                                          cGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGAc         <  2:639121/101‑1 (MQ=255)
                                                                                          cGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGAc         <  1:1875427/101‑1 (MQ=255)
                                                                                          cGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGAc         >  1:3288303/1‑101 (MQ=255)
                                                                                           gAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGAcg        >  2:4952029/1‑100 (MQ=255)
                                                                                           gAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACa        >  1:906668/1‑101 (MQ=255)
                                                                                           gAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACa        >  1:277675/1‑101 (MQ=255)
                                                                                           gAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACa        <  1:2952907/101‑1 (MQ=255)
                                                                                            aaaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGTCCGAGAAACAGGCGTCGCAATATTGCGGACAg       >  2:3708484/1‑101 (MQ=255)
                                                                                            aaaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAg       >  2:2375993/1‑101 (MQ=255)
                                                                                            aaaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAg       >  1:4180143/1‑101 (MQ=255)
                                                                                            aaaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAg       >  2:1762258/1‑101 (MQ=255)
                                                                                            aaaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGATTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAg       >  1:1690278/1‑101 (MQ=255)
                                                                                             aaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAgg      >  2:1032856/1‑101 (MQ=255)
                                                                                             aaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAgg      >  1:4421767/1‑101 (MQ=255)
                                                                                             aaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAgg      >  2:1244204/1‑101 (MQ=255)
                                                                                             aaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAgg      <  1:1564368/101‑1 (MQ=255)
                                                                                             aaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAgg      <  2:4460460/101‑1 (MQ=255)
                                                                                             aaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAgg      >  1:3656180/1‑101 (MQ=255)
                                                                                             aaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAgg      >  2:1473017/1‑101 (MQ=255)
                                                                                              aTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGc     >  2:37822/1‑101 (MQ=255)
                                                                                               tACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGGCAGGCg    >  2:2438752/1‑101 (MQ=255)
                                                                                               tACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCg    >  2:3382193/1‑101 (MQ=255)
                                                                                               tACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCg    <  1:3623490/101‑1 (MQ=255)
                                                                                                 cccccTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGtt  <  1:3899099/101‑1 (MQ=255)
                                                                                                 |                                                                                                    
GGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTT  >  NC_014034/1699870‑1700065

Base quality scores:  ATCG  < 3 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG