New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5957953 =18 (0.960)4 (0.230) 4/230 NT 18.3% coding (1301/2856 nt) cyaA adenylate cyclase
?NC_002947 5957991 = 19 (1.080)coding (1339/2856 nt) cyaA adenylate cyclase

CGCAGCCAATGGAAAGTGCAGCAGGTCGCCGTCGAGCGGCGTGAGCTGGTGGCCGAGCTTAACCTCAGTTATCGCTTCCTCAGCCAGTTCGCCCGCACCCAGAACGCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGAAGACACCCTGACGCTGGTCCAGGCACCCAACCGCAAGGAGCCGGGCAGCTACCACTGGGAGCTG  >  NC_002947/5957808‑5958098
                                                                                                                                                 |                                                                                                                                                 
agtgctgtacctcatgggtgctcaggttgccgctgtacagctcccagtggtagctgcccggctccttgcggttgggtgcctggaccagcgtcagggtgtcttcggcgaggtcaggggcgatgccggggttaatcacttcgaccttGCCgg                                                                                                                                               <  1:437747‑M1/5‑1 (MQ=255)
 gCAGTCAATGGAAAGTGCAGCAGGTCGCCGTCGAGCGGCGTGAGCTGGTGGCCGAGCTTAACCTCAGTTATCGCTTCCTCAGCCAGTTCGCCCGCACCCAGAACGCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGgc                                                                                                                                              <  1:213771/150‑1 (MQ=255)
 gCAGCCAATGGAAAGTGCAGCAGGTCGCCGTCGAGCGGCGTGAGCTGGTGGCCGAGCTTAACCTCAGTTATCGCTTCCTCAGCCAGTTCGCCCGCACCCAGAACGCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGgc                                                                                                                                              <  2:405550/150‑1 (MQ=255)
                       ggTCGCCGTCGAGCGGCGTGAGCTGGTGGCCGAGCTTAACCTCAGTTATCGCTTCCTCAGCCAGTTCGCCCGCACCCAGAACGCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAgcg                                                                                                                        >  2:207962/1‑150 (MQ=255)
                         tCGCCGTCGAGCGGCGTGAGCTGGTGGCCGAGCTTAACCTCAGTTATCGCTTCCTCAGCCAGTTCGCCCGCACCCAGAACGCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGa                                                                                                                           <  1:261853/145‑1 (MQ=255)
                         tCGCCGTCGAGCGGCGTGAGCTGGTGGCCGAGCTTAACCTCAGTTATCGCTTCCTCAGCCAGTTCGCCCGCACCCAGAACGCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGa                                                                                                                           >  2:261853/1‑145 (MQ=255)
                                        ctcccagtggtagctgcccggctccttgcggttgggtgcctggaccagcgtcaaggtgtcttcggcgaggtcaggggcgatgccggggttaatcacttcgaccttGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAgg                                                                                                          >  1:331532‑M1/106‑147 (MQ=255)
                                        ctcccagtggtagctgcccggctccttgcggttgggtgcctggaccagcgtcaaggtgtcttcggcgaggtcaggggcgatgccggggttaatcacttcgaccttGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAgg                                                                                                          <  2:331532‑M1/42‑1 (MQ=255)
                                            gCTGGTGGCCGAGCTTAACCTCAGTTATCGCTTCCTCAGCCAGTTCGCCCGCACCCAGAACGCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGt                                                                                                   <  1:192936/150‑1 (MQ=255)
                                                          cggctccttgcggttgggtgcctggaccagcgtcagggtgtcttcggcgaggtcaggggcgatgccggggttaatcacttcgaccttGCCGGCGCCTGTACGCCGCCTTCGAgcgccgcgcc                                                                                                                 >  1:326231‑M1/88‑122 (MQ=255)
                                                          cggctccttgcggttgggtgcctggaccagcgtcagggtgtcttcggcgaggtcaggggcgatgccggggttaatcacttcgaccttGCCGGCGCCTGTACGCCGCCTTCGAgcgccgcgcc                                                                                                                 <  2:326231‑M1/35‑1 (MQ=255)
                                                               cTCAGTCATCGCTTACTCACCCAGCTCGCCCGCACTCAGAACGCAAGCAGCCGCGCCGACCAGCGCGACCTGAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCCCCGGCCAGGTCGAAGTGATTCACCCCGGCATcccc                                                                                >  2:246879/1‑150 (MQ=255)
                                                                              cTCAGCCAGTTCGCCCGCACCCAGAACGCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGaa                                                                 <  2:72755/150‑1 (MQ=255)
                                                                              cTCAGCCAGTTCGCCCGCACCCAGAACGCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGaa                                                                 <  2:278522/150‑1 (MQ=255)
                                                                                     aGTTCGCCCGCACCCAGAACGCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGAAGACAccc                                                          >  1:171188/1‑150 (MQ=255)
                                                                                     aGTTCGCCCGCACCCAGAACGCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGAAGACAccc                                                          <  1:207962/150‑1 (MQ=255)
                                                                                                  ccAGAACGCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGAAGACACCCTGACGCTGGTCCa                                             >  1:75466/1‑150 (MQ=255)
                                                                                                         gCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGAAGACACCCTGACGCTg                                                  >  1:388402/1‑138 (MQ=255)
                                                                                                         gCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGAAGACACCCTGACGCTg                                                  <  2:388402/138‑1 (MQ=255)
                                                                                                                  cgcgCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGAAGACACCCTGACGCTGGTCCAGGCACCCAACCGCAAg                             <  2:75466/150‑1 (MQ=255)
                                                                                                                     gcCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGAAGACACCCTGACGCTGGTCCAGGCACCCAACCGCAAGGAg                          >  2:100785/1‑150 (MQ=255)
                                                                                                                              cgcgACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGAAGACACCCTGACGCTGGTCCAGGCACCCAACCGCAAGGAGCCGGGCAGc                 <  1:100785/150‑1 (MQ=255)
                                                                                                                                   ccTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGAAGACACCCTGACGCTGGTCCAGGCACCCAACCGCAAGGAGCCGGGCAGCTACCa            >  2:444576/1‑150 (MQ=255)
                                                                                                                                             cTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGAAGACACCCTGACGCTGGTCCAGGCACCCAACCGCAAGGAGCCGGGCAGCTACCACTGGGAGCTg  >  2:477130/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CGCAGCCAATGGAAAGTGCAGCAGGTCGCCGTCGAGCGGCGTGAGCTGGTGGCCGAGCTTAACCTCAGTTATCGCTTCCTCAGCCAGTTCGCCCGCACCCAGAACGCCAGCAGCCGCGCCGACCAGCGCGACCTCAATGTGCTGGGCCGGCGCCTGTACGCCGCCTTCGAGCGCCGCGCCGGCAAGGTCGAAGTGATTAACCCCGGCATCGCCCCTGACCTCGCCGAAGACACCCTGACGCTGGTCCAGGCACCCAACCGCAAGGAGCCGGGCAGCTACCACTGGGAGCTG  >  NC_002947/5957808‑5958098

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