New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 17570850 (0.000)76 (1.960) 28/130 0.1 100% coding (1388/1413 nt) pykF pyruvate kinase I
?NC_000913 = 2290113 NA (NA)noncoding (1/1195 nt) IS5 repeat region

TGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCAT  >  NC_000913/2290048‑2290180
                                                                 |                                                                   
tgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCAt                                                                   >  2:968885/1‑68 (MQ=33)
tgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCAt                                                                   >  1:246301/1‑68 (MQ=33)
tgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCAt                                                                   >  1:243970/1‑68 (MQ=33)
tgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCAt                                                                   >  2:1287089/1‑68 (MQ=33)
   aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAAc                                                                >  1:816526/1‑68 (MQ=37)
    cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACg                                                               >  2:1227399/1‑68 (MQ=39)
     ccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGc                                                              <  2:243968/68‑1 (MQ=40)
          tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                         <  1:165852/68‑1 (MQ=255)
          tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                         >  2:262383/1‑68 (MQ=255)
          tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                         >  1:532583/1‑68 (MQ=255)
             aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGcca                                                      <  1:262385/68‑1 (MQ=255)
             aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCAATAACGCTGTAGcca                                                      >  1:623150/1‑68 (MQ=255)
               ggCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCGCCTTCCATAACGCTGTAGccacc                                                    >  2:1432558/1‑68 (MQ=255)
               ggCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGccacc                                                    >  1:1351825/1‑68 (MQ=255)
                 cTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAg                                                  <  1:1260744/68‑1 (MQ=255)
                  tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGa                                                 >  2:130531/1‑68 (MQ=255)
                  tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGa                                                 <  1:862387/68‑1 (MQ=255)
                  tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGa                                                 <  1:566695/68‑1 (MQ=255)
                  tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGa                                                 <  1:1287093/68‑1 (MQ=255)
                        tGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGat                                           <  2:816521/68‑1 (MQ=255)
                                aTGATCTCATATCAGTGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGa                                   <  1:390031/68‑1 (MQ=255)
                                    tctcATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGa                               >  1:1309419/1‑68 (MQ=255)
                                      tcATATCAGGGACTTGTTCGCACATTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACa                             >  1:519914/1‑68 (MQ=255)
                                      tcATATCAGGGACTAGATCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAAAGACa                             >  1:1033661/1‑68 (MQ=38)
                                          atCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAaga                         >  2:688592/1‑68 (MQ=255)
                                          atCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAaga                         <  2:649043/68‑1 (MQ=255)
                                          atCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAaga                         <  2:1310801/68‑1 (MQ=255)
                                          atCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGAAATTGCGGAACGACAAaga                         >  2:1185007/1‑68 (MQ=255)
                                           tCAGGGGCTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAagag                        >  2:1445576/1‑68 (MQ=255)
                                           tCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAagag                        >  2:381820/1‑68 (MQ=255)
                                           tCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAagag                        >  1:1111039/1‑68 (MQ=255)
                                           tCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAagag                        >  1:830544/1‑68 (MQ=255)
                                              gggACTTGTTCGCACCTTCCATAACGCTGTAGcca                                                      >  1:284258/1‑35 (MQ=255)
                                              gggACTTGTTCGCACCTTCCATAACGCTGTAGcca                                                      <  2:284256/35‑1 (MQ=255)
                                              gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGaaa                     >  2:149197/1‑68 (MQ=255)
                                              gggACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGaaa                     <  2:1275387/68‑1 (MQ=255)
                                                gACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACa                   <  2:1309415/68‑1 (MQ=255)
                                                gACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACa                   >  1:1023722/1‑68 (MQ=255)
                                                  cTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGa                 <  1:11208/68‑1 (MQ=255)
                                                  cTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGa                 <  1:1416051/68‑1 (MQ=255)
                                                  cTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGa                 >  2:620689/1‑68 (MQ=255)
                                                   ttGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGaa                <  2:246299/68‑1 (MQ=255)
                                                    tGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAAc               >  1:677876/1‑68 (MQ=255)
                                                       tCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAg            <  2:1023717/68‑1 (MQ=255)
                                                        cGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGa           <  1:1088968/68‑1 (MQ=255)
                                                        cGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGa           <  2:1131289/68‑1 (MQ=255)
                                                           aCCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTg        >  2:766685/1‑68 (MQ=255)
                                                            ccTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGa       >  1:1442376/1‑68 (MQ=255)
                                                               accATACCGCTGTAGCCACCAGATCTGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGc    <  2:549918/67‑1 (MQ=255)
                                                                 cATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCAt  <  1:227610/68‑1 (MQ=255)
                                                                 |                                                                   
TGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCAT  >  NC_000913/2290048‑2290180

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