New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 392448 =25 (0.120)26 (0.120) 12/186 25.9 0.44 intergenic (‑96/+47) RCAP_rcc00363/trmD hypothetical protein/tRNA (guanine‑N1)‑methyltransferase
?NC_014034 1732717 = 42 (0.190)intergenic (‑141/‑262) RCAP_rcc01601/tgt FeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Rejected: Position hash score below cutoff.

GAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGA.CCGGGC..GGGCG..G..G..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGG  >  NC_014034/1732624‑1732817
                                                                                                      |                                                                                                    
gagGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGAGCCGGGC..CGGGG..G..G..GCGCCGGCCccccc                                                                                                >  2:3598854/1‑101 (MQ=255)
 agGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGC..GGGGG..G..G..GCGCCCCCCcccccc                                                                                               >  1:1317387/1‑100 (MQ=255)
  ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCAccccc                                                                                                <  2:4112580/101‑1 (MQ=25)
   gTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCg                                                                                               <  1:1866188/101‑1 (MQ=25)
   gTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCg                                                                                               <  2:4150718/101‑1 (MQ=25)
    tAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCgg                                                                                              <  2:2123651/101‑1 (MQ=25)
    tAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGC..AGG.G..G..G..CCCCCCGCCccggctcccc                                                                                           >  2:2207937/1‑93 (MQ=255)
     aaaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGG..GGGGG..G..G..GCCCCCCCCcccccccccc                                                                                           >  2:2846371/1‑96 (MQ=255)
     aaaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGC..GGGGGATC..C..GCCCCCGCACCCCCGGc                                                                                             <  1:3049891/101‑1 (MQ=37)
     aaaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGC..GGGGG..G..G..CCGCCCCCCcccccccccc                                                                                           >  2:452225/1‑96 (MQ=255)
      aaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCg                                                                                            <  1:2230478/101‑1 (MQ=25)
       aGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGt                                                                                           <  2:2039214/101‑1 (MQ=25)
        ggCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTa                                                                                          <  1:4291833/101‑1 (MQ=25)
         gCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTAt                                                                                         <  1:2822342/101‑1 (MQ=25)
           cGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTAttt                                                                                       <  2:2197933/101‑1 (MQ=25)
                 tGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCa                                                                                 <  1:3599102/101‑1 (MQ=21)
                    ggggCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCaaga                                                                              <  2:4033100/101‑1 (MQ=21)
                     gggCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCaagaa                                                                             <  2:937029/101‑1 (MQ=21)
                      ggCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCaagaag                                                                            <  2:1852620/101‑1 (MQ=18)
                      ggCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCaagaag                                                                            <  2:2109448/101‑1 (MQ=18)
                        cGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAgaagaa                                                                          <  1:2161466/101‑1 (MQ=18)
                         ggAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGaaa                                                                         <  2:1299150/101‑1 (MQ=18)
                          gAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAc                                                                        <  2:2078712/101‑1 (MQ=18)
                            aaGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACaa                                                                      <  1:3469757/101‑1 (MQ=18)
                            aaGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACaa                                                                      <  1:3918651/101‑1 (MQ=18)
                            aaGCCGCGTGACAAGCAGGGTTTCGCCGCGGCGGGGCAGCGGCGACCCGGGCAGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACaa                                                                      <  2:1058920/101‑1 (MQ=18)
                              gCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGc                                                                    <  1:3525803/101‑1 (MQ=18)
                              gCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGc                                                                    <  2:463680/101‑1 (MQ=18)
                                   gTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGAc                                                               <  2:3214114/101‑1 (MQ=18)
                                    tGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGAcc                                                              <  1:2781204/101‑1 (MQ=18)
                                     gCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCa                                                             <  1:3881522/101‑1 (MQ=18)
                                             gggTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGcccccc                                                     <  2:2977817/101‑1 (MQ=18)
                                              ggTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCt                                                    <  2:2144188/101‑1 (MQ=18)
                                              ggTTTCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCt                                                    <  2:3015952/101‑1 (MQ=18)
                                                tttCGCCGCGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGc                                                  <  2:1831561/101‑1 (MQ=18)
                                                 ttCGCCGCGGCGGGGCCGCGGCGACCCGGGC..GGGGT..GTCC..GCCGCCACGCCCCTGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCa                                                 <  2:3696035/101‑1 (MQ=25)
                                                 ctccccgcTGCGGGGCCGCCGCGACCCGGGCGGGGGCC..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCa                                                 <  2:2006285/97‑1 (MQ=21)
                                                    gccgcGGCGGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGcc                                              <  2:667191/101‑1 (MQ=18)
                                                       gcggggGGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCa                                           <  2:1317290/101‑1 (MQ=18)
                                                        cggcggGGCCGCGGCGACCCGGGG..GGGGG..G..GCCGCCCCCCCCCCCCCCGCGGATTTTCGCCCAGAAGAAACAAACGGGGCCAAGCCCCCCCCCACCCGccg                                          >  2:2082861/1‑99 (MQ=255)
                                                         ggcggGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATg                                         <  1:3258662/101‑1 (MQ=18)
                                                           cggGGCCGCGGCGACCCGGGCGGGGGCT..C..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGcc                                       <  1:298500/101‑1 (MQ=18)
                                                                                     cgggggctccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAgg                   <  1:2137218/91‑1 (MQ=255)
                                                                                      gggggctccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAggg                  >  1:2927042/10‑101 (MQ=255)
                                                                                      gggggctccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAACGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGGCAACAGGCAAggg                  >  1:4359490/10‑101 (MQ=25)
                                                                                        gggctccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGcc                <  1:345314/94‑1 (MQ=255)
                                                                                         ggctccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGccc               >  2:2670165/7‑101 (MQ=255)
                                                                                         ggctccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGccc               >  1:3446233/7‑101 (MQ=255)
                                                                                         ggctccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGccc               <  1:3384421/95‑1 (MQ=255)
                                                                                         ggctccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAAACCCCCCTGCCGCCGCATGCCGTCCGCGTCAACAGGCAAGGGccc               >  1:509692/7‑101 (MQ=37)
                                                                                          gctccgCCCCCGCCCCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGCCCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGcccc              <  1:2208076/96‑1 (MQ=255)
                                                                                          gctccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGcccc              <  1:2865157/96‑1 (MQ=255)
                                                                                          gctccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGcccc              >  2:833814/6‑101 (MQ=255)
                                                                                           ctccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCg             >  2:2334415/5‑101 (MQ=255)
                                                                                           ctccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCg             <  1:4119566/97‑1 (MQ=255)
                                                                                            tccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGc            >  2:2116198/4‑101 (MQ=255)
                                                                                            tccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGc            >  1:1967848/4‑101 (MQ=255)
                                                                                             ccgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCt           <  1:2889601/99‑1 (MQ=255)
                                                                                              cgCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCtt          >  2:500592/2‑101 (MQ=255)
                                                                                               gCCCCCGCCCCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCttt         <  2:4265896/101‑1 (MQ=255)
                                                                                                cccccGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCtttt        <  1:4404972/101‑1 (MQ=255)
                                                                                                 ccccGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTc       <  1:1913043/101‑1 (MQ=255)
                                                                                                  cccGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAAGCAAGGGCCCCGCTTTTcc      >  1:2528646/1‑101 (MQ=255)
                                                                                                    cGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTa    <  2:4141111/101‑1 (MQ=255)
                                                                                                     gCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAg   <  2:4359219/101‑1 (MQ=255)
                                                                                                      cACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAgg  <  2:3192492/101‑1 (MQ=255)
                                                                                                      |                                                                                                    
GAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGA.CCGGGC..GGGCG..G..G..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGG  >  NC_014034/1732624‑1732817

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG