New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 16653NA (NA)17 (0.130) 10/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                     <  2:2260750/100‑1 (MQ=11)
tACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATc                     <  1:1428439/100‑1 (MQ=11)
tACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATc                     >  2:3296493/1‑100 (MQ=11)
tACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATc                     <  2:842901/100‑1 (MQ=11)
tACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATc                     >  2:1216843/1‑100 (MQ=11)
 aCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCa                    >  2:2319823/1‑100 (MQ=31)
  cAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCAc                   <  1:3833547/100‑1 (MQ=32)
  cAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCAc                   >  1:910965/1‑100 (MQ=32)
  cAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCAc                   >  1:45462/1‑100 (MQ=32)
  cAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCAc                   <  2:3675559/100‑1 (MQ=32)
   aaaaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACt                  >  1:3594637/1‑100 (MQ=33)
     aaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCCTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTcc                >  2:799931/1‑100 (MQ=25)
     aaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTcc                >  2:1461898/1‑100 (MQ=35)
     aaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGGTCAACCCTTAAGTTAGCGCTTATGGGATCACTcc                >  1:2366761/1‑100 (MQ=25)
      aCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTccc               <  2:3877381/100‑1 (MQ=35)
      aCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTccc               <  1:499163/100‑1 (MQ=38)
      aCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTccc               >  1:854455/1‑100 (MQ=38)
       ctcAATTGAACAACTCGAGGCATCGTCTTAACGTTGCATCGAGGCTTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATc                     <  2:3933163/91‑1 (MQ=11)
       cGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTcccc              >  2:3533790/1‑100 (MQ=38)
       cGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTcccc              >  2:3284542/1‑100 (MQ=35)
        gCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccg             <  1:3803438/100‑1 (MQ=38)
         cAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgc            >  2:2647392/1‑100 (MQ=40)
         cAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgc            >  1:1691672/1‑100 (MQ=38)
         cAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgc            <  1:1001998/100‑1 (MQ=40)
          aaTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTCTGGGATCACTCCccgcc           >  2:434665/1‑100 (MQ=255)
           aTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgccg          >  1:1167344/1‑100 (MQ=40)
             tGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATTACTcccataagcg        <  2:1261075/100‑8 (MQ=11)
             tGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGtt        >  1:34243/1‑100 (MQ=255)
             tGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGtt        >  1:3117074/1‑100 (MQ=255)
             tGAAGAACGCGAGGCATCGTCTTAACGAGGAACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGtt        >  1:691554/1‑100 (MQ=255)
              gaagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTg       >  2:3949009/1‑100 (MQ=255)
               aagcaCGCGAGGCCTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCCGATGGTTAAACCATTAAGTTAGCGCTTATGGGATCCTTCCTCGCCGTTGc      <  2:1310954/100‑1 (MQ=11)
               aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGc      <  2:2239702/100‑1 (MQ=255)
               aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGc      >  1:341524/1‑100 (MQ=255)
                agaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGct     >  1:2637537/1‑100 (MQ=255)
                 gaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctc    <  2:3258080/100‑1 (MQ=255)
                 gaaCGCGAGGCAGCGTCTTAACGAGGCACCGAGGCGTCGCATGCTTCAGATGGTGCAACCCTTAATTTAGCGCTTATGGGATCACTCCCCGCCGTTGctc    >  2:2148485/1‑100 (MQ=17)
                  aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctct   <  1:2539956/100‑1 (MQ=255)
                  aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctct   <  1:2517469/100‑1 (MQ=255)
                   aCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCtt  >  2:2756320/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 < 27 ≤ 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.