New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 2880583 =2 (0.010)71 (0.350) 38/182 13.0 0.82 coding (1002/1053 nt) pdxA 4‑hydroxythreonine‑4‑phosphate dehydrogenase
?NC_014034 = 2880584 29 (0.130)coding (1001/1053 nt) pdxA 4‑hydroxythreonine‑4‑phosphate dehydrogenase
Rejected: Position hash score below cutoff.

CAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCC.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCG  >  NC_014034/2880484‑2880683
                                                                                                              |                                                                                                    
cAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCg.cg                                                                                                     >  2:2896521/1‑101 (MQ=255)
cAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCg.cg                                                                                                     >  2:2391693/1‑101 (MQ=255)
 aaGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCg.cgt                                                                                                    >  2:795088/1‑100 (MQ=255)
 aaGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCg.cgt                                                                                                    >  1:3201459/1‑100 (MQ=255)
   gggCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCg.cgtca                                                                                                  >  2:2185637/1‑98 (MQ=255)
    ggcggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCg.cgtcac                                                                                                 >  2:1042544/1‑97 (MQ=255)
    ggcggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCg.cgtcac                                                                                                 >  1:545642/1‑97 (MQ=255)
       ggcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCg.cgtcacagc                                                                                              >  1:1266681/1‑94 (MQ=255)
         cAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCg.cgtcacagccg                                                                                            >  2:1890414/1‑92 (MQ=25)
                                                 cATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATc                                                               >  2:2313433/1‑101 (MQ=9)
                                                 cATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATc                                                               >  2:2126126/1‑101 (MQ=9)
                                                 cATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATc                                                               >  2:3552354/1‑101 (MQ=9)
                                                 cATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGGCAAGGCCGCGATc                                                               >  2:3929105/1‑101 (MQ=2)
                                                  ctCTTCCTTCAAATACCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCa                                                              <  2:2120557/100‑1 (MQ=9)
                                                  aTCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCa                                                              >  2:2872736/1‑101 (MQ=9)
                                                  aTCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCa                                                              >  1:1486681/1‑101 (MQ=9)
                                                   tCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATccg                                                             >  2:2483284/1‑99 (MQ=9)
                                                    cTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGa                                                            >  1:2732462/1‑101 (MQ=9)
                                                    cTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGa                                                            >  2:1797565/1‑101 (MQ=9)
                                                     ttGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGAc                                                           >  2:2180277/1‑101 (MQ=9)
                                                     ttGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGAc                                                           >  1:3616674/1‑101 (MQ=9)
                                                       gCTTCAAATATCCGGGGGGTCGGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtg                                                         <  1:2185776/101‑1 (MQ=12)
                                                        cTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgt                                                        >  1:2041846/1‑101 (MQ=12)
                                                        cTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgt                                                        >  2:3713978/1‑101 (MQ=12)
                                                         ttCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgtg                                                       >  1:704137/1‑101 (MQ=12)
                                                         ttCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtctg                                                       >  2:1446721/1‑101 (MQ=9)
                                                          tCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCCTCGTCAAGGCCGCGATCAGACTGTGc                                                      >  1:1342240/1‑101 (MQ=9)
                                                           cAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGGCGCGATCAGACTGTGCg                                                     >  2:4010299/1‑101 (MQ=9)
                                                           cAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCg                                                     >  1:2577397/1‑101 (MQ=12)
                                                            aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCgg                                                    >  1:1920837/1‑101 (MQ=12)
                                                            aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCgg                                                    >  1:2788940/1‑101 (MQ=12)
                                                             aaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCggg                                                   >  1:693873/1‑101 (MQ=12)
                                                                                                    gtcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGGTCAGACTGTGCGGGCTGGCCTTGCCGGTGCTCGCAATGTCATAGGCCGTGCCa            >  2:738042/11‑101 (MQ=255)
                                                                                                    gtcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCa            >  1:3962463/11‑101 (MQ=255)
                                                                                                    gtcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCa            <  2:733111/91‑1 (MQ=255)
                                                                                                    gtcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCa            >  1:1323209/11‑101 (MQ=255)
                                                                                                     tcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCAt           <  1:2651174/92‑1 (MQ=255)
                                                                                                     tcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCAt           >  2:615361/10‑101 (MQ=255)
                                                                                                      cacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATg          <  2:1314101/93‑1 (MQ=255)
                                                                                                      cacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATg          <  1:4279224/93‑1 (MQ=255)
                                                                                                      cacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATg          >  1:3720755/9‑101 (MQ=255)
                                                                                                      cacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATg          >  2:760091/9‑101 (MQ=255)
                                                                                                       acagccgcgccgcCCCCATCTTCGCCGCTATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGa         >  1:2666075/8‑101 (MQ=255)
                                                                                                       acagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGa         >  1:1116793/8‑101 (MQ=255)
                                                                                                       acagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGa         >  1:604210/8‑101 (MQ=255)
                                                                                                       acagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGa         >  1:3622500/8‑101 (MQ=255)
                                                                                                        cagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGAt        >  1:4377941/7‑101 (MQ=255)
                                                                                                        cagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGAt        <  2:651285/95‑1 (MQ=255)
                                                                                                        cagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGAt        >  2:737644/7‑101 (MQ=255)
                                                                                                        cagccgcgccgcCCCCATCTACGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCAGGTGCCCGCAATGTCATAGGCCGTGCCATGAt        >  2:315174/7‑101 (MQ=255)
                                                                                                         agccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATc       >  2:2994822/6‑101 (MQ=255)
                                                                                                         agccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATc       <  1:1919893/96‑1 (MQ=255)
                                                                                                         agccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCGGGCCTTGCCGGAGCCCGCATTGTCATAGGCCGGGCCATGATc       >  2:1762576/6‑101 (MQ=255)
                                                                                                         agccgcgccgcCCCCATCTTCGCCCCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATc       >  1:4440299/6‑101 (MQ=255)
                                                                                                          gccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCg      >  2:3357095/5‑101 (MQ=255)
                                                                                                          gccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCg      >  1:998403/5‑101 (MQ=255)
                                                                                                           ccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCgg     >  2:2271686/4‑101 (MQ=255)
                                                                                                            cgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCggg    >  1:2928546/3‑101 (MQ=255)
                                                                                                            cgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCggg    >  1:3135495/3‑101 (MQ=255)
                                                                                                             gcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGc   >  1:338041/2‑101 (MQ=255)
                                                                                                             gcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGc   >  1:1236383/2‑101 (MQ=255)
                                                                                                              cgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCg  >  2:2457304/1‑101 (MQ=255)
                                                                                                              |                                                                                                    
CAAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCC.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCG  >  NC_014034/2880484‑2880683

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG