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

TAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATT  >  NC_000913/279836‑279951
                |                                                                                                   
tAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCTACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGaa                  <  1:407587/100‑1 (MQ=17)
tAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGAAGCATAACAGAGaa                  >  1:2973059/1‑100 (MQ=255)
 aGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAac                 <  1:3043299/100‑1 (MQ=33)
 aGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAac                 >  2:10831/1‑100 (MQ=33)
  gCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAaca                >  2:1188542/1‑100 (MQ=35)
   cTGAACAGTAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAacac               >  2:911310/1‑100 (MQ=25)
   cTGAACAGGAGGGACAGCTGATAGAAACAGTAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTTGCAGCATCACCGAGAacac               <  2:1148537/100‑1 (MQ=255)
     gAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACaa             >  1:1106799/1‑100 (MQ=38)
     gAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACaa             >  2:15781/1‑100 (MQ=38)
     gAACAGGAGGGACAGCTGATAGAAACAGAAGCAACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACaa             >  2:1875511/1‑100 (MQ=255)
      aaCAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCTGTAAGTTGGCAGCATCACCGAGAACACaaa            <  2:963727/100‑1 (MQ=255)
      aaCAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACaaa            >  1:2221954/1‑100 (MQ=40)
       aCAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAAt           <  1:1188553/100‑1 (MQ=40)
        cAGGAGGGACAGCTGATAGAAACAGAAACCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACaaa            >  2:1158293/1‑98 (MQ=255)
         aggaggGACAGTTGATAGAAACAGAAGCCACTTGTGCACCTCAAAAACACCATCATACACTAAAGCAGGAAGTTGGCAGCATCACCGAGAACACAAATg          <  1:1398647/99‑1 (MQ=14)
           gaggGACAGCTAATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAATTCAGTAAGTTGGCAGCATCAAAGAGAACACAAATGGGc       <  2:2609246/100‑1 (MQ=14)
            aggGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCt      <  1:2488948/100‑1 (MQ=255)
              ggACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTaa    <  2:1083414/100‑1 (MQ=255)
                aCAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAAtt  <  1:1644595/100‑1 (MQ=255)
                |                                                                                                   
TAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAGAACACAAATGGGCTAATT  >  NC_000913/279836‑279951

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