New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_014034	2880583 =	4 (0.020)	108 (0.500)	40/182	11.1	0.87	coding (1002/1053 nt)	pdxA	4‑hydroxythreonine‑4‑phosphate dehydrogenase
	?	NC_014034	= 2880584	27 (0.110)					coding (1001/1053 nt)	pdxA	4‑hydroxythreonine‑4‑phosphate dehydrogenase
Rejected: Position hash score below cutoff.

GGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........C..CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCGA  >  NC_014034/2880488‑2880684                                                                                                                      |                                                                                                       ggcggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcac                                                                                                    >  1:2053843/1‑97 (MQ=255) ggcggcAGCCCGTCAACCGCACCCAAAGCCCGCCCCCTTCCTTTTCATCTTTCTTCAAATATCCGGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcac                                                                                                    <  2:2985499/101‑5 (MQ=25)  gcggcAGCCCTTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATACCCGGGGGGTCCGGGGGCAGACCCCCCGGCG..........GC.........g..cgtcaca                                                                                                   <  2:4514289/101‑6 (MQ=25)  ctggcCGCCCCTCACCCGCAGCGAAGGCCCGCCCCCTACCTTTTCATCTTGCTTCAAATTTCCGGGGGGTCCGGGGGGAGAACCCCCGGCG..........GC.........g..cgtcaca                                                                                                   <  2:3098333/99‑6 (MQ=255)   cggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCCCCG..........GC.........g..cgtcacag                                                                                                  >  2:4207800/1‑95 (MQ=25)     gcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcacagcc                                                                                                >  1:1997346/1‑93 (MQ=255)      cAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........g..cgtcacagccg                                                                                               >  1:4159063/1‑92 (MQ=25)                                          ttttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCcgc                                                                    >  2:4501785/1‑101 (MQ=2)                                           tttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCcgcg                                                                   >  1:1037891/1‑101 (MQ=2)                                           tttCATCTTGCTTCAAATATCCCGGGGGCCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCcgcg                                                                   >  2:527606/1‑101 (MQ=2)                                             tCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGAt                                                                 >  2:123905/1‑101 (MQ=9)                                             tCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGAt                                                                 >  2:3495523/1‑101 (MQ=9)                                                tCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAg                                                              >  1:4722093/1‑101 (MQ=9)                                                  ttGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGGGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCCCCGCCATCGTCAAGGCCGCGATCCGAc                                                            >  2:2821630/1‑101 (MQ=9)                                                  ttGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCCTCTTCGCCGCCCTCGTCAAGGCCGCGATCAGAc                                                            >  1:2382809/1‑101 (MQ=9)                                                   tGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACt                                                           >  2:3135307/1‑101 (MQ=9)                                                   tGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACt                                                           >  2:4634201/1‑101 (MQ=9)                                                     cTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgt                                                         >  1:547718/1‑101 (MQ=12)                                                     cTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGCCAAGGCCGCGATCAGACtga                                                         >  1:1236281/1‑100 (MQ=9)                                                      ttCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGatcaga                                                        >  2:2840661/1‑96 (MQ=9)                                                       tCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGc                                                       >  1:3379932/1‑101 (MQ=12)                                                       tCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCACCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGc                                                       >  2:4729504/1‑101 (MQ=9)                                                         aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCgg                                                     >  2:3956869/1‑101 (MQ=12)                                                         aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCgg                                                     >  2:2475196/1‑101 (MQ=12)                                                         aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGGCACAGC.........CG.CGC..CGCCCCCCTCTTCGCCGCCATCGTCAAGGCCGCCATCAAACTGTGCgg                                                     >  2:3251469/1‑101 (MQ=2)                                                         aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGGCACAGC..CGCGGCGGCCCCATCTTCGCCGCCCTCGTCAAGGCCGCGCTCAGACTGCGCgg                                                     >  2:129286/1‑101 (MQ=2)                                                          aaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCggg                                                    >  1:3887630/1‑101 (MQ=12)                                                          aaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGC.........CG.CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCggg                                                    >  2:1398853/1‑101 (MQ=12)                                                           atatCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCGCGTCACAGCCGCGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGc                                                   >  1:1326942/1‑101 (MQ=12)                                                                                                           gtcacagccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCATTGTCATAGGCCGTGCCa             >  2:1654202/11‑101 (MQ=255)                                                                                                            gcaccgtcgcgc..cgcCCCCATCCTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCGGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCAt            <  1:773799/92‑1 (MQ=255)                                                                                                             cacagccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATg           <  2:1389644/93‑1 (MQ=255)                                                                                                             cacagccgcgc..cgcCCCCATCTTCGCCGCCATCCTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATg           <  1:4040707/93‑1 (MQ=255)                                                                                                              gcagccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGa          <  2:3065771/94‑1 (MQ=255)                                                                                                              acagccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGa          <  1:3060687/94‑1 (MQ=255)                                                                                                              acagccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGa          >  1:554441/8‑101 (MQ=255)                                                                                                               cagccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGAt         <  1:1555288/95‑1 (MQ=255)                                                                                                               cagccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGAt         <  2:3496094/95‑1 (MQ=255)                                                                                                               cagccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGAt         <  1:938158/95‑1 (MQ=255)                                                                                                                agccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATc        >  2:1141053/6‑101 (MQ=255)                                                                                                                agccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATc        >  1:4564696/6‑101 (MQ=255)                                                                                                                agccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATc        <  2:481709/96‑1 (MQ=255)                                                                                                                 gcctcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCg       <  2:3361049/97‑1 (MQ=255)                                                                                                                 gccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCg       <  2:2723093/97‑1 (MQ=255)                                                                                                                  ccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCgg      <  2:2574342/98‑1 (MQ=255)                                                                                                                  ccgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCgg      >  1:4141279/4‑101 (MQ=255)                                                                                                                   cgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCggg     <  1:1472135/99‑1 (MQ=255)                                                                                                                   cgcgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCggg     <  1:2257523/99‑1 (MQ=255)                                                                                                                     cgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCg   >  2:3778131/1‑101 (MQ=255)                                                                                                                     cgc..cgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCg   <  2:3006211/101‑1 (MQ=255)                                                                                                                      gc..cgccCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCGa  <  1:2255330/101‑1 (MQ=255)                                                                                                                      gc..cgccCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCGa  >  1:731672/1‑101 (MQ=255)                                                                                                                      gc..cgccCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCGa  >  1:592006/1‑101 (MQ=255)                                                                                                                      |                                                                                                       GGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GC.........C..CGC..CGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCGA  >  NC_014034/2880488‑2880684 Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG