New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 23256434 (0.760)6 (0.150) 5/262 5.2 16.0% intergenic (+15/+33) yafE/mltD putative S‑adenosyl‑L‑methionine‑dependent methyltransferase/putative membrane‑bound lytic murein transglycosylase D
?NC_000913_3_bme_pgi = 232571 32 (0.780)intergenic (+22/+26) yafE/mltD putative S‑adenosyl‑L‑methionine‑dependent methyltransferase/putative membrane‑bound lytic murein transglycosylase D
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TAGACGCTATTCGCATTTACCAACAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAAACAACGTCAGCTTATCGCCTGGTTGCAGATTCGCAGTGTCGCTGTTCCAGCGCATCACATCTTTGATGTTCACGC  >  NC_000913_3_bme_pgi/232428‑232703
                                                                                                                                        |                                                                                                                                           
tagACGCTATTCGCATTTACCAACAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTg                                                                                                                                 >  1:298067/1‑149 (MQ=255)
     gCTATTCGCATTTACCAACAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACATGGGGAATCGGTGCCttt                                                                                                                            >  2:308148/1‑149 (MQ=255)
       tATTCGCATTTACCAACAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCATCGGGGGAATCGGTGCCtttgt                                                                                                                          >  1:395934/1‑147 (MQ=255)
         ttCGCATTTACCAACAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTttat                                                                                                                        >  2:181343/1‑149 (MQ=255)
          tCGCATTTACCAACAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTttatt                                                                                                                       <  1:41186/149‑1 (MQ=255)
                    cAACAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTg                                                                                                             >  2:680653/1‑149 (MQ=255)
                     aaCAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGt                                                                                                            >  1:307090/1‑149 (MQ=255)
                            cgcATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAAt                                                                                                     >  1:622593/1‑149 (MQ=255)
                                      gagGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCAtgtt                                                                                           >  1:96558/1‑149 (MQ=255)
                                         gTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCgggg                                                                                                                                          <  1:659488/99‑1 (MQ=255)
                                         gTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCgggg                                                                                                                                          >  2:659488/1‑99 (MQ=255)
                                             gaactgaTTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCGGGTTTGTCAGGAATCTGGCATGTTGTTGttt                                                                                    <  1:64497/143‑1 (MQ=255)
                                              aaCGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGttttt                                                                                  >  2:421634/1‑149 (MQ=255)
                                                      ttGCCTCGCAGAATGATGGCTTTTTCCCCAGTGATATCATCATGGTAGATGCATATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCacaaac                                                                           <  2:333603/149‑1 (MQ=255)
                                                                   tgatgGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTca                                                                                >  2:532273/1‑131 (MQ=255)
                                                                   tgatgGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTca                                                                                <  1:532273/131‑1 (MQ=255)
                                                                                aCCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAGACAACGTCAGCTTATCGCCTGGTTGCAg                                                 >  1:740984/1‑149 (MQ=255)
                                                                                     tGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCa                                                                                               <  2:614769/98‑1 (MQ=255)
                                                                                     tGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCa                                                                                               >  1:614769/1‑98 (MQ=255)
                                                                                       atatCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAAACAACGTCAGCTTATCGCCTGGTTGCAGATTCGCa                                          <  2:96558/149‑1 (MQ=255)
                                                                                       atatCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGGGCCTTTTTATTATCTTGTTTGTCAGGCATCTGGCATGTTGTTGTTTTTCACAAACAACGTCAGCTTATCGCCTGGTTGCAGATTCGCa                                          <  2:515072/149‑1 (MQ=255)
                                                                                         atcatcATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAAACAACGTCAGCTTATCGCCTGGTTGCAGATTCGCAgt                                        <  2:307090/149‑1 (MQ=255)
                                                                                             tcatGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAAACAACATCAGCTTATCGCAAGATTGCAGGTTCGCAgtgtgc                                    >  2:204806/1‑147 (MQ=255)
                                                                                                 ggTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTACCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAAACAACGTCAGCTTATCGCCTGGTTGCAGATTCGCAGTCTCGCGGt                                >  2:562116/1‑149 (MQ=255)
                                                                                                  gTAGAGGGACAGAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAAACAACGTCAGCTTATCGCCTGGTTGCAGATTCGCAGTGTCGCTGtt                               <  1:226824/149‑1 (MQ=255)
                                                                                                   tAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAAACAACGTCAACTTATACACTGGTTGCAGATTCGCAGTGTCGCCGTTc                              >  2:39024/1‑149 (MQ=255)
                                                                                                      aTGCACATAAAGCGGCATAAATCAAAAAGGCACCGGGGGAATCGGTGCCTTTTGATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAAACAACGTCAGCTTATCGCCTGGTTGCAGATTCTCAGTGTCGCTGTTCCac                           >  2:365362/1‑148 (MQ=255)
                                                                                                            aTAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCAtgtt                                                                                           >  1:772742/1‑79 (MQ=255)
                                                                                                            aTAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCAtgtt                                                                                           <  2:772742/79‑1 (MQ=255)
                                                                                                                       taaataaaAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAAACAACGTCAGCTTATCGCCTGGTTGCAGATTCGCAGTGTCGCTGTTCCAGCGCATCACATCTTTGAt          >  1:644140/1‑149 (MQ=255)
                                                                                                                          ataaaAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAAACAACGTCAGCTTATCGCCTGGTTGCAGATGCGCAGTGTCGCTGTTCCAGCGCATCACATCTTTGATGtt       >  2:631596/1‑149 (MQ=255)
                                                                                                                           taaaAAAGGCACCGGGGGGATCGGTGCCTTTTTAGTATCTGGTTTGTCAGGAATCTGGTATGTTGTTGTTTTTCACAAACAACGTCAGCTTATCGCCTGGTTGTAGATTCGCAGTGTCGCTGTTGCAGCGCATCACATCTTTGATGtta      >  2:610092/1‑148 (MQ=255)
                                                                                                                             aaaaaGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAAACAACGTCAGCTTATCGCCTGGTTGCAGATTCGCAGTGTCGCTGTTCCAGCGCATCACATCTTTGATGTTCAc    >  2:10032/1‑149 (MQ=255)
                                                                                                                               aaaGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAAACAACGTCAGCTTATCGCCTGGTTGCAGATTCGCAGTGTCGCTGTTCCAGCCCATCACATCTTTGATGTTTACga  >  1:432756/1‑148 (MQ=255)
                                                                                                                                        |                                                                                                                                           
TAGACGCTATTCGCATTTACCAACAGAGCGCATCGACAGAGGTGAGAACGTATTTTGCCTTGCAGAATGATGGCTTTTTCACCAGTGATATCATCATGGTAGATGCACATAAAGCGGCATAAATAAAAAAGGCACCGGGGGAATCGGTGCCTTTTTATTATCTGGTTTGTCAGGAATCTGGCATGTTGTTGTTTTTCACAAACAACGTCAGCTTATCGCCTGGTTGCAGATTCGCAGTGTCGCTGTTCCAGCGCATCACATCTTTGATGTTCACGC  >  NC_000913_3_bme_pgi/232428‑232703

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.