New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 581208182 (1.100)3 (0.040) 3/244 NT 3.6% coding (1353/2220 nt) priA primosome assembly protein
?NC_002947 = 5812115 81 (1.130)coding (1387/2220 nt) priA primosome assembly protein
Rejected: Frequency below/above cutoff threshold.

TGCCACTGGACAGCGGCATCAGTGGCCCGCTGCAGCAGGCTATCCGACAGACCCTGGAGGCTGGCCAACAAGTGCTGGTGTTCCTCAACCGCCGCGGCTTTGCCCCAACCTTGCTGTGCCACGATTGCGGCTGGCTGTCCGAATGCCCACGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5811930‑5812081
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cacgcCGGAGCGCTGGTGCACGGTCATGCGCGCATCGCAGCGTGGGCATTCGGACAGCCAGCCGCAATCGTGGCACAGCAAGGTTGGGGCAAAGCCGCGGCGGTTGAGGAACACCAGCACTTGTTGGCCAGCCTCCAGGGTCTGTCGGATAGC  <  NC_002947/5812115‑5811968
                                                                                                                                                                                                                                                                                                            
TGCCACTGGACAGCGGCATCAGTGGCCCGCTGCAGCAGGCTATCCGACAGACCCTGGAGGCTGGCCAACAAGTGCTGGTGTGCCTCAACCGCCGCGGCTTTGCCCCAACCTTGCTGTGCCACGATTGCGGCTGGATGTCCGATTGCCCAC                                                                                                                                                        >  2:526339/1‑150
                                          TCCGACAGACCCTGGAGGCTGGCCAACAAGTGCTGGTGTTCCTCAACCGCCGCGGCTTTGCCCCAACCTTGCTGTGCCACGATTGCGGCTGGCTGTCCGAATGCCCACGCCGGAGCGCTGGTGCACGGTCATGCGCGCATCGCAGCGTGG                                                                                                              <  2:1275672/150‑1
                                                                                    TCAACCGCCGCGGCTTTGCCCCAACCTTGCTGTGCCACGATTGCGGCTGGCTGTCCGAATGCCCACGCCGGAGCGCTGGTGCACGGTCATGCGCGCATC                                                                                                                       <  1:923579/99‑1
                                                                                    TCAACCGCCGCGGCTTTGCCCCAACCTTGCTGTGCCACGATTGCGGCTGGCTGTCCGAATGCCCACGCCGGAGCGCTGGTGCACGGTCATGCGCGCATC                                                                                                                       >  2:923579/1‑99
                                                                                                                                                      GCCGGAGCGCTGGTGCACGGTTATACGCGCATCGCAGCGTGGGCATTCGGACAGCCAGCCGCAGTCGTGGCACAGCAAGGTTGGGGCAAAGCCGCGGCGGTTGAGGAACACCAGCACTTGTTGGCCAGCCTCCAGGGTCTGTCGGATAGC  <  2:962401/150‑1
                                                                                                                                                      GCCGGAGCGCTGGTGCACGGTCATGCGCGCATCGCAGCGTGGGCATTCGGACAGCCAGCCGCAATCGTGGCACAGCAAGGTTGGGGCAAAGCCGCGGCGGTTGAGGAACACCAGCACTTGTTGGCCAGCCTCCAGGGTCTGTCGGATAGC  >  2:608738/1‑150
                                                                                                                                                      GCCGGAGCGCTGGTGCACGGTCATGCGCGCATCGCAGCGTGGGCATTCGGACAGCCAGCCGCAATCGTGGCACAGCAAGGTTGGGGCAAAGCCGCGGCGGTTGAGGAACACCAGCACTTGTTG                             >  1:1637794/1‑123
                                                                                                                                                      GCCGGAGCGCTGGTGCACGGTCATGCGCGCATCGCAGCGTGGGCATTCGGACAGCCAGCCGCAATCGTGGCACAGCAAGGTTGGGGCAAAGCCGCGGCGGTTGAGGAACACCAGCACTTGTTG                             <  2:1637794/123‑1
                                                                                                                                                                                                                                                                                                            
TGCCACTGGACAGCGGCATCAGTGGCCCGCTGCAGCAGGCTATCCGACAGACCCTGGAGGCTGGCCAACAAGTGCTGGTGTTCCTCAACCGCCGCGGCTTTGCCCCAACCTTGCTGTGCCACGATTGCGGCTGGCTGTCCGAATGCCCACGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5811930‑5812081
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cacgcCGGAGCGCTGGTGCACGGTCATGCGCGCATCGCAGCGTGGGCATTCGGACAGCCAGCCGCAATCGTGGCACAGCAAGGTTGGGGCAAAGCCGCGGCGGTTGAGGAACACCAGCACTTGTTGGCCAGCCTCCAGGGTCTGTCGGATAGC  <  NC_002947/5812115‑5811968

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