New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 1699978 =63 (0.440)5 (0.040) 2/120 28.8 0.08 coding (1459/1467 nt) purF amidophosphoribosyltransferase
?NC_014034 = 2995646 46 (0.190)intergenic (+186/+34) RCAP_rcc02802/RCAP_rcc02803 membrane dipeptidase/hypothetical protein
Rejected: Position hash score below cutoff.

TCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTCGCGCCCGGCC  >  NC_014034/1699883‑1700078
                                                                                                 |                                                                                                    
tcGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCt                                                                                                 <  1:445212/101‑1 (MQ=255)
  ggCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTg                                                                                               <  1:3914188/101‑1 (MQ=255)
       aaCTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCCAGGCGAAATACCCCCTCCGCCCTTCATCTTGAAgac                                                                                          >  2:3405546/1‑99 (MQ=255)
            aCCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGc                                                                                     <  1:3782175/101‑1 (MQ=255)
                                       ggTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGAt                                                          <  2:3717718/101‑1 (MQ=255)
                                             cacgcccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCg                                                     <  2:540674/91‑1 (MQ=25)
                                              cccgccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCGCCGGATAAGCGc                                                    >  2:145965/10‑101 (MQ=25)
                                                gcccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGccc                                                  <  1:361603/94‑1 (MQ=37)
                                                aaGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCg                                                 >  1:2482818/1‑101 (MQ=255)
                                                   cgggtCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgag                                               >  1:2832272/5‑101 (MQ=37)
                                                   cgggtCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgag                                               >  2:4561554/5‑101 (MQ=37)
                                                      gtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGaaa                                            <  2:4432783/100‑1 (MQ=38)
                                                      gtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGaaa                                            >  1:2412615/2‑101 (MQ=38)
                                                       tCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACa                                          <  1:4507670/101‑1 (MQ=255)
                                                        cGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACa                                          >  2:1153614/1‑101 (MQ=38)
                                                           cTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAggag                                      >  2:1892201/1‑99 (MQ=255)
                                                            ctcggcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGc                                       >  1:1782454/11‑101 (MQ=255)
                                                             cGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTc                                    >  2:795182/1‑101 (MQ=255)
                                                              cggcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGt                                     >  2:3996622/9‑101 (MQ=255)
                                                               ggcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTc                                    >  2:2100789/8‑101 (MQ=255)
                                                               ggcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTc                                    >  1:761438/8‑101 (MQ=255)
                                                               ggcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTc                                    >  1:4472076/8‑101 (MQ=255)
                                                                gcgtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCg                                   >  1:767727/7‑101 (MQ=255)
                                                                 cgtcagcCGAGGCGAAATACCCCCTCCGCCCTTCATCTTGAAGCCCTTGCCGCTCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGc                                  >  1:1859364/6‑101 (MQ=255)
                                                                  gtcagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCa                                 <  2:1258817/97‑1 (MQ=255)
                                                                   tcagcCGAGGCGGAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCaa                                >  1:2523883/4‑101 (MQ=255)
                                                                    cagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAat                               <  1:1187055/99‑1 (MQ=255)
                                                                    cagcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCCGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAat                               >  1:3875477/3‑101 (MQ=255)
                                                                     agcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAata                              >  1:3474956/2‑101 (MQ=255)
                                                                     agcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAata                              <  1:2228524/100‑1 (MQ=255)
                                                                     agcCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGACCGAGAAACAGGCGTCGCAata                              >  2:878774/2‑101 (MQ=255)
                                                                       ccGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATAtt                            >  1:1358736/1‑101 (MQ=255)
                                                                        cGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTg                           >  1:1101084/1‑101 (MQ=255)
                                                                         taGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGc                          >  1:2888099/2‑101 (MQ=255)
                                                                         gAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGc                          <  1:2260709/101‑1 (MQ=255)
                                                                          aGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCg                         >  2:2106531/1‑101 (MQ=255)
                                                                          aGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCg                         >  1:637243/1‑101 (MQ=255)
                                                                            gCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGa                       <  1:2163703/101‑1 (MQ=255)
                                                                             cGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGAc                      >  1:2636678/1‑101 (MQ=255)
                                                                              gAAATACTCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACa                     >  1:876200/1‑101 (MQ=255)
                                                                               aaaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAg                    >  2:3717013/1‑101 (MQ=255)
                                                                               aaaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAg                    >  2:2758187/1‑101 (MQ=255)
                                                                               aaaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAg                    >  2:3005957/1‑101 (MQ=255)
                                                                               aaaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAg                    >  1:4093787/1‑101 (MQ=255)
                                                                               aaaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAg                    >  1:3812141/1‑101 (MQ=255)
                                                                                aaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAgg                   >  1:3994554/1‑101 (MQ=255)
                                                                                aaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAgg                   >  1:2335026/1‑101 (MQ=255)
                                                                                aaTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAgg                   >  1:2103416/1‑101 (MQ=255)
                                                                                 aTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGc                  <  2:1058464/101‑1 (MQ=255)
                                                                                 aTACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGc                  >  1:2042680/1‑101 (MQ=255)
                                                                                   tcccccTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGt                <  1:2359928/100‑1 (MQ=255)
                                                                                   aCCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGt                >  2:1877326/1‑101 (MQ=255)
                                                                                   aCCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGt                >  2:2220344/1‑101 (MQ=255)
                                                                                    cccccTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGtt               <  1:3778086/101‑1 (MQ=255)
                                                                                    cccccTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGtt               >  2:2105551/1‑101 (MQ=255)
                                                                                    cccccTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGtt               >  2:3276215/1‑101 (MQ=255)
                                                                                    cccccTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGta               >  2:2359766/1‑100 (MQ=255)
                                                                                     ccccTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTg              <  2:3811856/101‑1 (MQ=255)
                                                                                        cTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTcg           <  2:951170/101‑1 (MQ=255)
                                                                                        cTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTcg           >  1:2687949/1‑101 (MQ=255)
                                                                                          cAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTcgcg         >  2:1752518/1‑101 (MQ=255)
                                                                                          cAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTcgcg         >  2:1278956/1‑101 (MQ=255)
                                                                                          cAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTcgcg         >  2:1249735/1‑101 (MQ=255)
                                                                                          cAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTcgcg         >  1:3206366/1‑101 (MQ=255)
                                                                                          cAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTcgcg         >  2:4593170/1‑101 (MQ=255)
                                                                                          cAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTcgcg         >  1:106335/1‑101 (MQ=255)
                                                                                            gCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTCGCGcc       <  1:2567600/101‑1 (MQ=255)
                                                                                             cccTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTCGCGccc      >  2:3341703/1‑101 (MQ=255)
                                                                                             cccTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTCGCGccc      >  2:2040495/1‑101 (MQ=255)
                                                                                             cccTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTCGCGccc      <  1:2100929/101‑1 (MQ=255)
                                                                                              ccTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTCGCGCCCg     >  1:2301067/1‑101 (MQ=255)
                                                                                                ttCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTCGCGCCCGGc   >  1:4339979/1‑101 (MQ=255)
                                                                                                 tCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTCGCGCCCGGcc  >  2:2282741/1‑101 (MQ=255)
                                                                                                 tCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTCGCGCCCGGcc  >  2:3542219/1‑101 (MQ=255)
                                                                                                 |                                                                                                    
TCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGCAATATTGCGGACAGGCGTTGTCGCGCCCGGCC  >  NC_014034/1699883‑1700078

Base quality scores:  ATCG  < 3 ≤  ATCG  < 24 ≤  ATCG  < 31 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG