New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 442659134 (1.130)6 (0.230) 6/196 2.2 15.4% intergenic (+34/+37) rplI/yjfZ 50S ribosomal subunit protein L9/uncharacterized protein
?NC_000913 = 4426592 36 (1.360)intergenic (+35/+36) rplI/yjfZ 50S ribosomal subunit protein L9/uncharacterized protein
Rejected: Frequency below/above cutoff threshold.

CACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCG  >  NC_000913/4426479‑4426723
                                                                                                                 |                                                                                                                                   
cACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTa                                                                                                  >  2:139976/1‑149 (MQ=255)
 aCTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTAt                                                                                                 >  3:56632/1‑149 (MQ=255)
   tGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTTAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATtggtc                                                                                                                                >  6:171075‑M2/1‑111 (MQ=255)
     gCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGt                                                                                             >  6:45017/1‑149 (MQ=255)
        aaCACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTc                                                                                          <  1:139976/149‑1 (MQ=255)
        aaCACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTc                                                                                          <  5:45017/149‑1 (MQ=255)
                 tGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCtt                                                                                                                    <  4:135402/114‑1 (MQ=255)
                 tGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCtt                                                                                                                    >  3:135402/1‑114 (MQ=255)
                                 caGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATtctc                                                                 <  4:56632/149‑1 (MQ=255)
                                                      gATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCtt                                                               <  1:84785/130‑1 (MQ=255)
                                                      gATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCtt                                                               >  2:84785/1‑130 (MQ=255)
                                                      gATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATg                                            >  4:26418/1‑149 (MQ=255)
                                                       aTCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATtggtcggcgtttttacgtctcgttgaataacgaa                                                                                                   >  1:57564‑M2/1‑59 (MQ=255)
                                                       aTCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATtggtcggcgtttttacgtctcgttgaataacgaa                                                                                                   <  2:57564‑M2/93‑35 (MQ=255)
                                                        tCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATtctc                                                                 <  8:31074/126‑1 (MQ=255)
                                                        tCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGAAGATCCATTCACCATtctc                                                                 >  7:31074/1‑126 (MQ=255)
                                                          gTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCAt                                                                     >  6:119942/1‑120 (MQ=255)
                                                          gTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCAt                                                                     <  5:119942/120‑1 (MQ=255)
                                                           tAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATtggtcggcgtttttacgtctcgttgaataacgaattattcagctactacgtttacgatcactttcgcgaatacttcgctgtgaacctggaagct                                       <  2:16502‑M2/149‑95 (MQ=255)
                                                                    tagCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATtggtcggcgtttttacgtctcgttgaataacgaattattcagctactacgtttacgatcactttcgcgaatacttcgctgtgaacctggaagctcacttcgtg                              >  1:16502‑M2/1‑46 (MQ=255)
                                                                        tGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCa                                                        <  6:130589/119‑1 (MQ=255)
                                                                        tGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCa                                                        >  5:130589/1‑119 (MQ=255)
                                                                          aataatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGat                           <  4:11549/146‑1 (MQ=255)
                                                                          aataatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGat                           >  3:11549/1‑146 (MQ=255)
                                                                               ttCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAg                                                                                                     >  1:82641/1‑67 (MQ=255)
                                                                               ttCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAg                                                                                                     <  2:82641/67‑1 (MQ=255)
                                                                               ttCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATaa                                                     >  7:112235/1‑115 (MQ=255)
                                                                               ttCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATaa                                                     <  8:112235/115‑1 (MQ=255)
                                                                                         aaCGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGatacata                      >  2:167088/1‑136 (MQ=255)
                                                                                         aaCGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGatacata                      <  1:167088/136‑1 (MQ=255)
                                                                                               aCGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTg                                    >  5:151519/1‑116 (MQ=255)
                                                                                               aCGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTg                                    <  6:151519/116‑1 (MQ=255)
                                                                                                cGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCg  <  3:26418/149‑1 (MQ=255)
                                                                                                 gTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATtctc                                                                 >  8:28592/1‑85 (MQ=255)
                                                                                                 gTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATtctc                                                                 <  7:28592/85‑1 (MQ=255)
                                                                                                    aaaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATgg                                           >  4:79570/1‑104 (MQ=255)
                                                                                                    aaaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATgg                                           <  3:79570/104‑1 (MQ=255)
                                                                                                     aaaCGCCGACCATtggtcggcgtttttacgtctcgttgaataacgaattattcagctactacgtttacgatcactttcgcgaatacttcgctgtgaacctggaagctcacttcgtgttcgc                         >  5:89055‑M2/1‑13 (MQ=255)
                                                                                                     aaaCGCCGACCATtggtcggcgtttttacgtctcgttgaataacgaattattcagctactacgtttacgatcactttcgcgaatacttcgctgtgaacctggaagctcacttcgtgttcgc                         <  6:89055‑M2/121‑109 (MQ=255)
                                                                                                               cATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATc                                                                               <  4:149002/57‑1 (MQ=255)
                                                                                                               cATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATc                                                                               >  3:149002/1‑57 (MQ=255)
                                                                                                                 ttGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCAt                                                                            <  6:287/58‑1 (MQ=255)
                                                                                                                 ttGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCAt                                                                            >  5:287/1‑58 (MQ=255)
                                                                                                                 |                                                                                                                                   
CACTGGCGAACACGAAGTGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGGCGTTAGATACATACGAGCGCAACTGATGTCTCG  >  NC_000913/4426479‑4426723

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