New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_012967 3947791 =251 (0.880)10 (0.040) 8/132 16.9 6.6% coding (16/1656 nt) aslA arylsulfatase AslA
?NC_012967 3948049 = 50 (0.190)intergenic (‑243/+443) aslA/hemY arylsulfatase AslA/protoheme IX biogenesis protein HemY
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

Only 100 of 260 total aligned reads displayed.

TGACGGGGTATCTGCTGCGTTGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCATTGATAAT  >  NC_012967/3947720‑3947862
                                                                       |                                                                       
ggACGTGGTATCTGCTGCGTGGGCGATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGa                                                                         <  2:5236803/71‑1 (MQ=255)
 gagggggTAGCTGCGGCGTGGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGaa                                                                        <  2:6526943/69‑1 (MQ=255)
 gACTGGGTATCTGCTGCGTTGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGaa                                                                        <  1:65624/72‑1 (MQ=255)
  aCGGGGTATCTGCTGCGTTGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGaaa                                                                       <  1:9750082/72‑1 (MQ=255)
  aCGGGGTATCTGCTGCGTTGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGGTTGGGGGaaa                                                                       >  2:4612684/1‑72 (MQ=255)
   cGGGGTATCTGCTGCGTTGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGaaaa                                                                      >  1:2087288/1‑72 (MQ=255)
    ggggTATCTGCTGCGTTGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAAc                                                                     <  1:5381788/72‑1 (MQ=255)
           ctgctgCGTTGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAAtt                                                              <  1:5062169/72‑1 (MQ=255)
           ctgctgCGTTGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAAtt                                                              <  2:7772713/72‑1 (MQ=255)
            tgctgcGTTGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTc                                                             <  1:4298143/72‑1 (MQ=255)
            tgctgcGTTGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTc                                                             >  1:8766102/1‑72 (MQ=255)
             gctgcGTGGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTcc                                                            <  2:508906/72‑1 (MQ=255)
              ctgcGTTGGCCATGAGGGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCa                                                           <  1:5741135/72‑1 (MQ=255)
               tgcGTTGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCAt                                                          <  1:9564667/72‑1 (MQ=255)
               tgcGTTGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCAt                                                          <  1:8603678/72‑1 (MQ=255)
                   ttGGCCATGAGTGGAAGAGCTGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCaa                                                      >  1:7342045/1‑72 (MQ=255)
                   tgggCCGTGGGTGGAAGAGCGGGGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCaa                                                      <  2:8024072/70‑1 (MQ=255)
                    tGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCaaa                                                     <  1:8130075/72‑1 (MQ=255)
                       ccATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCaaaaag                                                  >  2:308468/1‑70 (MQ=255)
                       ccATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATg                                                  <  1:3669758/72‑1 (MQ=255)
                       ccATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATg                                                  <  1:2446634/72‑1 (MQ=255)
                       ccATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATg                                                  >  2:6030855/1‑72 (MQ=255)
                           gAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTcc                                              <  1:8472524/72‑1 (MQ=255)
                           gAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTcc                                              <  2:6891982/72‑1 (MQ=255)
                            aGTGGAAGAGCGGTGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCg                                             >  2:2583603/1‑72 (MQ=255)
                             gTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCgg                                            <  1:3885677/72‑1 (MQ=255)
                              tGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGt                                           <  1:4192212/72‑1 (MQ=255)
                              tGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAAATCCATGGAAAAAGCGcccgc                                           >  2:3661338/1‑69 (MQ=255)
                               ggAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGtt                                          <  2:1846033/72‑1 (MQ=255)
                               ggAAGAGCGGCGGCGACAGCTACAACAAGAACTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGtt                                          <  1:142786/72‑1 (MQ=255)
                                gAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGttt                                         <  2:9107090/72‑1 (MQ=255)
                                gAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGttt                                         <  1:7463995/72‑1 (MQ=255)
                                 aaGGGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTc                                        <  2:2495961/72‑1 (MQ=255)
                                  agagCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCa                                       <  2:9633426/72‑1 (MQ=255)
                                  agagCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCa                                       >  1:801442/1‑72 (MQ=255)
                                   gagCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCAt                                      <  2:8673595/72‑1 (MQ=255)
                                   gagCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCAt                                      <  1:1917574/72‑1 (MQ=255)
                                    agCGGCGGCGACAGCTACAACATGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATg                                     <  1:2716076/72‑1 (MQ=255)
                                    agCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATg                                     >  2:8819131/1‑72 (MQ=255)
                                    agCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATg                                     <  1:1482832/72‑1 (MQ=255)
                                     gcggcggcgACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATgt                                    <  1:3194375/72‑1 (MQ=255)
                                       ggcggcgACAGCTACAACNAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATgtcg                                  <  1:3824132/72‑1 (MQ=255)
                                       ggcggcgACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATgtcg                                  >  2:8451677/1‑72 (MQ=255)
                                          ggcgACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCa                               <  2:5409914/72‑1 (MQ=255)
                                          ggcgACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCa                               >  1:1039804/1‑72 (MQ=255)
                                           gcgACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCaa                              >  1:7834527/1‑72 (MQ=255)
                                            cgACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCaaa                             <  1:6513130/72‑1 (MQ=255)
                                             gACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCaaaa                            <  1:5199133/72‑1 (MQ=255)
                                             gACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCaaaa                            >  1:5536671/1‑72 (MQ=255)
                                              aCAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAAt                           >  1:7253073/1‑72 (MQ=255)
                                               cAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATg                          >  2:7152489/1‑72 (MQ=255)
                                                accTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGt                         >  1:3372101/3‑72 (MQ=255)
                                                 gCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGtt                        >  1:4786943/1‑72 (MQ=255)
                                                  cTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTg                       >  1:2773440/1‑72 (MQ=255)
                                                  cTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTg                       >  1:9250102/1‑72 (MQ=255)
                                                  cTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTg                       <  1:8891460/72‑1 (MQ=255)
                                                   tACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGa                      <  2:7088435/72‑1 (MQ=255)
                                                    acaacaAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGAc                     <  2:1969548/72‑1 (MQ=255)
                                                    acaacaAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGAc                     >  2:3441023/1‑72 (MQ=255)
                                                     caacaaGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACg                    >  2:1127984/1‑72 (MQ=255)
                                                     caacaaGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGNTGACg                    <  1:948971/72‑1 (MQ=255)
                                                      aacaaGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGt                   <  1:1978587/72‑1 (MQ=255)
                                                      aacaaGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGt                   <  1:4740563/72‑1 (MQ=255)
                                                      aacaaGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGt                   >  2:7254752/1‑72 (MQ=255)
                                                       acaaGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTa                  <  1:8021944/72‑1 (MQ=255)
                                                       acaaGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTa                  >  1:4410399/1‑72 (MQ=255)
                                                       acaaGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTa                  >  2:9388311/1‑72 (MQ=255)
                                                       acaaGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTa                  <  1:6487189/72‑1 (MQ=255)
                                                        caaGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTaa                 <  2:4198879/72‑1 (MQ=255)
                                                         aagaagTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAAt                <  1:5850437/72‑1 (MQ=255)
                                                         aagaagTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAAt                >  2:6992333/1‑72 (MQ=255)
                                                         aagaagTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAAt                >  2:4991942/1‑72 (MQ=255)
                                                          agaagTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAAtt               <  1:8969633/72‑1 (MQ=255)
                                                          agaagTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAAtt               <  1:9634538/72‑1 (MQ=255)
                                                            aagTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTaa             <  1:2505869/72‑1 (MQ=255)
                                                             agTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAg            <  2:2288354/72‑1 (MQ=255)
                                                              gTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGc           <  2:3278199/72‑1 (MQ=255)
                                                              gTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGc           <  2:2728636/72‑1 (MQ=255)
                                                              gTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGc           <  2:7611236/72‑1 (MQ=255)
                                                               tttGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCa          >  1:1720354/1‑72 (MQ=255)
                                                               tttGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCa          >  1:6870836/1‑72 (MQ=255)
                                                                ttGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCat         <  1:8243748/72‑1 (MQ=255)
                                                                ttGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCat         >  2:6448042/1‑72 (MQ=255)
                                                                 tGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCatt        <  1:3066308/72‑1 (MQ=255)
                                                                 tGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCatt        <  2:1498637/72‑1 (MQ=255)
                                                                  gggTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCattg       <  1:3653818/72‑1 (MQ=255)
                                                                  gggTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCattg       <  2:4068740/72‑1 (MQ=255)
                                                                   ggTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCattga      <  2:8273335/72‑1 (MQ=255)
                                                                    gTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCattgat     <  1:1119561/72‑1 (MQ=255)
                                                                    gTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCattgat     >  2:1450948/1‑72 (MQ=255)
                                                                    gTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCattgat     <  1:7652750/72‑1 (MQ=255)
                                                                     tGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCattgata    >  2:4004465/1‑72 (MQ=255)
                                                                     tGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCattgata    <  1:9948907/72‑1 (MQ=255)
                                                                     tGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCattgata    <  1:9388413/72‑1 (MQ=255)
                                                                     tGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTAACGTAATTAAGCattgata    <  2:1465618/72‑1 (MQ=255)
                                                                      gAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCattgataa   <  2:7724355/72‑1 (MQ=255)
                                                                       aaaaCGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCattgataat  >  2:6220890/1‑72 (MQ=255)
                                                                       aaaaCGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCattgataat  >  2:5744877/1‑72 (MQ=255)
                                                                       aaaaCGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCattgataat  <  2:352230/72‑1 (MQ=255)
                                                                       aaaaCGAAAATTCCATGCAAAATGCTCCGGTTTCATGNCGTCAAAATGTTGACGTNNNTAAGCattgataat  >  1:7029319/1‑72 (MQ=255)
                                                                       |                                                                       
TGACGGGGTATCTGCTGCGTTGGCCATGAGTGGAAGAGCGGCGGCGACAGCTACAACAAGAAGTTTGGGTGAAAACGAAAATTCCATGCAAAATGCTCCGGTTTCATGTCGTCAAAATGTTGACGTAATTAAGCATTGATAAT  >  NC_012967/3947720‑3947862

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ 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.