New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 16653NA (NA)30 (0.210) 15/164 NT NA noncoding (1267/1345 nt) IS186 repeat region
?NC_000913 = 16691 NA (NA)noncoding (1305/1345 nt) IS186 repeat region

TACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTT  >  NC_000913/16634‑16752
                   |                                                                                                   
tACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATc                     >  1:1867556/1‑100 (MQ=11)
tACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATc                     <  2:3272406/100‑1 (MQ=11)
 aCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCa                    <  2:2836707/100‑1 (MQ=31)
 aCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCa                    >  2:277090/1‑100 (MQ=31)
 aCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCa                    >  1:3077522/1‑100 (MQ=31)
 aCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCa                    <  2:3837013/100‑1 (MQ=31)
  cAAAACGCAATTGAAGATCGCGAGGCATCGTCTTAACGAGGCACCGATGTGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCAc                   <  1:800163/100‑1 (MQ=255)
   aaaaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACt                  >  2:1309998/1‑100 (MQ=33)
    aaaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTc                 >  1:649714/1‑100 (MQ=35)
      aCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTccc               <  2:1099194/100‑1 (MQ=35)
        gCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccg             <  2:2777876/100‑1 (MQ=38)
           aTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgccg          >  2:600996/1‑100 (MQ=255)
            ttGAAGAACGCGAAGCATCGTCTTAACGAGGCACCGAGACGTCGCATTCTTCAGATGGTTAAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGt         >  1:2642074/1‑100 (MQ=255)
             tGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGtt        <  2:77950/100‑1 (MQ=255)
              gaagaaCTCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATCCTTCAGCTGGTTCAACCCTTAGGTTAGCGCTTATGGGATCACTccc               <  1:2178611/92‑1 (MQ=11)
              gaagaaCGCGAGGCATCGTCTTAAGGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTg       <  1:1972289/100‑1 (MQ=255)
               aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGc      <  1:2176973/100‑1 (MQ=255)
               aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGc      <  2:1453672/100‑1 (MQ=255)
               aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGc      <  2:3261401/100‑1 (MQ=255)
               aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCGTAAGTTAGCGCTTACGGGATCACTCCCCGCCGTTGc      <  1:3125230/100‑1 (MQ=255)
                agaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTGGCGCTTATGGGATCACTCCCCGCCGTTGct     <  1:2803989/100‑1 (MQ=255)
                agaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGct     <  2:549906/100‑1 (MQ=255)
                agaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGct     <  2:2062394/100‑1 (MQ=255)
                agaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGct     >  2:2744488/1‑100 (MQ=255)
                agaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGct     <  2:60608/100‑1 (MQ=255)
                 gaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctc    <  2:939640/100‑1 (MQ=255)
                 gaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctc    >  1:639104/1‑100 (MQ=255)
                 gaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctc    <  1:2050147/100‑1 (MQ=255)
                  aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctct   <  1:4244431/100‑1 (MQ=255)
                  aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctct   >  1:1502392/1‑100 (MQ=255)
                  aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctct   <  1:1296204/100‑1 (MQ=255)
                   aCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCtt  >  1:2396104/1‑100 (MQ=255)
                   aCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCtt  >  1:1400015/1‑100 (MQ=255)
                   |                                                                                                   
TACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTT  >  NC_000913/16634‑16752

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