New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 1000980 =381 (0.880)77 (0.180) 2/188 41.4 0.29 intergenic (+24/+99) bhbA/bfr methylmalonyl‑CoA mutase/bacterioferritin
?NC_014034 = 2600943 1 (0.000)intergenic (+5/‑136) RCAP_rcc02424/surE short‑chain dehydrogenase/reductase family oxidoreductase/5'‑nucleotidase SurE
Rejected: Position hash score below cutoff.

GCCACCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACGGGGGGGGCCCCTGCC  >  NC_014034/2600843‑2600959
                                                                                                    |                
gCCACCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTgacggac                  <  2:1723917/101‑1 (MQ=255)
gCCACCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTgacagac                  >  1:4990749/1‑101 (MQ=255)
 ccACCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTgacggacg                 <  2:9177720/101‑1 (MQ=255)
 ccACCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTgacggacg                 >  2:7179372/1‑101 (MQ=255)
 ccACCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTgacggacg                 <  2:5421868/101‑1 (MQ=255)
 ccACCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTgacggacg                 <  1:3414847/101‑1 (MQ=255)
 ccACCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTgacggacg                 <  1:5665856/101‑1 (MQ=255)
  cACCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGacggacgg                >  2:4428334/1‑101 (MQ=255)
  cACCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGacggacgg                >  2:2124543/1‑101 (MQ=255)
  cACCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGacggacgg                <  1:1498104/101‑1 (MQ=255)
   cccgccgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggg               <  2:3591086/100‑1 (MQ=255)
   aCCGCCGTGCGGGCGGGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggg               >  1:7261521/1‑101 (MQ=255)
   aCCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggg               <  2:2052771/101‑1 (MQ=255)
   aCCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggg               >  1:4774714/1‑101 (MQ=255)
   aCCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggg               >  1:3115199/1‑101 (MQ=255)
   aCCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggg               >  1:9117853/1‑101 (MQ=255)
   aCCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggg               >  2:3415685/1‑101 (MQ=255)
      gccgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACgggggg            >  1:4917899/1‑101 (MQ=255)
      gccgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACgggggg            >  1:544550/1‑101 (MQ=255)
       ccgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggc           <  2:6986121/101‑2 (MQ=255)
       ccgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggc           <  1:8445940/101‑2 (MQ=255)
       ccgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggc           <  1:7585622/101‑2 (MQ=255)
        cgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggct          <  1:2785286/101‑3 (MQ=255)
        cgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggct          >  2:4772332/1‑99 (MQ=255)
        cgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggct          >  2:4795426/1‑99 (MQ=255)
        cgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggct          >  2:7627603/1‑99 (MQ=255)
        cgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggct          <  2:7673232/101‑3 (MQ=255)
        cgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggct          <  2:7953498/101‑3 (MQ=255)
        cgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggct          <  1:1431109/101‑3 (MQ=255)
        cgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggct          <  1:8651200/101‑3 (MQ=255)
        cgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggct          <  1:6271234/101‑3 (MQ=255)
        cgTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACcggggcgc          >  2:4204478/1‑98 (MQ=255)
         gTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctc         <  2:8421623/101‑4 (MQ=255)
         gTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctc         >  1:5877879/1‑98 (MQ=255)
         gTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctc         >  2:4882172/1‑98 (MQ=255)
         gTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACGGGGGGGGc         >  2:75949/1‑101 (MQ=255)
         gTGCGGGCGAGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGACGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctc         >  2:5954189/1‑98 (MQ=255)
          tGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctct        <  2:3207439/101‑5 (MQ=255)
          tGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctct        <  2:3559821/101‑5 (MQ=255)
          tGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctct        <  2:3749003/101‑5 (MQ=255)
          tGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctct        <  1:6311054/101‑5 (MQ=255)
          tGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctct        <  1:1655144/101‑5 (MQ=255)
           gCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctg       >  1:8425314/1‑96 (MQ=255)
            cGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgc      <  1:3773030/101‑7 (MQ=255)
            cGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgc      >  2:1004352/1‑95 (MQ=255)
             gggCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgcc     <  2:621595/101‑8 (MQ=255)
             gggCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgcc     >  1:5956494/1‑94 (MQ=255)
             gggCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgcc     >  1:2590293/1‑94 (MQ=255)
             gggCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgcc     >  2:7772175/1‑94 (MQ=255)
             gggCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgcc     <  2:8807838/101‑8 (MQ=255)
             gggCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgcc     >  1:472759/1‑94 (MQ=255)
              ggCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgccc    >  1:4464040/1‑93 (MQ=255)
               gcgcgcTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgcccc   <  2:7086638/101‑10 (MQ=255)
               gcgcgcTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgcccc   <  2:7381564/101‑10 (MQ=255)
               gcgcgcTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgcccc   <  2:4107371/101‑10 (MQ=255)
               gcgcgcTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgcccc   <  2:3901132/101‑10 (MQ=255)
                cgcgcTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgccccc  <  1:5002909/101‑11 (MQ=255)
                cgcgcTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACggggggctctgccccc  <  2:7828159/101‑11 (MQ=255)
                                                                                                    |                
GCCACCGCCGTGCGGGCGCGCTTCTTCCCGTCCGAGGACCGCACGCCGCTGATCGATTGCCGCGACGAAGCGAAACGCCTTCTTGCCGAGATCTGACGGACGGGGGGGGCCCCTGCC  >  NC_014034/2600843‑2600959

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