New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_014034	2880583 =	0 (0.000)	87 (0.450)	44/182	9.4	0.93	coding (1002/1053 nt)	pdxA	4‑hydroxythreonine‑4‑phosphate dehydrogenase
	?	NC_014034	= 2880584	14 (0.070)					coding (1001/1053 nt)	pdxA	4‑hydroxythreonine‑4‑phosphate dehydrogenase
Rejected: Position hash score below cutoff.

AAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........C..CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCGA  >  NC_014034/2880485‑2880684                                                                                                                         |                                                                                                     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                                                                  >  1:3976822/1‑101 (MQ=9)                                               ttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGa                                                                  >  2:1881732/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:2563696/1‑101 (MQ=12)                                                             aaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCggg                                                    >  1:3072557/1‑101 (MQ=12)                                                              atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGc                                                   >  2:2596713/1‑101 (MQ=12)                                                                                                              gtcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCa             >  2:3079274/11‑101 (MQ=255)                                                                                                              gtcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCa             >  1:2509069/11‑101 (MQ=255)                                                                                                              gtcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCa             <  1:168847/91‑1 (MQ=255)                                                                                                               tcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCAt            <  2:3357168/92‑1 (MQ=255)                                                                                                               tcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCAt            >  2:2175488/10‑101 (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)                                                                                                                         gccgccCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCGa  <  2:4050905/101‑1 (MQ=255)                                                                                                                         gccgccCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCGa  <  1:2697990/101‑1 (MQ=255)                                                                                                                         |                                                                                                     AAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........C..CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCGA  >  NC_014034/2880485‑2880684 Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG