New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 523868710 (0.620)4 (0.250) 4/250 NT 26.0% coding (151/204 nt) PP_5701 hypothetical protein
?NC_002947 = 5238796 13 (0.820)intergenic (+56/‑225) PP_5701/PP_4616 hypothetical protein/hypothetical protein

GAAGGCGACCGATTGCCGCAAGGCACCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCAC  >  NC_002947/5238540‑5238836
                                                                                                                                                   |                                                                                                                                                     
gAAGGCGACCGATTGCCGCAAGGCACCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCGACTTGAGTCCCCGCGGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGaa                                                                                                                                                     >  1:144768/1‑150 (MQ=255)
               ccGCAAGGCACCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTca                                                                                                                                      >  1:177049/1‑149 (MQ=255)
                                      cACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCg                                                                                                               <  2:326431/150‑1 (MQ=255)
                                               tGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGcaggtccggccctttcgcgggcacgcccgctcccacaggtactgtgcaa                                                                                                      >  2:21464‑M1/1‑101 (MQ=255)
                                                                         ttGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGcaggtccggccctttcgcgggcacgcccgctcccacaggtactgtgcaagtcttgactactgtggtgaccgctgt                                                                            >  2:17459‑M1/1‑75 (MQ=255)
                                                                             gTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGa                                                                        <  2:185672/150‑1 (MQ=255)
                                                                              tCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGccc                                                                                                                           <  1:266480/98‑1 (MQ=255)
                                                                              tCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGccc                                                                                                                           >  2:266480/1‑98 (MQ=255)
                                                                                 gcgcTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGcaggtccggccctttcgcgggcacgcccgctcccacaggtactgtgcaagtcttgactactgtggtgaccgctgtgggagcgg                                                                    <  1:17459‑M1/150‑84 (MQ=255)
                                                                                   gcTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCg                                                                                                               >  2:404538/1‑105 (MQ=255)
                                                                                   gcTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCg                                                                                                               <  1:404538/105‑1 (MQ=255)
                                                                                       gTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGcaggtccggccctttcgcgggcacgcccgctcccacaggtactgtgcaagtcttgactactgtggtgaccgctgtgggagcgggcgagc                                                              <  1:21464‑M1/150‑90 (MQ=255)
                                                                                       gTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGc                                                              <  2:144768/150‑1 (MQ=255)
                                                                                        ttCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGcc                                                             >  2:37870/1‑150 (MQ=255)
                                                                                                           tGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCACTCCCACAGCGGGAAACACAGAAGTCAAGACTTGCACAGTACCTATGGGAGCGGGCGTGCGCGCGAAAGGGCGGGCCCTg                                          >  1:53049/1‑150 (MQ=255)
                                                                                                            gAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGc                                         >  1:403400‑M2/1‑149 (MQ=255)
                                                                                                                                                   gcaggtccggccctttcgcgggcacgcccgctcccacaggtactgtgcaagtcttgactactgtggtgaccgctgtgggagcgggcgagcccgcgaagaggccggcccttccggcacgttaaccagccttagtgcgaagttcctgggctc  >  1:403400‑M1/1‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GAAGGCGACCGATTGCCGCAAGGCACCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCAC  >  NC_002947/5238540‑5238836

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