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

TACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCT  >  NC_000913/16634‑16751
                   |                                                                                                  
tACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATc                    >  2:436728/1‑100 (MQ=11)
tACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACACTTAAGTTAGCGATTATGGGATc                    >  2:709416/1‑100 (MQ=2)
 aCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCa                   <  2:611599/100‑1 (MQ=31)
 aCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCa                   >  1:1857189/1‑100 (MQ=31)
 aCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCa                   <  1:2747938/100‑1 (MQ=31)
  cAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCAc                  >  1:1645624/1‑100 (MQ=32)
   aaaaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACt                 >  2:2733374/1‑100 (MQ=33)
     aaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTcc               >  2:761824/1‑100 (MQ=35)
       cGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTcccc             >  2:3284611/1‑100 (MQ=38)
         cAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgc           >  1:3072306/1‑100 (MQ=40)
         cAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgc           <  2:3524136/100‑1 (MQ=40)
         cAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgc           >  1:2123859/1‑100 (MQ=40)
          aaTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgcc          >  1:1852725/1‑100 (MQ=40)
           aTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgccg         >  2:3463620/1‑100 (MQ=255)
           aTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgccg         >  1:2020062/1‑100 (MQ=255)
            ttGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGt        <  2:1143424/100‑1 (MQ=255)
            ttGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGt        >  2:2494725/1‑100 (MQ=255)
            ttGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGt        <  1:1691481/100‑1 (MQ=255)
             tGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGtt       <  1:715574/100‑1 (MQ=255)
              gaagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTg      >  2:2040490/1‑100 (MQ=255)
              gaagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTg      >  2:2863599/1‑100 (MQ=255)
              gaagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACACTTAAGTAAGCGATTATGGGATCACTCCCCGACGt        >  2:2272820/1‑98 (MQ=11)
                 gaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctc   >  2:2438665/1‑100 (MQ=255)
                 gaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctc   >  1:2642406/1‑100 (MQ=255)
                  aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctct  <  2:155472/100‑1 (MQ=255)
                  aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctct  <  1:709412/100‑1 (MQ=255)
                  aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctct  <  1:1429169/100‑1 (MQ=255)
                   aCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctct  <  2:3536926/99‑1 (MQ=255)
                   |                                                                                                  
TACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCT  >  NC_000913/16634‑16751

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