New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4166549 =NA (NA)3 (0.070) 3/592 6.5 4.9% intergenic (+264/‑110) murI/rrsB glutamate racemase/16S ribosomal RNA of rrnB operon
?NC_000913 4240228 = 59 (1.280)intergenic (‑173/‑97) yjbT/psiE putative periplasmic protein/phosphate starvation inducible protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below cutoff threshold.

TTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACG  >  NC_000913/4166250‑4166584
                                                                                                                                                                                                                                                                                                           |                                   
ttCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGGTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGc                                    <  1:50171/301‑1 (MQ=255)
                aaCTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGTCGTATTATGCACACCCCGCGCCGCTGAGTAAAAGCGAAGCGTGACTTATCTTTAACAAATGATCAGACACGCGGTGTGGTGTCTTGAAGACACGGAt                    >  1:341326/1‑301 (MQ=14)
                  cgtgCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCTGTCATAAAATTATATTAAATATCCTCGTATCATGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGAGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACAACCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTc                  <  1:302402/298‑1 (MQ=255)
                      cAGTTTTAGGTTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTCACAATTTATCAGACAAGCTGTGTGGGCACTCGACGATACGGATTCTTaa              >  1:151934/1‑301 (MQ=255)
                            tAGGCTGAGTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATGTCCTCTTGTCAGGCCGGAATAACTCCCGATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGc        <  1:408550/301‑1 (MQ=25)
                                cTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAg     <  2:375083/300‑1 (MQ=35)
                                   aTTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACg  <  2:149207/300‑1 (MQ=34)
                                                                                                                                                                                                                                                                                                           |                                   
TTCGAAACGCTCGAAAAACTGGCAGTTTTAGGCTGATTTGGTTGAATGTTGCGCGGTCAGAAAATTATTTTAAATTTCCTCTTGTCAGGCCGGAATAACTCCCTATAATGCGCCACCACTGACACGGAACAACGGCAAACACGCCGCCGGGTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAGGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCGAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACG  >  NC_000913/4166250‑4166584

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