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.

ATCAATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_012967/3947866‑3947791
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gttttAACGGCAATTAATATATATGTTATTAATTGAATGAATTTTATT  >  NC_012967/3948049‑3948091
                                                                                                                       
ATCAATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCG                                                 >  2:9008296/1‑72
 TCAATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGT                                                <  2:1677510/72‑1
  CAATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTT                                               <  1:2280055/72‑1
  CAATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTT                                               <  1:1124130/72‑1
  CAATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTT                                               <  2:7140981/72‑1
  CAATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTT                                               >  2:5354254/1‑72
  CAATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTT                                               <  1:7772508/72‑1
  CAATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTT                                               >  2:1294119/1‑72
   AATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCNGAGCATTTTGCATGGAATTTTCGTTT                                              >  1:674202/1‑72
   AATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTT                                              >  2:9527055/1‑72
   AATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTT                                              >  1:9757420/1‑72
    ATTATCAATGCTTANNNACGTCAACATTTTGACGNCATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTT                                             <  1:7029319/72‑1
    ATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTT                                             <  2:6220909/72‑1
    ATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTT                                             >  1:7372036/1‑72
    ATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTT                                             >  1:7670166/1‑72
    ATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTT                                             <  2:6220890/72‑1
    ATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTT                                             <  2:5744877/72‑1
    ATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTT                                             <  1:7745373/72‑1
    ATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTT                                             >  2:352230/1‑72
    ATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTT                                             <  1:8510369/72‑1
         CAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTAACGG                                        <  2:3399075/72‑1
           ATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTAACGGCA                                      <  2:6946082/72‑1
            TGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTAACGGCAA                                     <  2:4507269/72‑1
               TTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTAACGGCAATTA                                  <  2:3836881/72‑1
                        TCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTAACGGCAATTAATATATATG                         <  2:7534591/72‑1
                        TCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTAACGGCAATTAATATATATG                         <  2:7534609/72‑1
                          AACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTAACGGCAATTAATATATATGTT                       <  2:7378752/72‑1
                               TTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTTAACGGCAATTAATATATATGTTATTAA                  <  2:6295558/72‑1
                                     aACATGAAACCGGAGCATTTGGCATCGATTTTTCTTTTTATCGGCAATTAATATATATGTTATTAATTGAAT            <  2:5744642/71‑1
                                               cgaAGCATTTGGCATGAAATTTTCGTTTTCACGGCAATTAATATATATGTTATTAATTGAATGAATTTTATT  <  2:9298899/69‑1
                                                                                                                       
ATCAATTATCAATGCTTAATTACGTCAACATTTTGACGACATGAAACCGGAGCATTTTGCATGGAATTTTCGTTTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_012967/3947866‑3947791
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gttttAACGGCAATTAATATATATGTTATTAATTGAATGAATTTTATT  >  NC_012967/3948049‑3948091

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ 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.