New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_014034 2880583 =0 (0.000)77 (0.360) 44/182 12.3 0.91 coding (1002/1053 nt) pdxA 4‑hydroxythreonine‑4‑phosphate dehydrogenase
?NC_014034 = 2880584 15 (0.060)coding (1001/1053 nt) pdxA 4‑hydroxythreonine‑4‑phosphate dehydrogenase
Rejected: Position hash score below cutoff.

AAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCC.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCGA  >  NC_014034/2880485‑2880684
                                                                                                              |                                                                                                    
aaGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCg.cgt                                                                                                     >  1:4127771/1‑100 (MQ=255)
  gggCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCg.cgtca                                                                                                   >  2:1439754/1‑98 (MQ=255)
  gggCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCC..........GCg.cggcc                                                                                                   >  2:3221469/1‑98 (MQ=255)
       gcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCg.cgtcacagcc                                                                                              >  1:194009/1‑93 (MQ=255)
       gcAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCg.cgtcacagcc                                                                                              >  1:276504/1‑93 (MQ=255)
        cAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCC..........GCg.cggccccgccc                                                                                             >  1:1624015/1‑92 (MQ=25)
                                             tttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCcgcg                                                                   >  2:2947242/1‑101 (MQ=2)
                                              ttCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGa                                                                  >  1:4091452/1‑101 (MQ=9)
                                               tCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGAt                                                                 >  2:2107955/1‑101 (MQ=9)
                                                 aTCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCa                                                               >  2:4155814/1‑101 (MQ=9)
                                                   cTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCCTCTTCGCCGCCATCGTCAAGGCCGGGATCAGa                                                             >  2:4469453/1‑101 (MQ=2)
                                                    ttGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCCCGATCAGAc                                                            >  2:1138059/1‑101 (MQ=9)
                                                     tGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACt                                                           >  1:1559229/1‑101 (MQ=9)
                                                        ttCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACtgtg                                                        >  1:549121/1‑101 (MQ=12)
                                                         tCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGc                                                       >  1:291401/1‑101 (MQ=12)
                                                          cAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCg                                                      >  1:886195/1‑101 (MQ=12)
                                                           aaaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCgg                                                     >  2:2213265/1‑101 (MQ=12)
                                                            aaTATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCGGCGCGTCACAGCCGCGCCGCCCCCCTCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCggg                                                    >  1:428392/1‑101 (MQ=12)
                                                                                                   gtcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCa             >  2:3926327/11‑101 (MQ=255)
                                                                                                   gtcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCa             >  2:1822218/11‑101 (MQ=255)
                                                                                                   gtcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCa             <  2:4331255/91‑1 (MQ=255)
                                                                                                    tcacagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCAt            <  1:3968534/92‑1 (MQ=255)
                                                                                                      acagccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGa          >  2:1821586/8‑101 (MQ=255)
                                                                                                        agccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATc        <  1:4469748/96‑1 (MQ=255)
                                                                                                         gccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCg       <  2:367633/97‑1 (MQ=255)
                                                                                                          ccgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCgg      >  1:2921637/4‑101 (MQ=255)
                                                                                                           cgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCggg     >  2:683388/3‑101 (MQ=255)
                                                                                                           cgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCggg     >  1:3047942/3‑101 (MQ=255)
                                                                                                           cgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCggg     >  1:1855608/3‑101 (MQ=255)
                                                                                                           cgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCggg     >  1:1126798/3‑101 (MQ=255)
                                                                                                           cgcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCAATGATCggg     >  1:784799/3‑101 (MQ=255)
                                                                                                            gcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGc    <  2:1623910/100‑1 (MQ=255)
                                                                                                            gcgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGc    <  1:2798101/100‑1 (MQ=255)
                                                                                                             cgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCg   >  2:2155813/1‑101 (MQ=255)
                                                                                                             cgccgcCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCg   <  2:4091166/101‑1 (MQ=255)
                                                                                                              gccgccCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCGa  >  1:4583689/1‑101 (MQ=255)
                                                                                                              |                                                                                                    
AAGGGCGGCAGCCCGTCAATCGCAGCCAAAGCCCGTCCCCTTCCTTTTCATCTTGCTTCAAATATCCCGGGGGTCCGGGGGCAGAGCCCCCGGCG..........GCC.CGCCGCCCCCATCTTCGCCGCCATCGTCAAGGCCGCGATCAGACTGTGCGGGCTGGCCTTGCCGGTGCCCGCAATGTCATAGGCCGTGCCATGATCGGGCGA  >  NC_014034/2880485‑2880684

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