New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP034943 = 1328839323 (1.240)7 (0.030) 5/252 NT 2.1% coding (315/1359 nt) htrA serine protease HtrA
?NZ_CP034943 1328847 = 326 (1.250)coding (307/1359 nt) htrA serine protease HtrA
Rejected: Frequency below/above cutoff threshold.

AACAAAGAGCTGTTTGACTGAGCCTGAAGATTTGTAATGCCGACAATCGCCGGTGAAATGTCTTCCACCATATCTGAGATTTTAGAAGCATCCTCGCTCTTGAATGCAGATGAGCTGCTTGAGCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP034943/1328714‑1328839
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTTTGCTAGATTCAGAGGACGTGCTTTTAGCTGTAACCGATTGTGTTTGCTGCTGGCTGGATGATTGTTTTGAAGTGTCCTGAGCATCATGGTC  >  NZ_CP034943/1328847‑1328940
                                                                                                                                                                                                                            
AACAAAGAGCTGTTTGATTGTGCCTGAAGATTTGTAATACCGACAATCGCTGGTGAAAGGTCTTCTACCATATCTGAGATTTTAGAAGAGTCCTCGCTCTTGAATGCAGATGAGCTGCTTGAGCTTTTTT‑‑TAGATTCAGA                                                                                >  2:3191931/1‑140
AACAAAGAGCTGTTTGATTGTGCCTGAAGATTTGTAATACCGACAATCGCTGGTGAAAGGTCTTCTACCATATCTGAGATTTTAGAAGAGTCCTCGCTCTTGAATGCAGATGAGCTGCTTGAGCTTTTTT‑‑TAGATTCAGA                                                                                >  2:4251800/1‑140
                                              TCGCTGGTGAAAGGTCTTCTACCATATCTGAGATTTTAGAAGAGTCCTCGCTCTTGAATGCAGATGAGCTGCTTGAGCTTTTTT‑‑TAGATTCAGAGGAGGTGCTTGTTGCTGTAACAGATTGCGTTTGCTGCTGGCTGGAT                                  <  2:1004375/140‑1
                                                                      TATCTGAGATTTTAGAAGAGTCCTCGCTCTTGAATGCAGATGAGCTGCTTGAGCTTTTTT‑‑TAGATTCAGAGGAGGTGCTTGTTGCTGTAACAGATT                                                      <  1:603865/96‑1
                                                                      TATCTGAGATTTTAGAAGAGTCCTCGCTCTTGAATGCAGATGAGCTGCTTGAGCTTTTTT‑‑TAGATTCAGAGGAGGTGCTTGTTGCTGTAACAGATT                                                      >  2:603865/1‑96
                                                                              ATTTTAGAAGAGTCCTCGCTCTTGAATGCAGATGAGCTGCTTGAGCTTTTTT‑‑TAGATTCAGAGGAGGTGCTTGTTGCTGTAACAGATTGCGTTTGCTGCTGGCTGGATGATTGTTTTGCAGTGTCCTGAGAATCATGGTC  <  1:3191931/140‑1
                                                                              ATTTTAGAAGAGTCCTCGCTCTTGAATGCAGATGAGCTGCTTGAGCTTTTTT‑‑TAGATTCAGAGGAGGTGCTTGTTGCTGTAACAGATTGCGTTTGCTGCTGGCTGGATGATTGTTTTGCAGTGTCCTGAGAATCATGGTC  <  1:4251800/140‑1
                                                                                                                                                                                                                            
AACAAAGAGCTGTTTGACTGAGCCTGAAGATTTGTAATGCCGACAATCGCCGGTGAAATGTCTTCCACCATATCTGAGATTTTAGAAGCATCCTCGCTCTTGAATGCAGATGAGCTGCTTGAGCTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP034943/1328714‑1328839
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TTTTGCTAGATTCAGAGGACGTGCTTTTAGCTGTAACCGATTGTGTTTGCTGCTGGCTGGATGATTGTTTTGAAGTGTCCTGAGCATCATGGTC  >  NZ_CP034943/1328847‑1328940

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.