New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP034943 = 2653801463 (1.320)3 (0.010) 3/246 NT 0.67% coding (466/1161 nt) safA spore coat assembly protein SafA
?NZ_CP034943 2653808 = 440 (1.290)coding (459/1161 nt) safA spore coat assembly protein SafA
Rejected: Frequency below/above cutoff threshold.

GCTGAACAGCCGGCATATTCGGAACTGTGTGATGAACGTTTTCTTCTTCCATACCGCCTACCTCTGGCGCCTTTGGCATATTAGGCATATTCGGATAATTTGCATTTTCCATATTACTCATAGCCTCCTGTTGTGGAAATTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP034943/2653660‑2653801
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttgGTCTTGCATGTGATACATATGATCCATATTTCCGTCATGATGTTTTTTCGGCTCCTCCGGTTTAGGGGGCGACCAAGGCTGGAAAAGCTGTTTCATATCATAATAATCGTTCACATCAGCTTCAGGATACACATCTTCC  >  NZ_CP034943/2653808‑2653946
                                                                                                                                                                                                                                                                                         
GCTGAACAGCCGGCATATTCGGAACTGTGTGATGAACGTTTTCTTCTTCCATACCGCCTACCTCTGGCGCCTTTGGCATATTAGGCATATTCGGATAATTTGCATTTTCCATATTACTCATAGCCTCCTGTTGTGGAAAT                                                                                                                                               >  1:2420755/1‑140
GCTGAACAGCCGGCATATTCGGAACTGTGTGATGAACGTTTTCTTCTTCCATACCGCCTACCTCTGGCGCCTTTGGCATATTAGGCATATTCGGATAATTTGCATTTTCCATATTACTCATAGCCTCCTGTTGTGGAAAT                                                                                                                                               <  2:1646816/140‑1
GCTGAACAGCCGGCATATTCGGAACTGTGTGATGAACGTTTTCTTCTTCCATACCGCCTACCTCTGGCGCCTTTGGCATATTAGGCATATTCGGATAATTTGCATTTTCCATATTACTCATAGCCTCCTGTTGTGGAAAT                                                                                                                                               <  2:9286/140‑1
                                                                       TTTGGCATATTAGGCATATTCGGATAATTTGCATTTTCCATATTACTCATAGCCTCCTGTTGTGGAAATTGGTCTTGCATATGATACATATGATCCATATTTCCGTCATGATGTTTTTTCGGCTCCT                                                                                     >  1:2131088/1‑127
                                                                       TTTGGCATATTAGGCATATTCGGATAATTTGCATTTTCCATATTACTCATAGCCTCCTGTTGTGGAAATTGGTCTTGCATATGATACATATGATCCATATTTCCGTCATGATGTTTTTTCGGCTCCT                                                                                     <  2:2131088/127‑1
                                                                                              ATAATTTGCATTTTCCATATTACTCATAGCCTCCTGTTGTGGAAATTGGTCTTGCATATGATACATATGATCCATATTTCCGTCATGATGTTTTTTCGGCTCCTCCGGTTTAGGAGGCGACCAAGGCTGGAAAAGCTGTT                                                 >  2:1188648/1‑140
                                                                                                                                             GGTCTTGCATGTGATACATATGATCCATATTTCCGTCATGATGTTTTTTCGGCTCCTCCGGTTTAGGGGGCGACCAAGGCTGGAAAAGCTGTTTCATATCATAATAATCGTTCACATCAGCTTCAGGATACACATCTTCC  <  2:4618165/140‑1
                                                                                                                                             GGTCTTGCATGTGAAACATATGATCCATATTTCCGTCATGATGTTCTTTCGGCTCCTCCGGGTTAGGGGGCGACCAAGGCTGGAAAAGCTGTTTCATATCATAATAATCGTTCACATCAGCTTCAGGATACACATCTTCC  <  2:3453266/140‑1
                                                                                                                                                                                                                                                                                         
GCTGAACAGCCGGCATATTCGGAACTGTGTGATGAACGTTTTCTTCTTCCATACCGCCTACCTCTGGCGCCTTTGGCATATTAGGCATATTCGGATAATTTGCATTTTCCATATTACTCATAGCCTCCTGTTGTGGAAATTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP034943/2653660‑2653801
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttgGTCTTGCATGTGATACATATGATCCATATTTCCGTCATGATGTTTTTTCGGCTCCTCCGGTTTAGGGGGCGACCAAGGCTGGAAAAGCTGTTTCATATCATAATAATCGTTCACATCAGCTTCAGGATACACATCTTCC  >  NZ_CP034943/2653808‑2653946

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