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.

GAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCTAACGACAAATCATACAGGCGGAAACG  >  NC_000913/1243841‑1244110
                                                                                                                                              |                                                                                                                               
gaAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGACTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAacac                                                                                                                           >  6:16856/1‑149 (MQ=255)
 aaaaTAATCCGCTGTAATTGCCGCTCCACCGTCGGACTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACAcc                                                                                                                          <  5:16856/149‑1 (MQ=255)
                                    gCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGccc                                                                                                                 <  7:60331/123‑1 (MQ=255)
                                    gCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGGTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGccc                                                                                                                 >  8:60331/1‑123 (MQ=255)
                                             aaGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAgc                                                                                       >  3:28979/1‑140 (MQ=255)
                                             aaGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAgc                                                                                       <  4:28979/140‑1 (MQ=255)
                                                       aaGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCAtt                                                                                                                     <  8:40272/100‑1 (MQ=255)
                                                       aaGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCAtt                                                                                                                     >  7:40272/1‑100 (MQ=255)
                                                          cTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTaa                                                                       <  8:61554/143‑1 (MQ=255)
                                                          cTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTaa                                                                       >  7:61554/1‑143 (MQ=255)
                                                                                tttCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGtt                                                                                                           >  8:4937/1‑85 (MQ=255)
                                                                                tttCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGtt                                                                                                           <  7:4937/85‑1 (MQ=255)
                                                                                                aatcaCTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCt                            >  1:154866/1‑148 (MQ=255)
                                                                                                            ctcgctGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTg                                                                          <  3:118216/90‑1 (MQ=255)
                                                                                                            ctcgctGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTg                                                                          >  4:118216/1‑90 (MQ=255)
                                                                                                               gctGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCAt                                                                                   <  1:31563/78‑1 (MQ=255)
                                                                                                               gctGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCAt                                                                                   >  2:31563/1‑78 (MQ=255)
                                                                                                                 tGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCTAACGACAAATCATACAgg          >  2:72515/1‑149 (MQ=255)
                                                                                                                     ggAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCTAACGACAAAt                  >  4:64114/1‑137 (MQ=255)
                                                                                                                     ggAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCTAACGACAAAt                  <  3:64114/137‑1 (MQ=255)
                                                                                                                      gagaTTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTc                                                  <  4:18441/104‑1 (MQ=255)
                                                                                                                      gagaTTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTc                                                  >  3:18441/1‑104 (MQ=255)
                                                                                                                          ttAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCTAACGACAAATCATACAGGCGGAaacg  >  8:117376/1‑148 (MQ=255)
                                                                                                                                        aaGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCTAACGACaa                    <  6:966/116‑1 (MQ=255)
                                                                                                                                        aaGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCTAACGACaa                    >  5:966/1‑116 (MQ=255)
                                                                                                                                          gTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGc                                                             >  6:37998/1‑73 (MQ=255)
                                                                                                                                          gTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGc                                                             <  5:37998/73‑1 (MQ=255)
                                                                                                                                              |                                                                                                                               
GAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGGAGATTAACTGGGCGTCAAAGTGAAAAACACCCCATTGCCCCATGTTGACTTCAATCTGAGCGAAGCGCATGCCTTCTTGTAATTCGGCAGGCTTTGCTGTTTCCAGTAATGCGCCCATGCCGCCTAACGACAAATCATACAGGCGGAAACG  >  NC_000913/1243841‑1244110

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