New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 240111 =26 (0.860)5 (0.170) 5/276 3.0 15.9% coding (406/1128 nt) PP_0188 uroporphyrin‑III C‑methyltransferase
?NC_002947 240158 = 28 (0.970)coding (453/1128 nt) PP_0188 uroporphyrin‑III C‑methyltransferase
Rejected: Frequency below/above cutoff threshold.

GCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGCCTTCGCCGCCCGCGAGCAACTGGCCCGCAGCCTGGCCACGCTCAACA  >  NC_002947/240011‑240303
                                                                                                                                                   |                                                                                                                                                 
ggcacccgggtcgctctgctcgcgcaggatctcgtcggcgccttcgaccagggctttggcactggtgatgtcctgcaacgccgaaaggcgcagggtagccAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCggctg                                                                                                                                                   >  3:137636‑M1/101‑148 (MQ=255)
  cAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCggctg                                                                                                                                                   <  2:50208/146‑1 (MQ=255)
   acccgggtcgctctgctcgcgcaggatctcgtcggcgccttcgaccagggctttggcactggtgatgtcctgcaacgccgaaaggcgcagggtagccAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCggctggct                                                                                                                                                <  1:47668‑M1/51‑1 (MQ=255)
    cccgggtcgctctgctcgcgcaggatctcgtcggcgccttcgaccagggctttggcactggtgatgtcctgcaacgccgaaaggcgcagggtagccAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTa                                                                                                                                               >  2:47668‑M1/97‑148 (MQ=255)
        ccAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGACGCCGAGCACCTCCTGCGGCTGGCTACCCt                                                                                                                                           >  5:89791/1‑148 (MQ=255)
          aGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTgc                                                                                                                                         <  5:89500/148‑1 (MQ=255)
          aGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTgc                                                                                                                                         <  6:164960/148‑1 (MQ=255)
                  agaagaCCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCg                                                                                                                                 >  6:69162/1‑148 (MQ=255)
                   gaagaCCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCgg                                                                                                                                >  4:82498/1‑148 (MQ=255)
                    aagaCCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGc                                                                                                                               >  7:105877/1‑148 (MQ=255)
                          gccgGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCa                                                                                                                         <  8:50508/148‑1 (MQ=255)
                            cgGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAgg                                                                                                                       <  6:52941/148‑1 (MQ=255)
                                      tAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAg                                                                                                             >  6:52244/1‑148 (MQ=255)
                                      tAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAg                                                                                                             >  5:37364/1‑148 (MQ=255)
                                       aGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGt                                                                                                            <  3:129087/148‑1 (MQ=255)
                                       aGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGt                                                                                                            <  7:149035/148‑1 (MQ=255)
                                           cAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGcc                                                                                                         >  3:48604/1‑147 (MQ=255)
                                                 cAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCCCCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAACGccc                                                                                                  >  1:88150/1‑148 (MQ=255)
                                                      cGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTc                                                                                             <  6:108411/148‑1 (MQ=255)
                                                                           gcgcCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTg                                                                        >  1:10155/1‑148 (MQ=255)
                                                                                tGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGa                                                                   >  6:31986/1‑148 (MQ=255)
                                                                                  gAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGa                                                                   >  5:91356/1‑146 (MQ=255)
                                                                                              ggAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCgg                                                     >  7:3178/1‑148 (MQ=255)
                                                                                               tagccagccgcaggaggtgctcggcttcggccaggcgccactctttgcggctGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCggg                                                    >  2:47668‑M2/53‑148 (MQ=255)
                                                                                                agccagccgcaggaggtgctcggcttcggccaggcgccactctttgcggctGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGt                                                   <  1:47668‑M2/97‑1 (MQ=255)
                                                                                                   cagccgcaggaggtgctcggcttcggccaggcgccactctttgcggctGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGcc                                                >  3:137636‑M2/49‑148 (MQ=255)
                                                                                                    agccgcaggaggtgctcggcttcggccaggcgccactctttgcggctGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGCCt                                               >  6:56295‑M2/48‑148 (MQ=255)
                                                                                                     gccgcaggaggtgctcggcttcggccaggcgccactctttgcggctGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGCCtt                                              >  3:73875‑M2/47‑148 (MQ=255)
                                                                                                                 ggCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGCCTTCGCCGCCCGCGa                                  <  1:37526/148‑1 (MQ=255)
                                                                                                                       tGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGCCTTCGCCGCCCGCGAGCAACt                            <  6:37364/148‑1 (MQ=255)
                                                                                                                          ccGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCTCCGACGAGATCCTGCGCTAGCAGAGCGACCCGGGTGCCTTCGCCGCCCGCGAGCAACTg                           <  5:127824/146‑1 (MQ=255)
                                                                                                                           cGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGCCTTCGCCGCCCGCGAGCAACTGGcc                        >  8:80536/1‑148 (MQ=255)
                                                                                                                           cGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGCCTTCGCCGCCCGCGAGCAACTGGcc                        >  4:85752/1‑148 (MQ=255)
                                                                                                                                    gCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGCCTTCGCCGCCCGCGAGCAACTGGCCCGCAGCCTg               <  6:91356/148‑1 (MQ=255)
                                                                                                                                        ctcctGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGCCTTCGCCGCCCGCGAGCAACTGGCCCGCAGCCTGGCCa           >  5:135013/1‑148 (MQ=255)
                                                                                                                                         tcctGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGCCTTCGCCGCCCGCGAGCAACTGGCCCGCAGCCTGGCCAc          >  1:10245/1‑148 (MQ=255)
                                                                                                                                            tgcggctGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGCCTTCGCCGCCCGCGAGCAACTGGCCCGCAGCCTGGCCACGc        >  5:35506‑M2/8‑147 (MQ=255)
                                                                                                                                                 ctGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGCCTTCGCCGCCCGCGAGCAACTGGCCCGCAGCCTGGCCACGCTCAACa  <  5:52244‑M2/146‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCCAGCCGCCAGCGAGCTAGAAGACCGCCGGCGTTTGGTAGCCCAGTTGCAGGGCGACCAGCAGCGTCTCAGCCAGCGCCTGGAAACCGTGCTGGGAGAAAGCCGCAAAGAGTGGCGCCTGGCCGAAGCCGAGCACCTCCTGCGGCTGGCTACCCTGCGCCTTTCGGCGTTGCAGGACATCACCAGTGCCAAAGCCCTGGTCGAAGGCGCCGACGAGATCCTGCGCGAGCAGAGCGACCCGGGTGCCTTCGCCGCCCGCGAGCAACTGGCCCGCAGCCTGGCCACGCTCAACA  >  NC_002947/240011‑240303

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