New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 39246233 (0.140)26 (0.110) 2/198 38.6 0.43 intergenic (‑110/+33) RCAP_rcc00363/trmD hypothetical protein/tRNA (guanine‑N1)‑methyltransferase
?NC_014034 1253896 = 37 (0.160)intergenic (+54/+37) RCAP_rcc01181/hisH hypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Rejected: Position hash score below cutoff.

TGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGAG....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCC  >  NC_014034/1253796‑1253993
                                                                                                             |                                                                                                 
tGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cg                                                                                                   >  2:826320/1‑101 (MQ=255)
cgggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cg                                                                                                   <  2:3855967/100‑1 (MQ=255)
 gggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgg                                                                                                  >  1:3535510/1‑100 (MQ=255)
  ggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgga                                                                                                 <  1:3904558/101‑3 (MQ=255)
  ggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgga                                                                                                 <  2:4385871/101‑3 (MQ=255)
  ggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgga                                                                                                 >  2:1550662/1‑99 (MQ=255)
   gCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggag                                                                                                <  2:391647/101‑4 (MQ=255)
    taaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagc                                                                                               <  1:2121340/100‑5 (MQ=255)
     aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcc                                                                                              <  2:4255640/101‑6 (MQ=255)
     aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcc                                                                                              >  2:3299303/1‑96 (MQ=255)
     aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcc                                                                                              <  2:2834656/101‑6 (MQ=255)
     aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagca                                                                                              >  2:2121190/1‑96 (MQ=255)
      aGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccc                                                                                             <  2:2022593/101‑7 (MQ=255)
      aGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccc                                                                                             >  2:1028516/1‑95 (MQ=255)
        cTCGGCTTCCAGCCCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccccc                                                                                           <  2:3792155/101‑9 (MQ=38)
        cTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccccc                                                                                           >  1:2540634/1‑93 (MQ=255)
        cTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccccc                                                                                           <  2:4254239/101‑9 (MQ=255)
         tCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcccccg                                                                                          >  2:3379154/1‑92 (MQ=38)
             cTTCCAGCCCGGCTGCGCCTTGGAGGCGAAATAACCCGATCCCGCCCCGCTCCCATCTTCCCCAAAATCCCCCGGGGGTGCG.....GGGG....CGGCgcccccgccgg                                                                                      <  2:2902753/101‑11 (MQ=19)
                                     ggCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCg                                                                       >  1:783322/1‑101 (MQ=12)
                                        gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaaa                                                                    >  2:2113155/1‑101 (MQ=12)
                                        gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaaa                                                                    >  1:2234569/1‑101 (MQ=9)
                                        gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaaa                                                                    >  2:3503396/1‑101 (MQ=12)
                                        gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGGGCGGGGGCGGAGCCCCCGCCGGGTAGCCTCGGCGAAAGCCGAGGCGaaa                                                                    >  2:2855155/1‑101 (MQ=9)
                                         aaaTAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAAt                                                                   >  1:3902527/1‑101 (MQ=9)
                                           aTAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATcc                                                                 >  2:3645734/1‑101 (MQ=9)
                                            tAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCct                                                                >  1:4270036/1‑101 (MQ=9)
                                            tAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCct                                                                >  1:4383723/1‑101 (MQ=9)
                                            tAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCct                                                                >  1:3380305/1‑101 (MQ=9)
                                              aaCCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCctct                                                              >  2:2745360/1‑101 (MQ=9)
                                                ccGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCa                                                            >  2:128056/1‑101 (MQ=9)
                                                 cGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCAc                                                           >  2:2252298/1‑101 (MQ=9)
                                                 cGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGAGGAGCCCCCGCCGGGTCGCCGCGGCGAAAGCAGCGGCGAAATCCTCTCAc                                                           >  1:1160051/1‑101 (MQ=12)
                                                  gATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACg                                                          >  1:2968515/1‑101 (MQ=9)
                                                  gATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGGCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCCCg                                                          >  2:2784583/1‑101 (MQ=2)
                                                   aTCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACgg                                                         >  2:3473969/1‑101 (MQ=9)
                                                   aTCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACgg                                                         >  2:2940972/1‑101 (MQ=9)
                                                    tCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGc                                                        >  1:776976/1‑101 (MQ=9)
                                                    tCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGc                                                        >  1:4090564/1‑101 (MQ=9)
                                                     cccgccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCt                                                       >  2:2918898/1‑101 (MQ=9)
                                                        gccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTcc                                                    >  2:621781/1‑101 (MQ=9)
                                                        gccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTcc                                                    >  2:1351082/1‑101 (MQ=9)
                                                         ccccgcTCTCATCTTCCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCCCGGCTTCCa                                                   >  1:1748467/1‑101 (MQ=9)
                                                         ccacgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAGGCCGAGGCGTAATCCTCTCACGGCTTccc                                                   >  1:2658816/4‑100 (MQ=2)
                                                          cccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAg                                                  >  1:13246/1‑101 (MQ=9)
                                                            cgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGtt                                                >  2:4562557/1‑101 (MQ=9)
                                                            cgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGtt                                                >  1:589026/1‑101 (MQ=9)
                                                             gcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTc                                               >  1:1072279/1‑101 (MQ=9)
                                                              ctctCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCCCGGCTTCCAGTTCa                                              >  2:1374257/1‑101 (MQ=9)
                                                               tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAg                                             >  2:3098563/1‑101 (MQ=9)
                                                               tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAg                                             >  1:3239096/1‑101 (MQ=9)
                                                                 tcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGa                                           >  1:1852214/1‑101 (MQ=9)
                                                                 tcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGa                                           >  2:1546584/1‑101 (MQ=9)
                                                                   aTCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGaaa                                         >  2:2896212/1‑101 (MQ=9)
                                                                    tCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAAt                                        >  2:2188324/1‑101 (MQ=9)
                                                                                                    ggagcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAt        >  2:3873057/9‑101 (MQ=255)
                                                                                                    ggagcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAt        >  2:2257756/9‑101 (MQ=255)
                                                                                                     gcgccccggcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAtt       <  1:2919127/93‑1 (MQ=255)
                                                                                                      agcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAttt      >  2:1102081/7‑101 (MQ=255)
                                                                                                      agcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGAttt      >  1:523554/7‑101 (MQ=255)
                                                                                                       gcccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtc     >  1:778416/6‑101 (MQ=255)
                                                                                                        cccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtct    >  2:4127384/5‑101 (MQ=255)
                                                                                                        cccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtct    >  2:4169726/5‑101 (MQ=255)
                                                                                                         ccccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTtctc   >  1:3548841/4‑101 (MQ=255)
                                                                                                          cccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTcc  <  1:693072/99‑1 (MQ=255)
                                                                                                          cccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTcc  <  1:3177504/99‑1 (MQ=255)
                                                                                                          cccgcCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTcc  <  1:2262720/99‑1 (MQ=255)
                                                                                                             |                                                                                                 
TGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGAG....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAATTCGAGATGATCCGCAGCCCGGTGCGCTGCGATTTCTCC  >  NC_014034/1253796‑1253993

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