New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 318460 =125 (0.570)33 (0.150) 7/194 30.5 0.29 intergenic (+88/+67) rpsO/RCAP_rcc00282 30S ribosomal protein S15/aldehyde dehydrogenase
?NC_014034 1253898 = 35 (0.150)intergenic (+56/+35) RCAP_rcc01181/hisH hypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Rejected: Position hash score below cutoff.

GGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGAG.........CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGT  >  NC_014034/1253833‑1253997
                                                                               |                                                                                                   
ggCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCg                                                                           >  1:275792/1‑101 (MQ=12)
  cGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaa                                                                         >  2:4104994/1‑101 (MQ=12)
   gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaaa                                                                        >  2:627028/1‑101 (MQ=12)
      aTAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATcc                                                                     >  2:4045528/1‑101 (MQ=9)
        aaaCCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGCGGCGAAATCctc                                                                   >  2:103769/1‑101 (MQ=9)
        aaaCCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGCCGCCTCGGCGAAAGCCGAGGCGAAATCctc                                                                   >  1:2395901/1‑101 (MQ=9)
         aaCCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCctct                                                                  >  2:3266944/1‑101 (MQ=9)
         aaCCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCctct                                                                  >  2:3242456/1‑101 (MQ=9)
          aCCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCtctc                                                                 >  1:1827775/1‑101 (MQ=9)
            cGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCCGCGAAAGCCGAGGCGAAATCCTCCCAc                                                               >  1:2325378/1‑101 (MQ=2)
             gATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACg                                                              >  2:1324756/1‑101 (MQ=9)
              aTCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACgg                                                             >  2:425547/1‑101 (MQ=9)
              aTCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACgg                                                             >  2:2454519/1‑101 (MQ=9)
                cccgccccgcTCTCATCTTGCCCTAAATCCCTGGGGGGTGGG.....GGGGCGGACCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCt                                                           <  2:3004441/101‑1 (MQ=9)
                cccgccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCt                                                           >  2:2511506/1‑101 (MQ=9)
                  cgccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTc                                                         >  1:2709680/1‑101 (MQ=9)
                  cgccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTc                                                         >  1:1068876/1‑101 (MQ=9)
                   gccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTcc                                                        >  1:2598650/1‑101 (MQ=9)
                    ccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTGCa                                                       >  2:1015137/1‑101 (MQ=2)
                    ccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCa                                                       >  1:1164426/1‑101 (MQ=9)
                     cccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAg                                                      >  1:388773/1‑101 (MQ=9)
                       cgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGtt                                                    >  1:3429484/1‑101 (MQ=9)
                          tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAg                                                 >  1:2908617/1‑101 (MQ=9)
                           ctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAgg                                                >  1:1722335/1‑101 (MQ=9)
                           ctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAgg                                                >  1:3149085/1‑101 (MQ=9)
                            tcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGa                                               >  1:4357286/1‑101 (MQ=9)
                              aTCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGaaa                                             >  1:729813/1‑101 (MQ=9)
                               tCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCC.....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAAt                                            >  1:3488224/1‑101 (MQ=9)
                                                                  ggcatctcctcgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCg              <  2:579824/91‑1 (MQ=255)
                                                                   cggagcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGa             >  2:2647070/10‑101 (MQ=255)
                                                                   ccccccccgggcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGa             <  1:2399965/91‑1 (MQ=255)
                                                                    ggcccccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAt            <  1:4331141/93‑1 (MQ=255)
                                                                    ggagcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAt            >  2:1639168/9‑101 (MQ=255)
                                                                    ggagcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAt            >  2:764841/9‑101 (MQ=255)
                                                                     gagcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAtt           >  1:3383119/8‑101 (MQ=255)
                                                                      gccccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAttt          <  2:1451238/95‑1 (MQ=255)
                                                                      agcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAttt          >  2:1863516/7‑101 (MQ=255)
                                                                      agcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAttt          >  2:2231343/7‑101 (MQ=255)
                                                                       gcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtc         >  2:3051593/6‑101 (MQ=255)
                                                                        cccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtct        <  1:1656276/97‑1 (MQ=255)
                                                                        cccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtcc        >  1:3323764/5‑100 (MQ=255)
                                                                         ccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtctc       >  2:3016736/4‑101 (MQ=255)
                                                                         ccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtctc       >  2:2001306/4‑101 (MQ=255)
                                                                         ccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtctc       >  2:1961816/4‑101 (MQ=255)
                                                                         ccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtctc       >  2:1723697/4‑101 (MQ=255)
                                                                           ccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCg     >  2:1200690/2‑101 (MQ=255)
                                                                            cgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCgg    >  1:2942023/1‑101 (MQ=255)
                                                                              ccGGGTCGCCTCGGCGAAAGCCGAGGCTGAATCCTCTCACGGCTTCCAGTTGAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGt  <  1:735073/101‑1 (MQ=255)
                                                                              ccGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCgggg  >  2:767233/1‑100 (MQ=255)
                                                                              ccGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGt  >  2:3238078/1‑101 (MQ=255)
                                                                              ccGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGt  >  1:19171/1‑101 (MQ=255)
                                                                               |                                                                                                   
GGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGAG.........CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCCGGGT  >  NC_014034/1253833‑1253997

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG