New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 334914534 (0.820)4 (0.100) 4/276 5.7 11.2% coding (665/729 nt) mtgA biosynthetic peptidoglycan transglycosylase
?NC_000913 = 3349179 30 (0.740)coding (631/729 nt) mtgA biosynthetic peptidoglycan transglycosylase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCGCCAAATTCGGCG  >  NC_000913/3349037‑3349318
                                                                                                                                              |                                                                                                                                           
cAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCga                                                                                                                                             <  1:608120/143‑1 (MQ=255)
 aaCTTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgag                                                                                                                                            <  1:268467/143‑1 (MQ=255)
 aaCCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgag                                                                                                                                            <  1:648363/143‑1 (MQ=255)
        ggagaGCAAAACGCTGGCCCTCGCTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTGGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACtt                                                                                                                                     <  2:410520/142‑1 (MQ=255)
         gagaGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACttt                                                                                                                                    >  1:419634/1‑143 (MQ=255)
                 aaCGCTGGCCCTCGCTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGaa                                                                                                                            <  1:418011/143‑1 (MQ=255)
                 aaCGCTGGCCCTCACTTCGCGCGCAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGCCGGCTACACACGTACCGCGATGACGAGGAGACGTTCAAACGat                                                                                                                            >  2:79276/1‑142 (MQ=255)
                           cTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATACACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTCGCCCGATGGCGAGGCGACGTTGAAACGAAGAGGACTAgc                                                                                                                  >  2:75638/1‑142 (MQ=255)
                                ttCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGCGACTTTGAAACGAAGTGGATTAGGTAATa                                                                                                             >  1:125792/1‑143 (MQ=255)
                                           cTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGCCTTTGAAACGACGTGGATTAGGGAATAGAGCTGCCAGt                                                                                                  >  2:426972/1‑143 (MQ=255)
                                           cTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGt                                                                                                  >  1:271429/1‑143 (MQ=255)
                                           cTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGCGGATTAGGTAATACAGCTGCCAGt                                                                                                  >  1:96225/1‑143 (MQ=255)
                                               gtCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATg                                                                                              <  1:723584/142‑1 (MQ=255)
                                                   aGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGCTGCTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAgg                                                                                          >  1:45606/1‑142 (MQ=255)
                                                                   tAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTCGCCCGATGGCGACGCGACTTTGAAACGAAGTGGATTAGGTAATACATCTGCCAGTAATGCAGCTTCCGACCGGGaaaat                                                                          >  1:529917/1‑141 (MQ=255)
                                                                   tAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAGTACAGCTGCCAGTAATGCAGCTTCCGCCCGGGTAAGt                                                                          >  1:600919/1‑143 (MQ=255)
                                                                    aaaCGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGtt                                                                         >  1:741933/1‑143 (MQ=255)
                                                                        ggatCACCGCCTAACTCATACATCTGCCGTAAATTCCCCGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCt                                                                     <  2:33767/140‑1 (MQ=255)
                                                                             aCCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGCCTTTGAAACGAAGCGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCCcgcgc                                                                >  2:133756/1‑142 (MQ=255)
                                                                             aCCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCgg                                                                >  1:230421/1‑143 (MQ=255)
                                                                             aCCGCCTAACTGATACAGCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCgg                                                                <  2:709888/143‑1 (MQ=255)
                                                                                     aCTGATACATCTGCCGTAAAATACACGCATGACGGGTACGCACGTAGCACGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCGAGTAATGCAGCTTAGGAGCGGGTAAGTTTGCTAGCGGGAtggtgg                                                        >  1:257511/1‑143 (MQ=255)
                                                                                             aTCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCGGTAATGCAGCTACCGACAGGGGAAGTTAGCTCGCGGGGTTGTGGAaacaaca                                                >  2:514860/1‑142 (MQ=255)
                                                                                                  ccGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTtgtg                                           <  2:52328/143‑1 (MQ=255)
                                                                                                     tAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCATTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCg                                        >  2:373970/1‑143 (MQ=255)
                                                                                                       aagtgcgcGCCTGAGGGCTGCGCAGGTAGCGCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATGCAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCa                                      <  2:340565/136‑1 (MQ=255)
                                                                                                        aacccAGGCGTGGCGGGTACGCAGGTAGGCCGAGGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAg                                     <  1:305414/140‑1 (MQ=255)
                                                                                                        aaTCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTTTGGAACTACCGTTGTGCCGCAg                                     >  2:494191/1‑143 (MQ=255)
                                                                                                        aaTCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAg                                     >  1:252977/1‑143 (MQ=255)
                                                                                                                    cgggCTACGTGCGTAGCCCGATGGCGAGGAGACTTTGGACCGAAGTGGCGTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTATGTTTGCGCGCGGGATTGTGGAAATAACGTGGTGACGCAGCTTCGACGCCaa                         >  2:124080/3‑143 (MQ=255)
                                                                                                                      ggCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAac                       >  2:457926/1‑143 (MQ=255)
                                                                                                                                         gggCGAGGAGACCTTGAAACGAATGGGGTGAGGTAATGCGGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTGGCTCGCGGGTTGGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCGCCAAATTCgg    <  1:659699/142‑1 (MQ=255)
                                                                                                                                           gCGAGGGGACTTTGGAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCGCCAAATTCGGCg  <  2:208092/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
CAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCGCCAAATTCGGCG  >  NC_000913/3349037‑3349318

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