New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 2765252 =26 (0.710)8 (0.230) 7/260 3.1 23.9% noncoding (15/199 nt) PP_mr33 Cobalamin regulatory region
?NC_002947 2765283 = 26 (0.740)intergenic (‑17/+15) PP_mr33/PP_2419 Cobalamin regulatory region/membrane protein
Rejected: Coverage evenness skew score above cutoff.

CGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACACCGCTGCAATACC  >  NC_002947/2765134‑2765417
                                                                                                                                                     |                                                                                                                                      
cGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCGCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCa                                                                                                                                        <  2:6076/150‑1 (MQ=255)
      gcAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGaa                                                                                                                                  >  2:813765/1‑150 (MQ=255)
                     tCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAg                                                                                                                   <  2:61245/150‑1 (MQ=255)
                             gCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTACGTGGTGGAAACCGGCGGACCGCATGGTGGACCATGGGGACGCCAGGTGATAACGCTTTTCACAAGTACTTCaa                                                                                                           >  2:302539/1‑149 (MQ=255)
                                 gCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCaa                                                                                                       >  1:836639/1‑150 (MQ=255)
                                  cTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAAc                                                                                                      <  2:836639/150‑1 (MQ=255)
                                                ttGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGa                                                                                        >  1:179275/1‑150 (MQ=255)
                                                            ttCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGaa                                                                                                                    <  2:660441/110‑1 (MQ=255)
                                                            gtCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGaa                                                                                                                    >  1:660441/2‑110 (MQ=255)
                                                                 aTTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAggcgggc                                                                       <  2:535946/150‑1 (MQ=255)
                                                                        tGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTAGCGGCAACGCTGCTTCAGGAGGGCGCACCcg                                                                >  1:639127/1‑150 (MQ=255)
                                                                         gTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCcgc                                                               <  1:813765/150‑1 (MQ=255)
                                                                               gCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACgg                                                                                                    <  1:238583/107‑1 (MQ=255)
                                                                               gCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCAAACgg                                                                                                    >  2:238583/1‑107 (MQ=255)
                                                                                cAGGCCTGGGATTCCTTTTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGGCACCATGTTAAAACGATTTTCAGAAGCACTTCAGTCACCAGAAGTCGCGGCGACGCTTCTTCAGGCGGGCGAACCCGCGCACCgg                                                        <  1:882768/150‑1 (MQ=255)
                                                                                   gCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAg                                                     <  2:65496/150‑1 (MQ=255)
                                                                                        ggATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACa                                                >  1:292811/1‑150 (MQ=255)
                                                                                                  ttAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGc                                      >  1:894092/1‑150 (MQ=255)
                                                                                                              tggAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACcgc                          >  1:519048/1‑150 (MQ=255)
                                                                                                                           cccGCATGGTAGACCATCGCGCCGCCAGGTGTAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACAcc             >  2:558117/1‑150 (MQ=255)
                                                                                                                              gCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACACCGCt          <  2:519048/150‑1 (MQ=255)
                                                                                                                                   gTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACACCGCTGCAAt     >  1:615258/1‑150 (MQ=255)
                                                                                                                                      gACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACACCGCTGCAATAcc  >  2:310184/1‑150 (MQ=255)
                                                                                                                                      gACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACACCGCAACCGCCAGTAGCACCGCGCCCAGACACACCGCTGCAATAcc  <  2:292811/150‑1 (MQ=255)
                                                                                                                                         cATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAgc                                              >  1:627505/1‑103 (MQ=255)
                                                                                                                                         cATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAgc                                              <  2:627505/103‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                      
CGGGCGGCAGGCCCAGTGGCATCGATGGAGCAGGCTCGGCACCTACAGTTGCGGGGGCAGTTCCGATTGGCCTGTATTGGCAGGCCTGGGATTCCCTGTTAGTTCCGTGGTGGAAACCGGCGGCCCGCATGGTAGACCATCGCGCCGCCAGGTGAAAACGCTTTTCAGAAGTACTTCAGCCAACGGATGTCGCGGCGACGCTGCTTCAGGCGGGCGAACCCGCGCACCGGCAGGTACAGCGCAACCGCCAGCAGCACCGCGCCCAGCCACACCGCTGCAATACC  >  NC_002947/2765134‑2765417

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