New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4621160 =55 (0.720)3 (0.040) 3/162 NT 5.8% coding (278/720 nt) deoD purine nucleoside phosphorylase 1; nicotinamide 1‑beta‑D‑riboside phosphorylase
?NC_000913 4621189 = 47 (0.670)coding (307/720 nt) deoD purine nucleoside phosphorylase 1; nicotinamide 1‑beta‑D‑riboside phosphorylase
Rejected: Frequency below cutoff threshold.

TAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACT  >  NC_000913/4621061‑4621258
                                                                                                   |                                                                                                  
tAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTgg                                                                                                    <  2:1618444/100‑1 (MQ=255)
tAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTgg                                                                                                    <  1:755609/100‑1 (MQ=255)
tAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTgg                                                                                                    <  1:611646/100‑1 (MQ=255)
tAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTgg                                                                                                    <  1:1644632/100‑1 (MQ=255)
caccatgtcgaagtcagcgatagcggcaaagtcatggtctttaaaacggatgcggttaactttggaatcggtgcaggcacccataccgataacgacgtcg                                                                                                    <  2:1450037‑M1/1‑1 (MQ=255)
 accatgtcgaagtcagcgatagcggcaaagtcatggtctttaaaacggatgcggttaactttggaatcggtgcaggcacccataccgataacgacgtcgc                                                                                                   <  2:1300617‑M1/2‑1 (MQ=255)
    gggTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCa                                                                                                <  2:480104/100‑1 (MQ=255)
    atgtcgaagtcagcgatagcggcaaagtcatggtctttaaaacggatgcggttaactttggaatcggtgcaggcacccataccgataacgacgtcGCGCa                                                                                                >  1:1965493‑M1/96‑100 (MQ=255)
     tgtcgaagtcagcgatagcggcaaagtcatggtctttaaaacggatgcggttaactttggaatcggtgcaggcacccataccgataacgacgtcGCGCAg                                                                                               <  1:29430‑M1/6‑1 (MQ=255)
     ggTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAg                                                                                               <  2:2026958/100‑1 (MQ=255)
     ggTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAg                                                                                               <  1:386197/100‑1 (MQ=255)
     ggTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAg                                                                                               >  1:1955982/1‑100 (MQ=255)
      gtcgaagtcagcgatagcggcaaagtcatggtctttaaaacggatgcggttaactttggaatcggtgcaggcacccataccgataacgacgtcGCGCAGt                                                                                              >  1:894416‑M1/94‑100 (MQ=255)
       tcgaagtcagcgatagcggcaaagtcatggtctttaaaacggatgcggttaactttggaatcggtgcaggcacccataccgataacgacgtcGCGCAGtt                                                                                             >  1:1277505‑M1/93‑100 (MQ=255)
        cgcGGTATGGGTAGCCCGTCCTGCTCCATCTACACCAAAGAAGTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTc                                                                                            >  1:1405410/3‑100 (MQ=255)
        cACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGtt                                                                                             >  2:433013/1‑99 (MQ=255)
           ggTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGc                                                                                         <  1:286998/100‑1 (MQ=255)
            gTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGcc                                                                                        >  1:595014/1‑100 (MQ=255)
              aTGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGc                                                                                      <  2:739413/100‑1 (MQ=255)
              aTGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGc                                                                                      <  2:1059057/100‑1 (MQ=255)
                gggTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCAc                                                                                    <  2:827806/100‑1 (MQ=255)
                gggTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCAc                                                                                    >  2:140363/1‑100 (MQ=255)
                 ggTATCCCGTCCTGCTCCATCTACACAAAAGAACTGATCACCGATTTCGGCGTGAAGAGAATTATCCGCGTGGGTTCCTGTGGCGCAGTTATGCCGCACg                                                                                   >  2:661256/1‑100 (MQ=255)
                    aTCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTaa                                                                                <  1:923084/100‑1 (MQ=255)
                     tCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTaaa                                                                               >  1:1552689/1‑100 (MQ=255)
                         gTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTg                                                                           >  1:503079/1‑100 (MQ=255)
                           ccTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTgcg                                                                         >  1:1300623/1‑100 (MQ=255)
                            cTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTgcgc                                                                        <  2:503075/100‑1 (MQ=255)
                             tctcatggtctttaaaacggatgcggttaactttggaatcggtgcaggcacccataccgataacgacgtcGCGCAGTTCTGCCGCACGTAAAACTGcgcg                                                                       >  2:1000439‑M1/71‑100 (MQ=255)
                             tGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGcgcg                                                                       >  1:319197/1‑100 (MQ=255)
                              gCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGa                                                                      >  2:1724173/1‑100 (MQ=255)
                                tCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGAcg                                                                    >  1:1348717/1‑100 (MQ=255)
                                catggtctttaaaacggatgcggttaactttggaatcggtgcaggcacccataccgataacgacgtcGCGCAGTTCTGCCGCACGTAAAACTGCGcgaga                                                                    <  1:1000442‑M1/33‑3 (MQ=255)
                                  cATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGAcgtc                                                                  >  1:778251/1‑100 (MQ=255)
                                    tCTACACGGAAGAGCTGATCACCGATTTGGGCGTGAAGAAACTTATCCGCGTGGGTTCCTGTGGCGCCGTTCTGCCGCACGTAAAACTGCGCGAcgtc                                                                  >  2:1199854/1‑98 (MQ=255)
                                     cTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGtt                                                               >  1:390842/1‑100 (MQ=255)
                                          ccAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCgg                                                          <  1:424006/100‑1 (MQ=255)
                                          ccAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCgg                                                          >  1:133143/1‑100 (MQ=255)
                                          ccAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCgg                                                          >  1:84798/1‑100 (MQ=255)
                                            aaaGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTa                                                        <  2:824617/100‑1 (MQ=255)
                                             aaGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTAt                                                       <  1:253905/100‑1 (MQ=255)
                                                aaCTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATggg                                                    <  2:366458/100‑1 (MQ=255)
                                                aaCTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATggg                                                    <  1:22183/100‑1 (MQ=255)
                                                  cTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTg                                                  >  1:1792343/1‑100 (MQ=255)
                                                    gATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGcc                                                >  2:1056775/1‑100 (MQ=255)
                                                    gATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGcc                                                >  2:548803/1‑100 (MQ=255)
                                                    gATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCACAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGcc                                                <  2:105969/100‑1 (MQ=255)
                                                     aTCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCt                                               >  1:1170402/1‑100 (MQ=255)
                                                       cACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGc                                             >  2:1168903/1‑100 (MQ=255)
                                                         aactttggaatcggtgcaggcacccataccgataacgacgtcGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCAc                                           >  1:820617‑M1/43‑100 (MQ=255)
                                                          cGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCAcc                                          <  1:182219/100‑1 (MQ=255)
                                                          cGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCAcc                                          <  1:418412/100‑1 (MQ=255)
                                                               tCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTc                                     <  1:1116986/100‑1 (MQ=255)
                                                               tCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCAGGCACCCAtac                                     <  2:522389/100‑3 (MQ=255)
                                                                cGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTATGCCGCACGTAAAACTGCGCGACGTAGTTATCGGTATGGGTGCATGCACCGAttac                                    >  2:23342/1‑98 (MQ=255)
                                                                  gCGTGAAGAAAATTATCCGCGTGGGTCCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCaa                                  <  1:759853/100‑1 (MQ=255)
                                                                    gTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATTGGTGCCTGCACCGATTCCAAAg                                <  2:1208854/100‑1 (MQ=255)
                                                                    gTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAg                                <  2:51892/100‑1 (MQ=255)
                                                                     tGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGt                               <  1:893803/100‑1 (MQ=255)
                                                                     tGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGt                               <  2:79283/100‑1 (MQ=255)
                                                                      gaagaaAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGtt                              >  2:739207/1‑100 (MQ=255)
                                                                             aTTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCa                       >  1:1590776/1‑100 (MQ=255)
                                                                                     cgTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGtttt               <  2:1488799/100‑1 (MQ=255)
                                                                                     cgTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGtttt               <  2:609714/100‑1 (MQ=255)
                                                                                      gggggTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTa              <  1:1241853/98‑1 (MQ=255)
                                                                                       tGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTaa             <  2:1590471/100‑1 (MQ=255)
                                                                                        gggTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTaaa            >  1:337517/1‑100 (MQ=255)
                                                                                           ttCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGAc         <  1:1978847/100‑1 (MQ=255)
                                                                                             ccTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCa       >  1:1125060/1‑100 (MQ=255)
                                                                                              cTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCa       >  1:297178/1‑99 (MQ=255)
                                                                                                gtgGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGa    <  1:94572/100‑1 (MQ=255)
                                                                                                gtgGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGa    <  1:94543/100‑1 (MQ=255)
                                                                                                gtgGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGa    <  2:736817/100‑1 (MQ=255)
                                                                                                gtgGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGa    >  2:867078/1‑100 (MQ=255)
                                                                                                  ggCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGAACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACt  <  2:778248/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACT  >  NC_000913/4621061‑4621258

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.