New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 390790526 (0.700)4 (0.120) 4/238 4.3 14.4% coding (462/774 nt) pstB phosphate ABC transporter ATPase
?NC_000913 = 3907917 23 (0.660)coding (450/774 nt) pstB phosphate ABC transporter ATPase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3907822‑3907905
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caccagagCGGTTACTCTCTCTCTGGTGGTCAGCAACAGCGTCTGTGTATTGCGCGTGGTATCGCCATTCGCCCGGAAGTGCTGCTGCTCGACGAACCGTGTTCGGCGCTCGACCCTATCTCTACCGGGCGTATTGAAGAGCTGATCAC  <  NC_000913/3907917‑3907777
                                                                                                                                                                                                                                 
GAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGCGGTTACTCTCTCTCTGGTGGTCAGCAACA                                                                                                                 <  1:152244/114‑1
GAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGCGGTTACTCTCTCTCTGGTGGTCAGCAACA                                                                                                                 >  2:152244/1‑114
              GAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGCGGTTACTCTCTCTCTGGT                                                                                                                            <  3:55029/89‑1
              GAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAGCGGTTACTCTCTCTCTGGT                                                                                                                            >  4:55029/1‑89
                                                                            CACCAGAGCGGTTACTCTCTCTCTGGTGGTCAGCAACAGCGTCTGTGTATTGCGCGTGGTATCGCCATTCGCCCGGAAGTGCTGCTGCTCGACGAACCGTGTTCGGCGCTCGACCCTATCTCTACCGGGCGTATTGAAGAGCTGATCAC  <  2:132416/149‑1
                                                                               CAGAGCGGTTACTCTCTCTCTGGTGGTCAGCAACAGCGTCTGTGTATTGCGCGTGGTATCGCCATTCGCCCGGAAGTGCTGCTGCTCGACGAACCGTGTTCGGCGCTCGACCCTATCTCTACCGGGCGTAT                 >  3:93702/1‑131
                                                                               CAGAGCGGTTACTCTCTCTCTGGTGGTCAGCAACAGCGTCTGTGTATTGCGCGTGGTATCGCCATTCGCCCGGAAGTGCTGCTGCTCGACGAACCGTGTTCGGCGCTCGACCCTATCTCTACCGGGCGTAT                 <  4:93702/131‑1
                                                                                                                                                                                                                                 
GAACACGGTTCGTCGAGCAGCAGCACTTCCGGGCGAATGGCGATACCACGCGCAATACACAGACGCTGTTGCTGACCACCAGAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3907822‑3907905
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caccagagCGGTTACTCTCTCTCTGGTGGTCAGCAACAGCGTCTGTGTATTGCGCGTGGTATCGCCATTCGCCCGGAAGTGCTGCTGCTCGACGAACCGTGTTCGGCGCTCGACCCTATCTCTACCGGGCGTATTGAAGAGCTGATCAC  <  NC_000913/3907917‑3907777

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