New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 394682420 (0.230)4 (0.050) 4/278 6.9 6.2% intergenic (+5/‑48) rrfC/aspT 5S ribosomal RNA of rrnC operon/tRNA‑Asp
?NC_000913 = 4171801 102 (1.200)intergenic (+45/‑256) rrfB/murB 5S ribosomal RNA of rrnB operon/UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAG  >  NC_000913/3946677‑3946966
                                                                                                                                                   |                                                                                                                                              
aaGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAAt                                                                                                                                                >  2:322330/1‑148 (MQ=32)
     gTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAg                                                                                                                                            >  7:51428/1‑147 (MQ=32)
     gTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCa                                                                                                                                          <  4:264502/149‑1 (MQ=32)
      tCTGATGAATCAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAg                                                                                                                                         <  2:164341/149‑1 (MQ=17)
       cTGATGAAACAGAATTTGCCTTGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTTTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGt                                                                                                                                        <  6:291462/149‑1 (MQ=17)
         gatgaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGt                                                                                                                                        >  7:108354/1‑147 (MQ=32)
          atgaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTa                                                                                                                                       <  5:220281/147‑1 (MQ=32)
           tgaAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAGc                                                                                                                                    >  7:330333/1‑149 (MQ=31)
                                                                                                                             tAGGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCAcc                  >  8:5399/1‑149 (MQ=11)
                                                                                                                                    cTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTa           <  7:171681/149‑1 (MQ=35)
                                                                                                                                       ccAGGCATCAAATTAAGCAGTAAGCCGGGCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGGCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAggg        <  3:245616/149‑1 (MQ=255)
                                                                                                                                       ccAGGCAGCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGTTGGAGCGGTAGTTCATTCTGTTAGAATAACTGCCTGTCACGCAGGTGTTCGCGTGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAggg        <  1:95710/149‑1 (MQ=11)
                                                                                                                                         aGGCATCAAATTAAGCAGTAAGCCGGTCATAAACCCGGTGCTTGTAAAAGAAATCGGTGGAGCGGAAGTTCAGTCGGTAAGAATACCTGCCAGTCACCCACCGCGTCGCGGGTTCGACTCCCGTCCCATCCGCCACCCTAAATAGCGGc      >  8:173541/1‑149 (MQ=11)
                                                                                                                                         aGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGc      >  8:6129/1‑149 (MQ=40)
                                                                                                                                           gCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCg     >  7:125773/1‑148 (MQ=40)
                                                                                                                                           gCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGt    >  5:121695/1‑149 (MQ=255)
                                                                                                                                            cATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTa   <  2:920/149‑1 (MQ=255)
                                                                                                                                            cATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTa   <  8:330333/149‑1 (MQ=255)
                                                                                                                                             aTCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCTCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGTGGAGTAg  >  4:280987/1‑149 (MQ=255)
                                                                                                                                             aTCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCCCCCTAATTAGGGGCGTAg  >  2:67770/1‑149 (MQ=255)
                                                                                                                                             aTCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAg  >  3:139760/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAAGCAGTAAGCCGGTCATAAAACCGGTGGTTGTAAAAGAATTCGGTGGAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAG  >  NC_000913/3946677‑3946966

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.