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)	38 (0.820)	36/488	0.3	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)	27 (0.580)	24/488	0.8	100%	noncoding (1/1258 nt)	IS3	repeat region
	?	NC_000913	= 577690	0 (0.000)					coding (78/498 nt)	rrrD	lysozyme

AAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGGTG  >  NC_000913/316246‑316730                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                        aaaaGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGCAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCCGTCCGGAACCATTTAAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGGAGGATCATAAAGTCCTCggg                                                                                                                                                                                                                                            >  1:143591/1‑251 (MQ=255) aaaaGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGCACAATTTGAAAAACAGAACCTCGCTTAGGCCTGTGTCCCTATTACGTGGGTAGGATCATAAAATCCTCggg                                                                                                                                                                                                                                            >  2:158020/1‑251 (MQ=255)      gTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTTTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGTTCGtt                                                                                                                                                                                                                                       <  1:70929/251‑1 (MQ=25)                  aTAATGCCTGCGTGGAAAGCTTCTTTCAGTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCtttt                                                                                                                                                                                                                           <  2:165326/251‑1 (MQ=255)                                            cATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATGCACATACCg                                                                                                                                                                                                  <  1:157764/250‑1 (MQ=255)                                                    gAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGACAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGACCATCAAGGCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATAATCACATCCCGTAACAGCa                                                                                                                                                                                          >  2:404089/1‑250 (MQ=17)                                                                       ggAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGGTGGGGTGGCGGCCTCCGGCAGGAAAAAATTGAAAACAAGAACCTCGCTTAGGCCTGCGTCCATATTACGGGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCGACCGTGAGCa                                                                                                                                                                       >  1:302312/1‑250 (MQ=21)                                                                          gaACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAGTCGGTGGCGGCGGCTCAGTTGGTGTGGCGGTCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCTTCATAGTCTCATACCGTTCCAGCAAAAGTTTCAGCCGTGAGCAAATc                                                                                                                                                                   <  1:78109/251‑1 (MQ=18)                                                                          gaACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATc                                                                                                                                                                   <  1:109337/251‑1 (MQ=255)                                                                                   aTCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGGGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGTGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGAttt                                                                                                                                                          <  2:357166/251‑1 (MQ=255)                                                                                            gAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATaaa                                                                                                                                                                                                                                                     >  2:112223/1‑151 (MQ=32)                                                                                            gAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATaaa                                                                                                                                                                                                                                                     <  1:112223/151‑1 (MQ=32)                                                                                              aataatGCGGGCAACGGTGTTTAATTATAGCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGtt                                                                                                                                               <  2:465625/251‑1 (MQ=255)                                                                                               ataatGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTc                                                                                                                                              <  1:365340/251‑1 (MQ=255)                                                                                                                   taattagtgcgaatggtaTTACAGTCTGTGGCGGCGACACAGTTGGTGTGGCGGACTCAGTCCGGATCAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCGTTCTTTTGGCATGTCAGCATATTCACATAGCGTACCAGCAAAAGGTTCAACCGTGAGCAAGTCGACCGATTATGTTGACGGTTCAAACGTGGAAACGAAATAg                                                                                                                           <  2:305079/233‑1 (MQ=25)                                                                                                                       tataTCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCGCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGatat                                                                                                                       <  2:466553/250‑1 (MQ=255)                                                                                                                            cGAATGTGATTACAATCGGTGCCGGCGGCACAGTTGGTGTGGCGGGCTCAGTCCGGAACCATTTGAACACAAAAACCTCGCTTAGGCCTGTGTCCAATTTACGTGGGTAGGAGCATCAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGGTTCCACCGTGAGCAAATCGACCGAGTTTGTTGCCGGATCAAAAGTGGGAACGGAATAGatttcaggct                                                                                                                 >  2:41709/1‑243 (MQ=255)                                                                                                                              aaTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCtt                                                                                                                >  1:372835/1‑250 (MQ=255)                                                                                                                               aTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTa                                                                                                               >  1:388867/1‑250 (MQ=255)                                                                                                                                 gtgAATACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCTACCGATTTTGTTGACGGTTCAAACGTGGAAAGGAAATAGATATCCGTCTTaaaa                                                                                                            <  2:434537/251‑1 (MQ=255)                                                                                                                                      tgcAATCGGTGGCGGCGGCACAGGTGTTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATAATACGTGGGTATGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCAGATTCACATACCGTACCAGCAAAAGTTTCAACCGTGGGCAAATCGACCGATTTTGGTGACGGTTCAAACGTGGAAACGAAATAGAGATCCGGCTTAAAAttct                                                                                                        <  2:411340/248‑1 (MQ=37)                                                                                                                                       aCAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAAttcttc                                                                                                      <  2:128899/251‑1 (MQ=255)                                                                                                                                                   ggcggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCt                                                                                          <  1:165682/251‑1 (MQ=255)                                                                                                                                                       gcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGc                                                                                                                                                                        >  1:94552/1‑169 (MQ=255)                                                                                                                                                       gcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGc                                                                                                                                                                        <  2:94552/169‑1 (MQ=255)                                                                                                                                                                       ggcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGttt                                                                      >  2:147058/1‑251 (MQ=255)                                                                                                                                                                            cAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGActt                                                                 >  1:92601/1‑251 (MQ=255)                                                                                                                                                                                                               aGGCCTGTGTCCATATTACGGTGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACAAACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATTCATCGGGCTAAATGATGACGTCGTAGTTTGCCGTCTTCATCCGGTGATTTATTCGACCCTTCAATGTCa                              <  2:302312/251‑1 (MQ=255)                                                                                                                                                                                                                 gCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCata                            <  1:404089/251‑1 (MQ=255)                                                                                                                                                                                                                         ccATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATg                    <  2:388867/251‑1 (MQ=255)                                                                                                                                                                                                                                           aTCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGgtg  >  2:418367/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                        AAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGACGGTTCAAACGTGGAAACGAAATAGATATCCGTCTTAAAATTCTTCATACATCTGGCTAAATGATGACGTCGTAGTTTGACTTCTTCATCCGGTGATTTATTCGACACTTCAATGTCATATTTAAATGAGTTGCCAAAAATAGGTG  >  NC_000913/316246‑316730

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