New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 291401NA (NA)57 (1.270) 27/134 0.1 53.3% noncoding (1/768 nt) IS1 repeat region
?NC_000913 873161 = 50 (1.110)coding (183/2349 nt) pdeI putative c‑di‑GMP phosphodiesterase PdeI

AGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGA  >  NC_000913/291338‑291466
                                                               |                                                                 
aGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAAt                                                               <  1:50173/68‑1 (MQ=33)
   aGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAAtttt                                                            >  1:3213618/1‑68 (MQ=39)
     ccACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCa                                                          >  1:2978688/1‑68 (MQ=255)
            aGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGa                                                   >  1:1003309/1‑68 (MQ=255)
              cACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATg                                                 <  1:580197/68‑1 (MQ=255)
               aCCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGt                                                >  1:2277637/1‑68 (MQ=255)
               aCCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGt                                                >  1:1189402/1‑68 (MQ=255)
                 cTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGttt                                              <  1:3227086/68‑1 (MQ=255)
                 cTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGttt                                              >  1:2688230/1‑68 (MQ=255)
                  tCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTc                                             <  1:2599085/68‑1 (MQ=255)
                    aaaaaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCAc                                           <  1:404059/68‑1 (MQ=255)
                       aaCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGCTTCACGTa                                        <  1:2206332/68‑1 (MQ=255)
                         cacCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTAcc                                      <  1:2591736/68‑1 (MQ=255)
                          acCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCt                                     >  1:3392951/1‑68 (MQ=255)
                          acCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCt                                     >  1:3030904/1‑68 (MQ=255)
                            catcatACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTAt                                   <  1:303883/68‑1 (MQ=255)
                              tcatACACTAAATCAGTAAGTTGGCANCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGt                                 <  1:2860882/68‑1 (MQ=255)
                              tcatACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGt                                 <  1:787634/68‑1 (MQ=255)
                              tcatACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGt                                 <  1:584982/68‑1 (MQ=255)
                               catACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTc                                >  1:2907225/1‑68 (MQ=255)
                               catACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTc                                >  1:1832194/1‑68 (MQ=255)
                                  acacTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGaa                             >  1:2590960/1‑68 (MQ=255)
                                    acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAata                           >  1:57952/1‑68 (MQ=255)
                                    acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAata                           >  1:1847781/1‑68 (MQ=255)
                                    acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAata                           >  1:2695765/1‑68 (MQ=255)
                                    acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAata                           >  1:780220/1‑68 (MQ=255)
                                    acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAata                           >  1:2227602/1‑68 (MQ=255)
                                    acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAata                           >  1:1312741/1‑68 (MQ=255)
                                    acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAata                           >  1:1677307/1‑68 (MQ=255)
                                    acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAata                           >  1:934991/1‑68 (MQ=255)
                                    acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGACGAata                           >  1:1361662/1‑68 (MQ=255)
                                    acTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCAGATGTTTCACGTACCTATGTCGAata                           >  1:2875884/1‑68 (MQ=255)
                                     cTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAatat                          >  1:508785/1‑68 (MQ=255)
                                          tCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGc                     >  1:1715749/1‑68 (MQ=255)
                                          tCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGc                     >  1:1150495/1‑68 (MQ=255)
                                          tCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGc                     >  1:3217557/1‑68 (MQ=255)
                                          tCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGc                     >  1:764028/1‑68 (MQ=255)
                                          tCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGc                     >  1:2918130/1‑68 (MQ=255)
                                            aGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCg                   >  1:2404138/1‑68 (MQ=255)
                                              tAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGAt                 >  1:1442934/1‑68 (MQ=255)
                                                aGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATaa               >  1:1062278/1‑68 (MQ=255)
                                                 gTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAg              <  1:632966/68‑1 (MQ=255)
                                                  ttGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAgg             <  1:74109/68‑1 (MQ=255)
                                                   tGGCAGCATCGCCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGt            <  1:2120721/68‑1 (MQ=255)
                                                   tGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGt            <  1:2224914/68‑1 (MQ=255)
                                                   tGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGt            <  1:1310459/68‑1 (MQ=255)
                                                        gcaTCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCgg       >  1:936973/1‑68 (MQ=255)
                                                        gcaTCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCgg       >  1:306120/1‑68 (MQ=255)
                                                        gcaTCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCgg       >  1:2749075/1‑68 (MQ=255)
                                                        gcaTCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCgg       >  1:1258084/1‑68 (MQ=255)
                                                        gcaTCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTCGCgg       >  1:2428545/1‑68 (MQ=255)
                                                             aCCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCga  >  1:1448186/1‑68 (MQ=255)
                                                               |                                                                 
AGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGAATTTTCAGCCCGGATGTTTCACGTACCTATGTCGAATATCTGGCCGATAAGGTTGCGGATCGA  >  NC_000913/291338‑291466

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