New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 788008 =125 (0.560)59 (0.260) 29/188 15.2 0.49 intergenic (+80/+56) pntB/RCAP_rcc00717 pyridine nucleotide transhydrogenase subunit beta/hypothetical protein
?NC_014034 = 3701428 NA (NA)intergenic (+5/+44) grpE/mutS GrpE protein/DNA mismatch repair protein MutS
Rejected: Position hash score below cutoff.

GGTGACATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGG.......G.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGAAGGGCGCGCAGGGCGTAAAG  >  NC_014034/3701328‑3701504
                                                                                                    |                                                                                          
ggTGACATCATTGAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcag                                                                                            >  1:2566335/1‑101 (MQ=255)
ggTGACATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGTCGTGAGTTCGATGCCGAAAGCCTGaggcag                                                                                            <  1:4227051/101‑1 (MQ=255)
ggTGACATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcag                                                                                            <  2:3036955/101‑1 (MQ=255)
 gTGACATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcagg                                                                                           >  1:2074299/1‑101 (MQ=255)
 gTGACATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcagg                                                                                           >  1:4091984/1‑101 (MQ=255)
  tGACATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcaggg                                                                                          >  2:2173870/1‑100 (MQ=255)
  tGACATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcaggg                                                                                          >  2:2397842/1‑100 (MQ=255)
   gACATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcagggg                                                                                         >  2:379667/1‑99 (MQ=255)
   gACATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcagggg                                                                                         >  1:496890/1‑99 (MQ=255)
   gACATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGcc                                                                                         <  2:1575652/101‑1 (MQ=255)
   gACATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGcc                                                                                         <  2:3179486/101‑1 (MQ=255)
    aCATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCg                                                                                        <  2:343347/101‑1 (MQ=255)
    aCATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCg                                                                                        <  2:3383217/101‑1 (MQ=255)
     catcatTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcagggggg                                                                                       >  1:2267633/1‑97 (MQ=255)
     catcatTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCgg                                                                                       <  2:1200766/101‑1 (MQ=255)
     catcatTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCgg                                                                                       <  1:859169/101‑1 (MQ=255)
      atcatTCAGGTGATGGACGACGGCTTCCTGGGGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcagggggcg                                                                                      >  1:4227407/1‑96 (MQ=255)
      atcatTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcaggggggg                                                                                      >  2:2757119/1‑96 (MQ=255)
         atTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcagggggggggg                                                                                   >  1:4429885/1‑93 (MQ=255)
         atTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcagggggggggg                                                                                   >  2:4609191/1‑93 (MQ=255)
         atTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCgg.......gggc                                                                            <  1:1077816/101‑4 (MQ=255)
          ttCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcaggggggggggg                                                                                  >  1:4630088/1‑92 (MQ=255)
          ttCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcaggggggggggg                                                                                  >  2:1170659/1‑92 (MQ=255)
          ttCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGaggcaggggggggggg                                                                                  >  2:2945920/1‑92 (MQ=255)
           tCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCgg.......gggctt                                                                          <  2:2516687/101‑6 (MQ=255)
            cAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCgg.......gggcttc                                                                         <  2:1453788/101‑7 (MQ=255)
            cAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCgg.......gggcttc                                                                         <  2:1025309/101‑7 (MQ=255)
               gTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCgg.......gggcttcgcc                                                                      <  1:3297190/101‑10 (MQ=255)
                tGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCgg.......gggcttcgccc                                                                     <  2:1234503/101‑11 (MQ=255)
                tGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCgg.......gggcttcgccc                                                                     <  1:180029/101‑11 (MQ=255)
                                    ctgcACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATAtttt                                                 <  2:3130132/101‑1 (MQ=21)
                                    ctgcACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATAtttt                                                 <  1:2576188/101‑1 (MQ=21)
                                     tgcACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTc                                                <  2:3662048/101‑1 (MQ=21)
                                      gcACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCa                                               <  1:48815/101‑1 (MQ=21)
                                         cGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCa                                            <  2:2074178/101‑1 (MQ=21)
                                         cGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCa                                            <  1:3000666/101‑1 (MQ=21)
                                          gACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCaa                                           <  2:1159207/101‑1 (MQ=21)
                                           aCCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAg                                          <  2:2103242/101‑1 (MQ=21)
                                            ccGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGa                                         <  1:4428370/101‑1 (MQ=21)
                                             cGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGAt                                        <  1:1203388/101‑1 (MQ=21)
                                             cGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGAt                                        <  1:4609444/101‑1 (MQ=21)
                                                ctgctgCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATaaa                                     <  1:16531/101‑1 (MQ=21)
                                                ctgctgCGTCCCGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATaaa                                     <  1:1930486/101‑1 (MQ=21)
                                                   ctgcGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGAt                                  <  2:976552/101‑1 (MQ=21)
                                                   ctgcGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGAt                                  <  2:3819192/101‑1 (MQ=21)
                                                     gcGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCa                                <  2:3945348/101‑1 (MQ=21)
                                                     gcGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCa                                <  1:4662126/101‑1 (MQ=21)
                                                      cGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAg                               <  2:4429635/101‑1 (MQ=21)
                                                      cGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAg                               <  2:2201751/101‑1 (MQ=21)
                                                       gTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGt                              <  2:4457076/101‑1 (MQ=21)
                                                       gTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGt                              <  2:1099706/101‑1 (MQ=21)
                                                            gCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTc                         <  1:752606/101‑1 (MQ=21)
                                                            gCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTc                         <  1:2311689/101‑1 (MQ=21)
                                                            gCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGG.......GGGCTTCTCCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTc                         <  2:774777/101‑1 (MQ=25)
                                                             ggAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCa                        <  2:1683700/100‑1 (MQ=21)
                                                              gAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAg                       <  2:3282769/101‑1 (MQ=21)
                                                               aaaGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGa                      <  1:1111416/101‑1 (MQ=21)
                                                               aaaGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGa                      <  2:2295580/101‑1 (MQ=21)
                                                                   tCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGAAggg                  <  1:2461314/101‑1 (MQ=21)
                                                                     ggCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGAAGGgcg                <  2:4361228/101‑1 (MQ=21)
                                                                     ggCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGAAGGgcg                <  2:4657781/101‑1 (MQ=21)
                                                                       cGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGAAGGgcgcg              <  1:1318645/101‑1 (MQ=21)
                                                                               cGATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGAAGGGCGCGCAGGGCGt      <  2:4629835/101‑1 (MQ=21)
                                                                                gATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGAAGGGCGCGCAGGGCGTa     <  2:1434016/101‑1 (MQ=21)
                                                                                gATGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGAAGGGCGCGCAGGGCGTa     <  2:477106/101‑1 (MQ=21)
                                                                                 aTGCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGAAGGGCGCGCAGGGCGTaa    <  1:489742/101‑1 (MQ=21)
                                                                                 aTGCCGAAAGCCTGAGGCAGGCCGGGGCCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGAAGGGCGCGCAGGGCGTaa    <  1:3934272/101‑1 (MQ=21)
                                                                                   gCCGAAAGCCTGAGGCAGGCCGGGGGCTTCG.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGAAGGGCGCGCAGGGCGTAAAg  <  2:4247543/101‑1 (MQ=21)
                                                                                                    |                                                                                          
GGTGACATCATTCAGGTGATGGACGACGGCTTCCTGCTGCACGACCGTCTGCTGCGTCCGGCGAAAGTCGGCGTGAGTTCGATGCCGAAAGCCTGAGGCAGGCCGG.......G.......CCCCCGGACCCCCAGGATATTTTCAGCAAGATAAAGATCAGTCTTTCAGAAGGGCGCGCAGGGCGTAAAG  >  NC_014034/3701328‑3701504

Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG