New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 95322326 (0.960)6 (0.230) 6/288 2.3 19.8% coding (29/624 nt) cyoC cytochrome bo terminal oxidase subunit III
?NC_002947 = 953241 23 (0.870)coding (47/624 nt) cyoC cytochrome bo terminal oxidase subunit III
Rejected: Frequency below/above cutoff threshold.

CGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCGGGTCATGAC  >  NC_002947/953092‑953389
                                                                                                                                                     |                                                                                                                                                    
cGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATgg                                                                                                                                                      <  1:445114/150‑1 (MQ=255)
 gaagaaGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGg                                                                                                                                                     >  2:237456‑M2/1‑149 (MQ=255)
      agTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATga                                                                                                                                                >  1:545013/1‑150 (MQ=255)
               cATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGaccacca                                                                                                                                       >  2:83508/1‑150 (MQ=255)
                   gAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACcaccaccac                                                                                                                                   <  1:522997/150‑1 (MQ=255)
                          aGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCgg                                                                                                                            <  1:483023/150‑1 (MQ=255)
                          aGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCgg                                                                                                                            <  1:483027/150‑1 (MQ=255)
                              tATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCa                                                                                                                        >  1:521524/1‑150 (MQ=255)
                                tGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGa                                                                                                                      >  2:571563/1‑150 (MQ=255)
                                             cAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCTGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTggg                                                                                                         >  1:69916/1‑150 (MQ=255)
                                                          ctgcctgcCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTAc                                                                                            <  1:313931/150‑1 (MQ=255)
                                                          ctgcctgcCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTAc                                                                                            <  1:313932/150‑1 (MQ=255)
                                                                   ggCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGAcc                                                                                   >  2:93109/1‑150 (MQ=255)
                                                                   ggCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGAcc                                                                                   >  2:93107/1‑150 (MQ=255)
                                                                                              gtgtAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTc                                                        >  2:528983/1‑150 (MQ=255)
                                                                                              gtgtAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTc                                                        >  2:528982/1‑150 (MQ=255)
                                                                                                   aTCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCAc                                                   <  1:413658/150‑1 (MQ=255)
                                                                                                       aTGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTAc                                               >  1:589884/1‑150 (MQ=255)
                                                                                                          tCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTAcgcg                                            >  1:3039/1‑150 (MQ=255)
                                                                                                             aGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTg                                         >  2:534164/1‑150 (MQ=255)
                                                                                                              gTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGgcagcgccgtgcattacttgactggacattgattacacccgctc                                                                                                          <  1:447118‑M2/84‑45 (MQ=255)
                                                                                                              gTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGgcagcgccgtgcattacttgactggacattgattacacccgctc                                                                                                          >  2:447118‑M2/1‑40 (MQ=255)
                                                                                                                 aaGTACTGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGt                                     >  2:570761/1‑150 (MQ=255)
                                                                                                                      aTGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGc                                >  1:212627/1‑150 (MQ=255)
                                                                                                                             gcgcTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTg                         >  1:88171/1‑150 (MQ=255)
                                                                                                                               gctgctCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGcc                       <  1:391485/150‑1 (MQ=255)
                                                                                                                                 tgctCATGGTCACGACCATGGgcagcgccgtgcattacttgactggacattcattacacccgctcaaacgattagacacgtgcgccggcaggcagagcacctgctttggccggcgctttcatacgctcgccttcgatgcgcgccacttct                     <  1:205507‑M2/150‑130 (MQ=255)
                                                                                                                                  gctCATGGTCACGACCATGGgcagcgccgtgcattacttgactggacattgattacacccgctcaaacgattcgacacgtgcgccggcaggcagagcacctgctttggccagcgctttcatacgctcgccttcgatg                                 >  2:451413‑M2/1‑20 (MQ=255)
                                                                                                                                  gctCATGGTCACGACCATGGgcagcgccgtgcattacttgactggacattgattacacccgctcaaacgattcgacacgtgcgccggcaggcagagcacctgctttggccagcgctttcatacgctcgccttcgatg                                 <  1:451413‑M2/137‑118 (MQ=255)
                                                                                                                                    tCATGGTCACGACCATGAgcagcgccgtgcattacttgactggacattgattacacccgctcaaacgattcgacacgtgcgccggcaggcagagcacctgctttggccagcgctttcatacgctcgccttcgatgcgcgccacttcttcg                  <  1:527312‑M2/150‑133 (MQ=255)
                                                                                                                                                    gggCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCGGGTCATGAc  >  1:452607/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGAAGAAGTGGCGCGCATCGAAGGCGAGCGTATGAAAGCGCTGGCCAAAGCAGGTGCTCTGCCTGCCGGCGCACGTGTCGAATCGTTTGAGCGGGTGTAATCAATGTCCAGTCAAGTAATGCACGGCGCTGCTCATGGTCACGACCATGGGCATGACGACCACCACCACGACTCGGGCCAGATGACCGTACTGGGCTTCTGGCTGTACCTGATGACCGACTGCATCCTGTTTGCGTCGCTCTTCGCCACCTACGCGGTGCTGTCCGGCAGTTTTGCCGGCGGCCCGTCGGGTCATGAC  >  NC_002947/953092‑953389

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