Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	577,688	IS3 (+) +3 bp	coding (76‑78/498 nt)	rrrD →	lysozyme

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 316483	NA (NA)	37 (1.090)	36/486	0.1	100%	noncoding (1255/1255 nt)	IS3	repeat region
	?	NC_000913	577688 =	0 (0.000)					coding (76/498 nt)	rrrD	lysozyme
*	?	NC_000913	566777 =	NA (NA)	14 (0.410)	13/486	1.2	100%	noncoding (1/1258 nt)	IS3	repeat region
	?	NC_000913	= 577690	0 (0.000)					coding (78/498 nt)	rrrD	lysozyme

GGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGGTG  >  NC_000913/316250‑316730                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        ggTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTAGCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGCCCGGAACAATGTGAAAACAAGAACCTCGCTTATGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGt                                                                                                                                                                                                                                        <  1:32756/251‑1 (MQ=255)  gTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGtt                                                                                                                                                                                                                                       <  2:260178/251‑1 (MQ=38)        gctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTGTAATTAGATCGAATGTGATTACAATCGGTGGCGGCGGCACAGGGGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGTATCATAAAGTCCTCGGGTCGTTGTCCa                                                                                                                                                                                                                                  <  2:345677/250‑1 (MQ=255)           aCGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCActct                                                                                                                                                                                                                              <  1:196330/251‑1 (MQ=255)                       gCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTTCGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCa                                                                                                                                                                                                                  <  1:246138/251‑1 (MQ=255)                             aaaGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTGTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAAGCGGTGGCGGCGGCACAGTTGGTGTGGCGGTTGCTGTCCGGAACAATTTGAAAACAAGAACCTCGCTGAGGCCTTTGTCCATAGTACGTGGTGAGGATCAGAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCatat                                                                                                                                                                                                            <  1:331366/251‑1 (MQ=17)                               aGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTc                                                                                                                                                                                                          <  2:91440/251‑1 (MQ=255)                                             gCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAg                                                                                                                                                                                            <  1:289133/251‑1 (MQ=255)                                             gCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAg                                                                                                                                                                                            <  2:226446/251‑1 (MQ=255)                                                    gTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGt                                                                                                                                                                                     <  1:146499/251‑1 (MQ=255)                                                                                     cGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTttgttg                                                                                                                                                     >  1:114095/1‑250 (MQ=255)                                                                                             aatGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGGCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCa                                                                                                                                             <  1:229519/250‑1 (MQ=255)                                                                                                                          aaTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGAGTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTa                                                                                                               <  2:118568/251‑1 (MQ=255)                                                                                                                                      aTCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGAcc                                                                                                                                                               <  1:79311/191‑1 (MQ=255)                                                                                                                                      aTCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGAcc                                                                                                                                                               >  2:79311/1‑191 (MQ=255)                                                                                                                                      aTCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGCTTTTGTTGACGGGTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCAt                                                                                                    >  1:190667/1‑250 (MQ=255)                                                                                                                                                                ggcggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCaaa                                                                                                                                                                     >  1:5842/1‑159 (MQ=255)                                                                                                                                                                ggcggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCaaa                                                                                                                                                                     <  2:5842/159‑1 (MQ=255)                                                                                                                                                                                          aaaCACGAAACTCTCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAATTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATGCACATACCGTACCAGCAAAAGTTGCATCCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAGACGAAATAGATCTCCGTCGTAAAATTCTTTATACATCTGGCTAAAGGATGACGTCGGAGTTTGCCTTCTGCATCCGCTGATTTAt                                               <  2:346341/251‑1 (MQ=255)                                                                                                                                                                                                                                       aTCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATTACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTGCAATGTCATATGTAAATGAGTTGCCAAAAATAGgtg  <  2:40637/251‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        GGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGGTG  >  NC_000913/316250‑316730

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