New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 291401NA (NA)137 (2.700) 43/128 0.0 70.3% noncoding (1/768 nt) IS1 repeat region
?NC_000913 873161 = 58 (1.140)coding (183/2349 nt) pdeI putative c‑di‑GMP phosphodiesterase PdeI

AGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGAT  >  NC_000913/291341‑291463
                                                            |                                                              
aGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAAt                                                            >  1:668068/1‑65 (MQ=33)
 agcACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAAtt                                                           <  1:3690871/63‑1 (MQ=25)
     cTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCa                                                       >  1:3574734/1‑65 (MQ=255)
          gCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCg                                                  >  1:2625684/1‑65 (MQ=255)
              cTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATg                                              >  1:34523/1‑65 (MQ=255)
                cAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGtt                                            <  1:3377252/65‑1 (MQ=255)
                 aaaaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGttt                                           <  1:445550/65‑1 (MQ=255)
                  aaaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTc                                          <  1:2721600/65‑1 (MQ=255)
                  aaaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTc                                          <  1:2408462/65‑1 (MQ=255)
                    aaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCAc                                        <  1:3065376/65‑1 (MQ=255)
                     acacCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACg                                       >  1:1963493/1‑65 (MQ=255)
                       acCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCActgt                                     >  1:1044185/1‑62 (MQ=255)
                       acCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCAcctg                                     >  1:2724391/1‑62 (MQ=255)
                       acCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTa                                     <  1:1882236/65‑1 (MQ=255)
                       acCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTa                                     <  1:2323028/65‑1 (MQ=255)
                       acCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTa                                     <  1:4033429/65‑1 (MQ=255)
                       acCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTa                                     <  1:439056/65‑1 (MQ=255)
                         catcatACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTAcc                                   <  1:3660733/65‑1 (MQ=255)
                          atcatACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCt                                  >  1:3338622/1‑65 (MQ=255)
                          atcatACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCt                                  >  1:2475237/1‑65 (MQ=255)
                           tcatACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTa                                 >  1:1737304/1‑65 (MQ=255)
                             atACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATg                               >  1:2743384/1‑65 (MQ=255)
                             atACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATg                               >  1:1427141/1‑65 (MQ=255)
                               acacTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTc                             <  1:2976298/65‑1 (MQ=255)
                                 acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTGCCTATGTCGa                           >  1:2597601/1‑65 (MQ=255)
                                 acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACTTATGTCgc                           >  1:2495103/1‑64 (MQ=255)
                                 acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGa                           >  1:2887438/1‑65 (MQ=255)
                                 acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGa                           >  1:1111397/1‑65 (MQ=255)
                                 acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGa                           >  1:3562153/1‑65 (MQ=255)
                                 acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGa                           >  1:3528676/1‑65 (MQ=255)
                                 acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGa                           >  1:2865014/1‑65 (MQ=255)
                                 acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGa                           >  1:2490554/1‑65 (MQ=255)
                                 acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGa                           >  1:1748818/1‑65 (MQ=255)
                                 acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGa                           >  1:1247499/1‑65 (MQ=255)
                                 acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGa                           >  1:2016009/1‑65 (MQ=255)
                                       tCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCt                     >  1:3396813/1‑65 (MQ=255)
                                       tCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCt                     >  1:1425251/1‑65 (MQ=255)
                                       tCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCt                     >  1:1619077/1‑65 (MQ=255)
                                       tCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCt                     >  1:122332/1‑65 (MQ=255)
                                       tCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCt                     >  1:2508890/1‑65 (MQ=255)
                                         aGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTgg                   >  1:242087/1‑65 (MQ=255)
                                         aGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTgg                   >  1:2833213/1‑65 (MQ=255)
                                          gTAGGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTGCCTATGTCGAATATCTGGc                  <  1:3695392/65‑1 (MQ=255)
                                          gTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGc                  <  1:1181295/65‑1 (MQ=255)
                                           tAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGcc                 <  1:1833571/65‑1 (MQ=255)
                                             aGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGa               >  1:2359964/1‑65 (MQ=255)
                                                 ggCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAg           <  1:3310559/65‑1 (MQ=255)
                                                 ggCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAg           <  1:389415/65‑1 (MQ=255)
                                                   cagcaTCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGt         <  1:2871676/65‑1 (MQ=255)
                                                   cagcaTCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGt         <  1:3596918/65‑1 (MQ=255)
                                                   cagcaTCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGt         <  1:3098924/65‑1 (MQ=255)
                                                     gcaTCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATACGGGTg       >  1:913621/1‑65 (MQ=255)
                                                     gcaTCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTg       >  1:1360358/1‑65 (MQ=255)
                                                     gcaTCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTg       >  1:434432/1‑65 (MQ=255)
                                                     gcaTCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTg       >  1:3554674/1‑65 (MQ=255)
                                                     gcaTCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTg       >  1:655640/1‑65 (MQ=255)
                                                          aCCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGAt  >  1:1891308/1‑65 (MQ=255)
                                                            |                                                              
AGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGAT  >  NC_000913/291341‑291463

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