New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1243939 =17 (0.610)4 (0.150) 4/224 3.1 17.1% coding (589/735 nt) ycgR flagellar velocity braking protein, c‑di‑GMP‑regulated
?NC_000913 1243983 = 23 (0.890)coding (545/735 nt) ycgR flagellar velocity braking protein, c‑di‑GMP‑regulated
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTTTTCACTTTGACGCCCAGTTAATCTCCATCAGCGAGCGCAAAGTGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1243986‑1243939
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caaagtgaAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCTAACGACAA  >  NC_000913/1243983‑1244092
                                                                                                                                                              
TTTTTCACTTTGACGCCCAGTTAATCTCCATCAGCGAGCGCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGC                                               <  5:86491/113‑1
TTTTTCACTTTGACGCCCAGTTAATCTCCATCAGCGAGCGCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGC                                               >  6:86491/1‑113
           GACGCCCAGTTAATCTCCATCAGCGAGCGCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAG                             <  5:30639/120‑1
           GACGCCCAGTTAATCTCCATCAGCGAGCGCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAG                             >  6:30639/1‑120
                                          AAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCTAACGACAA  >  5:966/1‑116
                                          AAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCTAACGACAA  <  6:966/116‑1
                                            GTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGC                                           <  5:37998/73‑1
                                            GTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGC                                           >  6:37998/1‑73
                                                                                                                                                              
TTTTTCACTTTGACGCCCAGTTAATCTCCATCAGCGAGCGCAAAGTGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1243986‑1243939
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caaagtgaAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCTAACGACAA  >  NC_000913/1243983‑1244092

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