New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4169328 =	76 (0.560)	71 (0.600) +GGAATTAATT	50/130	1.1	51.9%	noncoding (113/120 nt)	rrfB	5S ribosomal RNA
	?	NC_000913	= 4169500	NA (NA)					intergenic (+165/‑136)	rrfB/murB	5S ribosomal RNA/UDP‑N‑acetylenolpyruvoylglucosamine reductase

GGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCT  >  NC_000913/4169472‑4169575                             |                                                                            gggCAGCACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCACCCCGACCGATGGCCTTTCTGCGtt                              >  2:1457615/1‑76 (MQ=11)   gCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTc                            >  2:145362/1‑76 (MQ=38)    cAGGACCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCt                           <  1:4537417/76‑1 (MQ=255)     aGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTa                          >  2:891445/1‑76 (MQ=40)      ggACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCtcc                         >  2:1939924/1‑74 (MQ=255)      ggACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTAc                         >  2:4362241/1‑76 (MQ=255)       gACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACa                        >  2:589041/1‑76 (MQ=255)         cccccgccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACaaa                      <  1:3758563/74‑1 (MQ=255)          ccccgccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAc                     <  1:4528122/75‑1 (MQ=255)          ccccgccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAc                     <  1:2032672/75‑1 (MQ=255)           cccgccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAct                    >  1:718046/1‑76 (MQ=255)           cccgccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAct                    <  2:2247105/76‑1 (MQ=255)           cccgccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAct                    >  2:1494704/1‑76 (MQ=255)            ccgccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAActc                   >  1:4269601/1‑76 (MQ=255)             cgccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAActct                  >  1:867753/1‑76 (MQ=255)             cgccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAActct                  >  1:3037840/1‑76 (MQ=255)              gccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCtt                 <  2:3453570/76‑1 (MQ=255)              gccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTCACGGACGGCCTTTGTGCGTTTCTACAAActctc                 >  2:4598477/1‑75 (MQ=17)               ccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTc                >  2:838799/1‑76 (MQ=255)               ccATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTc                <  1:1070/76‑1 (MQ=255)                cATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTcc               >  2:4696149/1‑76 (MQ=255)                 aTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCt              <  2:3502636/76‑1 (MQ=255)                 aTAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCt              <  2:3646875/76‑1 (MQ=255)                  tAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTg             <  2:867753/76‑1 (MQ=255)                    aaCTGCCCGGCATCAAATTACGCAGAAGGCCATCCTGACGGATGGTCTTTTTGCGTCTCTACCAACTCCTCCtctc           >  2:3133184/1‑76 (MQ=14)                    aaCTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgtc           >  2:3129204/1‑76 (MQ=255)                    aaCTGCCAGGCATCAAATTAAGCAGAAGGCCAGCCTGACGGATGGCATTTCTGCGGTTCTACAAACTCGTCCTgtc           >  1:4147382/1‑76 (MQ=17)                      cTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgtcgt         >  2:1765389/1‑76 (MQ=255)                      cTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgtcgt         >  2:3479444/1‑76 (MQ=255)                       tGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgtcgtc        >  1:4798603/1‑76 (MQ=255)                       tGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgtcgtc        <  1:3248545/76‑1 (MQ=255)                       tGCCAGGCATCAAATTAAGCAGAAGGCCATCATGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgtcgtc        >  2:4380660/1‑76 (MQ=255)                         ccAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCat      >  2:1771024/1‑76 (MQ=255)                           aGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCatat    >  2:2635740/1‑76 (MQ=255)                           aGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTGTGAGTTTCTACAAACTCTTCCTGTCGGCatat    >  1:4301074/1‑76 (MQ=21)                            ggCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATc   <  2:281970/76‑1 (MQ=255)                             gCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCt  <  1:3013548/76‑1 (MQ=255)                             gCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCt  <  1:1330390/76‑1 (MQ=255)                             |                                                                            GGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCT  >  NC_000913/4169472‑4169575

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