New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 20474NA (NA)18 (0.150) 6/168 NT NA noncoding (90/768 nt) IS1 repeat region
?NC_000913 = 20485 NA (NA)noncoding (79/768 nt) IS1 repeat region

AGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATG  >  NC_000913/20467‑20573
       |                                                                                                   
aGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTACATCAGTAAGTTGGCAGCATCaccta          >  1:618670/1‑99 (MQ=11)
aGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCacctac         >  2:1747008/1‑100 (MQ=31)
 gTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCGTCATACACTAAATCAGTAAGTTGGCAGCATCacctacc        <  1:2083310/100‑1 (MQ=17)
 gTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCacctacc        <  1:1021387/100‑1 (MQ=32)
 gTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCacctacc        >  2:3205300/1‑100 (MQ=32)
 gTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGCAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCacctacc        <  2:135998/100‑1 (MQ=17)
  tggCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCacctacct       >  2:2280863/3‑100 (MQ=18)
   atcTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTc      <  1:129401/98‑1 (MQ=25)
    gCTGGACAGGAGGAACAGCTGATAGAAACAGAAGCAAATGGAGCAGCTAAAAAACACCATCATACACTAAATCAGTAAGTTGGAAGCATGACCTACCTCa     >  1:161561/1‑100 (MQ=11)
    gCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCa     >  2:2224654/1‑100 (MQ=35)
    gCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCa     <  2:2967717/100‑1 (MQ=35)
     cTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCaa    <  2:385061/100‑1 (MQ=38)
      tGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATCCACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAAt   <  1:826056/100‑1 (MQ=255)
      tGAACAGGAGCGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAAt   <  1:265321/100‑1 (MQ=255)
       gAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCagat  >  2:129401/1‑97 (MQ=25)
       |                                                                                                   
AGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTACCTCAATG  >  NC_000913/20467‑20573

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