New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 279841NA (NA)14 (0.160) 5/168 NT NA noncoding (90/768 nt) IS1 repeat region
?NC_000913 = 279852 NA (NA)noncoding (79/768 nt) IS1 repeat region

GTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAAT  >  NC_000913/279835‑279950
                 |                                                                                                  
gTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCgaga                  <  1:2347429/100‑1 (MQ=31)
  aGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAac                <  1:895657/100‑1 (MQ=33)
   gCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAaca               >  2:598839/1‑100 (MQ=35)
   gCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAaca               >  2:973095/1‑100 (MQ=35)
     tGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAAcaca             >  2:722490/1‑100 (MQ=38)
        aCAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAAt          <  1:1651542/100‑1 (MQ=40)
        aCAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAAt          >  1:438471/1‑100 (MQ=40)
         cAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTACGTTGGCAGCATCACCGAGAACACAAAt          >  2:1328866/1‑99 (MQ=255)
         cAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATg         <  1:1884695/100‑1 (MQ=255)
           ggaggGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATggg       >  1:1639979/1‑100 (MQ=255)
             aggGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCt     <  2:1326498/100‑1 (MQ=255)
                gACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAATTTGGCAGCATCACCGCGAACACAAATGGGCTAAt  >  2:773159/1‑100 (MQ=255)
                gACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAAt  >  1:340071/1‑100 (MQ=255)
                gACAGCTGATAGAAACAGAAGCCACAAAAAAAAAAAAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAAt  <  2:446181/100‑1 (MQ=17)
                 aCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAACTCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAAt  >  2:1665359/1‑99 (MQ=255)
                 |                                                                                                  
GTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAAT  >  NC_000913/279835‑279950

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.