New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 = 39246821 (0.110)55 (0.300) 32/174 12.5 0.84 intergenic (‑116/+27) RCAP_rcc00363/trmD hypothetical protein/tRNA (guanine‑N1)‑methyltransferase
?NC_014034 2880583 = 0 (0.000)coding (1002/1053 nt) pdxA 4‑hydroxythreonine‑4‑phosphate dehydrogenase
Rejected: Position hash score below cutoff.

AAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........C..CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCG  >  NC_014034/2880485‑2880683
                                                                                                                       |                                                                                                    
aaGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgt                                                                                                    >  2:4033388/1‑100 (MQ=255)
aaGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgg                                                                                                    >  2:2860165/1‑100 (MQ=255)
    gcggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cggcaca                                                                                                >  1:4144279/1‑96 (MQ=255)
        cAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcacagccg                                                                                            >  2:452169/1‑92 (MQ=25)
                                            ttttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCcgc                                                                   >  2:3435108/1‑101 (MQ=2)
                                              ttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGa                                                                 >  2:1881732/1‑101 (MQ=9)
                                              ttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGa                                                                 >  1:3976822/1‑101 (MQ=9)
                                                   cTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGa                                                            >  2:4056388/1‑101 (MQ=9)
                                                       cTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgt                                                        >  2:1285584/1‑101 (MQ=12)
                                                       cTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgg                                                        >  2:1522827/1‑100 (MQ=9)
                                                           aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCgg                                                    >  1:1396561/1‑101 (MQ=12)
                                                           aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCgg                                                    >  2:2697784/1‑101 (MQ=12)
                                                           aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCgg                                                    >  2:231165/1‑101 (MQ=12)
                                                            aaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCggg                                                   >  1:3072557/1‑101 (MQ=12)
                                                            aaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCggg                                                   >  1:2563696/1‑101 (MQ=12)
                                                             atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGc                                                  >  2:2596713/1‑101 (MQ=12)
                                                                                                             gtcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCa            >  1:2509069/11‑101 (MQ=255)
                                                                                                             gtcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCa            <  1:168847/91‑1 (MQ=255)
                                                                                                             gtcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCa            >  2:3079274/11‑101 (MQ=255)
                                                                                                              tcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCAt           >  2:2175488/10‑101 (MQ=255)
                                                                                                              tcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCAt           <  2:3357168/92‑1 (MQ=255)
                                                                                                               cacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATg          <  2:4144017/93‑1 (MQ=255)
                                                                                                                ccagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGa         <  2:3098356/94‑1 (MQ=255)
                                                                                                                 cagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGAt        <  2:3386951/95‑1 (MQ=255)
                                                                                                                  agccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATc       <  2:3689305/96‑1 (MQ=255)
                                                                                                                  agccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATc       >  1:3078262/6‑101 (MQ=255)
                                                                                                                  agccgcgccccCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATc       <  2:822651/96‑1 (MQ=255)
                                                                                                                    ccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCgg     <  2:2783971/98‑1 (MQ=255)
                                                                                                                       cgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCg  >  2:1483640/1‑101 (MQ=255)
                                                                                                                       cgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCg  <  1:448025/101‑1 (MQ=255)
                                                                                                                       cgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCg  >  1:1344502/1‑101 (MQ=255)
                                                                                                                       |                                                                                                    
AAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........C..CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCG  >  NC_014034/2880485‑2880683

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