New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 525739927 (1.170)3 (0.140) 3/234 NT 11.7% coding (451/909 nt) PP_4634 alpha/beta family hydrolase
?NC_002947 = 5257429 20 (0.920)coding (421/909 nt) PP_4634 alpha/beta family hydrolase

GCGTGGTCCAGTCGTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCTGGCCAGCCGATGGTGCCAGGCCTACTGCACAACCGCCGAAACTGTGCCCCACCACGTCCACCGGCTGACCTGGAAAATGCGCTGCCG  >  NC_002947/5257261‑5257542
                                                                                                                                          |                                                                                                                                               
gcgTGGTCCAGTCGTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAgc                                                                                                                                      >  1:133735/1‑150 (MQ=255)
  gTGGTCCAGTCGTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAccgc                                                                                                                                    <  2:57764/150‑1 (MQ=255)
     gTCCAGTCGTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAg                                                                                                                                 <  1:425304/150‑1 (MQ=255)
        cAGTCGTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAg                                                                                                                              >  2:242651/1‑150 (MQ=255)
             gTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGaa                                                                                                                         >  2:273771/1‑150 (MQ=255)
             gTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGaa                                                                                                                         <  1:242651/150‑1 (MQ=255)
             gTGCAACACGACAGATGGGGTATCTTACAGCTACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACTTGCCACTTGCCGAACAGCTGCCCGCGCTGGGTGGCCGCATAGTAGCGCCAGTAGGCGaa                                                                                                                         >  2:195172/1‑150 (MQ=255)
                     cGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGc                                                                                                                            >  2:60000/1‑139 (MQ=255)
                     cGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGc                                                                                                                            <  1:60000/139‑1 (MQ=255)
                          gATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGACGAACAGCTTCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCgccaacc                                                                                                            <  2:391408‑M2/150‑8 (MQ=255)
                           ctgacgtctCTGCCAGCCACCCCGGGGGCTGTCCGTGAAAATAGCCGAATACCCCGGTCAGCAGCGGCATGACCACATGCCACGTGCCGAACAGCTGCCATCGCTGGTTGGCgcgcagttcgcctactggcgctactatgcggccaacca                                                                                                           <  2:86190‑M1/146‑39 (MQ=255)
                              gCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCgcgcagttcgcctactggcgctactatgcggccaaccagag                                                                                                        >  1:86190‑M1/1‑109 (MQ=255)
                               cGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCAc                                                                                                       >  2:143362/1‑150 (MQ=255)
                                gaggcttccAGCAACCCCCGGGGCTGGCCCGGAAAATAGCCGGATAACCGGGTCAGCAGCGGCATGCCCACATGCCACCTGCCTAACAGCTGCCAGCGCTGGTTGGCgcgcagttcgcctaattgcgctactatgcggcaaaccagcgct                                                                                                      <  2:47956‑M1/146‑44 (MQ=255)
                                     ttCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGAc                                                                                                    <  2:408597/147‑1 (MQ=255)
                                     ttCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGAc                                                                                                    >  1:408597/1‑147 (MQ=255)
                                            cACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACAAGATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATc                                                                                          <  1:6551/150‑1 (MQ=255)
                                             aCCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCt                                                                                         <  1:143362/150‑1 (MQ=255)
                                                            cGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGAccc                                                                                                  <  1:403786/126‑1 (MQ=255)
                                                            cGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGAccc                                                                                                  >  2:403786/1‑126 (MQ=255)
                                                            cGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCTGGCCAGCCGATGGTg                                                                          <  1:457317/150‑1 (MQ=255)
                                                            cGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCTGGCCAGCCGATGGTg                                                                          >  1:253475/1‑150 (MQ=255)
                                                                 aaaTGGCCGAATCGCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCTGGCCAGCCGATGGTGCCAgg                                                                     <  2:436734/150‑1 (MQ=255)
                                                                       ccGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCTGGCCAGCCGATGGt                                                                           <  1:345737/138‑1 (MQ=255)
                                                                       ccGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCTGGCCAGCCGATGGt                                                                           >  2:345737/1‑138 (MQ=255)
                                                                       ccGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCTGGCCAGCCGATGGTGCCAGGCCTACt                                                               <  2:253475/150‑1 (MQ=255)
                                                                                          gcGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCTGGCCAGCCGATGGTGCCAGGCCTACTGCACAACCGCCGAAACtgt                                            <  2:542165/150‑1 (MQ=255)
                                                                                          gcGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCTGGCCAGCCGATGGGGCCAGGCCTACTGCACAACCGCCGAAACtgt                                            <  1:273771/150‑1 (MQ=255)
                                                                                                             cTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCTGGCCAGCCGATGGTGCCAGGCCTACTGCACAACCGCCGAAACTGTGCCCCACCACGTCCACCgg                         >  1:202789/1‑150 (MQ=255)
                                                                                                                              aGCGCTGGTTTGCCGCATAGTAGCGCAATTAGGCGAACTGCGCgccaaccagcgctggcagctgttaggcaggtggcatgtgggcatgccgctgctgacccggttatccggctattttccgggccagccccgggggttgctggaagcctc        <  2:47956‑M2/150‑108 (MQ=255)
                                                                                                                               gcgcTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCTGGCCAGCCGATGGTGCCAGGCCTACTGCACAACCGCCGAAACTGTGCCCCACCAAGTCCACCGGCTGACCTGGAAAATgcgc       <  1:391408/150‑1 (MQ=255)
                                                                                                                                ctctgGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCgccaaccagcgctggcagctgttcggcaagtggcatgtggtcatgccgctgctgacccgggtattcggctattttcccggcaagcgcctggggtggctggaagacacgc      >  1:86190‑M2/3‑41 (MQ=255)
                                                                                                                                   tGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCgccaaccagcgatggcagctgttcggcacgtggcatgtggtcatgccgctgctgaccggggtattcggctattttcacggacagcccccggggtggctggcagagacgtcag   <  2:86190‑M2/150‑113 (MQ=255)
                                                                                                                                   tGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCTGGCCAGCCGATGGTGCCAGGCCTACTGCACAACCGCCGAAACTGTGCCCCACCACGTCCACCGGCTGACCTGGAAAATGCGCTGcc   >  1:337088/1‑150 (MQ=255)
                                                                                                                                    ggTTGGCgcgcagttcgcctactggcgctactatgcggccaaccagcgctggaagctgttcgtcaagtggcatgtggtcatgccgctgctgacccgggtattcggctattttcccggcaagcgcctggggtggctggaagacacgccatc  <  2:391408‑M1/150‑144 (MQ=255)
                                                                                                                                          |                                                                                                                                               
GCGTGGTCCAGTCGTGCACCACGCCGGATGGCGTGTCTTCCAGCCACCCCAGGCGCTTGCCGGGAAAATAGCCGAATACCCGGGTCAGCAGCGGCATGACCACATGCCACTTGCCGAACAGCTGCCAGCGCTGGTTGGCCGCATAGTAGCGCCAGTAGGCGAACTGCGCGCCAACCATCACGACCCGACGGATCTGGCCAGCCGATGGTGCCAGGCCTACTGCACAACCGCCGAAACTGTGCCCCACCACGTCCACCGGCTGACCTGGAAAATGCGCTGCCG  >  NC_002947/5257261‑5257542

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.