New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	577820 =	24 (0.940)	3 (0.120)	3/484	2.8	11.1%	coding (208/498 nt)	rrrD	lysozyme
	?	NC_000913	= 1469240	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region
Rejected: Frequency below/above cutoff threshold.

ATTAAAGAAGATAAGCGCAAGGCTGCGAGAGGTGAATAATGCCTCCATCATTACGAAAAGCCGTTGCTGCTGCTATTGGTGGCGGAGCAATTGCTATAGCATCAGTGTTAATCACTGGCCCAAGTGGTAACGATGGTCTGGAAGGTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCGCGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCTGGCAATTTTAGAACATC‑GACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCT  >  NC_000913/577575‑578058                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                aTTAAAGAATATAAGCGCAAGGCTGCGAGAGGTGAATAATGCCTCCATCATTACGAAAAGCCGTTGCTGCTGCTATTGGTGGCGGAGCAATTGCTATAGCATCAGTGTTAATCACTGGCCCAAGTGGTAACGATGGTCTGGAAGGTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGTTAAAACGTATACCAAAGCAGAATGc                                                                                                                                                                                                                                            <  2:284610/251‑1 (MQ=255)    agagaagaTAAACGAGAGGCTGCGAGAGGAGAAGAATGCCTCCCTCATTACGAAAAGCCGTTGCTGCTGCTATTGGAGGCGGAGAAATTGCTATAGCAACAGTGTTAAGCACTGGCCCAAGTGGTAACGAAGGTATGGCAGGTGTCAGCTACATAACATACAAGGATATGGTTGGTGGATGGGCTGGATGTCACGGGCACACCCGAAAAGACATCATGCTCGGTACAGCGGTTACCAAAACAGAATGCaa                                                                                                                                                                                                                                          >  2:36584/3‑250 (MQ=255)    aaaGAAGATAAGCGCAAGGCTGCGAGAGCTGAATAATGCCTCCATCATTACGAAAAGCCGTTGCTGCTGCTAGTGGTTGCGGATCAATTGCTATAGCATCAGTGTGAAGCACATGCCCAAGGGGTAACGATGGTCTGGAAGTTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCAGGCTCGGTAATACGTATACCAAAGGAGAATGCaaa                                                                                                                                                                                                                                         <  2:144693/251‑1 (MQ=255)                 gcAAGGCTGCGAGAGGTGAATAATGCCTCCATCATTACGAAAAGCCGTTGCTGCTGCTATTGGTGGCGGAGCAATTGCTATAGCATCAGTGTTAATCACTGGCCCAAGTGGTAACGATGGTCTGGAAGGTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATa                                                                                                                                                                                                                            <  1:11375/251‑1 (MQ=255)                         gCGAGAGGTGCATAATGCATCCATCATTACGAATAGCCGTTGCTGCTGCTATTGGTGGCGGAGCAATTGCTATAGAATCAGTGTTAATCACTGGCCCAAGTGGTACGGATGGTCTGGAATGGGTCAGCGACATACCAAAAAAGGATATTGTTGGTGTATGGACTGTATACCACGGCCACACCGGAAAAGACATCATGCTCGGTAAAACAGATACTAAAGAAGAATGCAAAGCACTCTTGCATAAAgcactt                                                                                                                                                                                                                    >  2:191640/1‑246 (MQ=255)                               ggTGAATAATGCCTCCATCATTACGAAAAGCCGTTGCTGCTGCTATTGGTGGCGGAGCAATTGCTATAGCATCAGTGTTAATCACTGGCCCAAGTGGTAACGATGGTCTGGAAGGTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACt                                                                                                                                                                                                              >  1:44014/1‑251 (MQ=255)                                   aataatGTCTCCATCATTACGAAAAGCCGTTGCTGCTGCTATTGGTGGCGGAGCAATTGCTATAGCATCAGTGTTAATCACTGGCCCAAGTGGTAACGATGGTCTGGAAGGTGTCAGCTACATACCATACAAAGATATTGTTTGTGTATGGACTGTATGTCACGGACACACCGCAAAAGACATCTTGCTCGGCAAAACGTATCCCAAAGCAGAATGCAAATCCCTCTTGAATAAAGACCTTGCCACTGTCg                                                                                                                                                                                                          <  2:80252/251‑1 (MQ=255)                                   aataatGCCTCCATCATTACGAAAAGCCGTTGCTGCTGCTATTGGTGGCGGAGCAATTGCTATAGCATCAGTGTTAATCACTGGCCCAAGTGGTAACGATGGTCTGGAAGGTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGACTAAAGACCTTGCCACTGTCg                                                                                                                                                                                                          >  1:15107/1‑251 (MQ=255)                                           cTCCATCATTACGCAAAGCCGTTGCTGCTGCTATTGGTGGCGGAGCAATTGCTATAGCATCAGTGTTAATCACTGGCCCAGGTGGTAACGATGGTCTGGAAGGTGTCAGATACATACCATACAAAGATATTGTTGGTGGATGGACTGTATGTCACGGACACACCGGAAAAGAAATCATGCTCGGTAAAACGTATACCAACGCAGAATGCAAAGCACTCTTGAATAAAGACCCTGCCACTGTCGCCAGaaa                                                                                                                                                                                                   >  2:19076/1‑248 (MQ=255)                                                             ccGTTGCTGCTGATATTGGTGGCGGAGCACTTGCCATAGCATCAGTGTTCATCAATGGCCCAAGTGGTAACGATGGTCTGGAAGGTGTCAGCTACATACCATACAAAGGTATTGTTGGTGTATGGCCTGCATGTCACGGACAAACCGGAAAAGACATCATGCTCGGAAAAAAGTATACAAACGCAGAATGAAAAGCACTCTTGAATAAAGACCTTGCAACTGCCGCCAGACAAATTACCCCGTATATCaa                                                                                                                                                                                 >  2:97423/1‑250 (MQ=255)                                                                        gctATTGGTGGCGGAGCAATTGCTATAGCATCAGTGTTAATCACTGGCCCAAGTGGTAACGATGGTCTGGAAGGTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATAc                                                                                                                                                                      <  1:5402/250‑1 (MQ=255)                                                                                    ggAGCAATTGCTATAGCATCAGTGTTAATCACTGGCCCAAGTGGTAACGATAGTCTGGAAGGTGTCAGATACATACCATACAACGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAAcgcg                                                                                                                                                         >  1:111358/1‑251 (MQ=255)                                                                                        cAATTGCTATAGCATCAGTGTTAATCACTGGCCCAAGTGGTAACGATGGTCTGGAAGGTGGCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCgcggc                                                                                                                                                     <  2:71746/251‑1 (MQ=255)                                                                                                                  ccTGGCCCAAGTGGTAACGATGGTCTGGAAGGTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCGCGGCGCTCTTTACGCATTCGTTTACAACGt                                                                                                                           <  1:254868/250‑1 (MQ=255)                                                                                                                         cAAGTGGTAACGATGGTCTGGAAGGTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCGCGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCt                                                                                                                    >  1:1796/1‑251 (MQ=255)                                                                                                                                tttctccggccctaagacatcaatcatctgttctccaatgactagtctaaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCGCGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCTGGCAAtt                                                                                                             <  1:83246‑M1/133‑1 (MQ=255)                                                                                                                                    gATGGTCTGGAAGGTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGAAATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCGAGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCTGGCAAtttt                                                                                                           >  2:241609/1‑249 (MQ=255)                                                                                                                                        gccctaagacatcaatcatctgttctccaatgactagtataaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCGCGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCTGGCAATTTTAGAAc                                                                                                      <  2:62114‑M1/140‑1 (MQ=255)                                                                                                                                                     cAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCGCGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCTGGCAATTTTAGAACATC‑GACGCttcttc                                                                                       >  1:47832/1‑251 (MQ=255)                                                                                                                                                                              ataaaaactagtattaagactatcacttatttaagtgatattggttgtctggagattcagggggccagtctaGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCGCGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCTGGCAATTTTAGAACATC‑GACGCTTCTTCGCAAAATAAACCAGGGCGATATCa                                                               >  1:62114‑M1/73‑250 (MQ=255)                                                                                                                                                                                    gtATGGACTGCATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAACCGTATACCAAAGCAGAATGCAAAGCACTCTTGACTAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGCAACAACGCGCGGCGCTCTTTACTCATTCGTTTACAACGTGGGCGCTGGCAATTTTAGAACATC‑GACGCTTCTTCGCAAAATAAACCAGGGCGATCTCaacgacgc                                                        >  2:103689/1‑245 (MQ=255)                                                                                                                                                                                                         cacacCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCGCGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCTGGCAATTTTAGAACATC‑GACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTg                                   >  1:218484/1‑251 (MQ=255)                                                                                                                                                                                                                                       aCGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTTACCCGTATATCAAAGTCGATATACCGGAAACAACGCGCGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCTGGCAATTTTAGAACATC‑GACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAgg     <  1:10742/251‑1 (MQ=255)                                                                                                                                                                                                                                       aCGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCGCGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCTGGCAATTTTAGAACATC‑GACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAgg     >  1:186321/1‑251 (MQ=255)                                                                                                                                                                                                                                         gTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCGCGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCTGGCAATTTTAGAACATC‑GACGCTTCTTCGCAAAATAAACCAGGGCGATATCACAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGGAAg      >  2:99058/1‑248 (MQ=255)                                                                                                                                                                                                                                         gTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCGCGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCTGGCAATTTTAGAACATC‑GACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAg      <  1:99058/248‑1 (MQ=255)                                                                                                                                                                                                                                          tataCCAAAGCAGAATGCAAAGCACTCTGGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCACGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCGAGCAATTTAAGAACATCTGA‑GCTTCTTCGAAAGATAAACCAGGACGATATCAAAGGCGCAGGTGATCAGCTGCGTCTCTGGACATACTATAGCGGTAAGCAATGGAAAAGCCt  >  1:122047/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                ATTAAAGAAGATAAGCGCAAGGCTGCGAGAGGTGAATAATGCCTCCATCATTACGAAAAGCCGTTGCTGCTGCTATTGGTGGCGGAGCAATTGCTATAGCATCAGTGTTAATCACTGGCCCAAGTGGTAACGATGGTCTGGAAGGTGTCAGCTACATACCATACAAAGATATTGTTGGTGTATGGACTGTATGTCACGGACACACCGGAAAAGACATCATGCTCGGTAAAACGTATACCAAAGCAGAATGCAAAGCACTCTTGAATAAAGACCTTGCCACTGTCGCCAGACAAATTAACCCGTATATCAAAGTCGATATACCGGAAACAACGCGCGGCGCTCTTTACTCATTCGTTTACAACGTGGGTGCTGGCAATTTTAGAACATC‑GACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCT  >  NC_000913/577575‑578058

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