Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,908,726	IS2 (–) +5 bp	coding (37‑41/144 nt)	mgrB ←	PhoQ kinase inhibitor

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1651548	NA (NA)	33 (0.680)	30/494	0.6	100%	noncoding (706/706 nt)	IS2	repeat region
	?	NC_000913	= 1908730	0 (0.000)					coding (37/144 nt)	mgrB	PhoQ kinase inhibitor
*	?	NC_000913	1908726 =	0 (0.000)	40 (0.820)	35/494	0.4	100%	coding (41/144 nt)	mgrB	PhoQ kinase inhibitor
	?	NC_000913	= 2070271	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

GCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTG  >  NC_000913/1651303‑1651776                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      gCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAAGCCAC‑AGGt                                                                                                                                                                                                                                 <  1:160152/251‑1 (MQ=14)                 aaCGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCtg                                                                                                                                                                                                                 <  2:93484/251‑1 (MQ=38)                 aaCGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCtg                                                                                                                                                                                                                 <  2:171773/251‑1 (MQ=38)                  aCGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCATGATCAAGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCGTGCAGAGGCTGTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAATTCCAC‑AGGTGAACTACGCTACTCtgt                                                                                                                                                                                                                <  2:476742/251‑1 (MQ=11)                  aCGGAATAGCAGAGAACTTCGTGAAAACGATAAAGCGTGACTACATCAGTAACATGCCCAAACCAGACGGGTTAACTGCAGCAAAAAACCTTGCAGAGGCGTTCGAGCATTATAACGAATCGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGCATATCTGCGGCAGCGGGCTTGTCATGGGTTAAGTGATACCAGATGTCTGGAACTATAGGGGCAACTC‑ACGCGGTGAACTACTCTCCTCtgt                                                                                                                                                                                                                >  2:374512/1‑251 (MQ=11)                    ggAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCtgttt                                                                                                                                                                                                              <  1:381445/251‑1 (MQ=40)                         aGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAAGGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑ATGTGAACTACGCTCCTCTGTTTGTTAc                                                                                                                                                                                                         <  1:374512/251‑1 (MQ=255)                         aGCAGAGAGCTGCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGTGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTa                                                                                                                                                                                                          <  1:16042/250‑1 (MQ=255)                           cAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCGGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGTGCTTGTTATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGGGAACTACGTTCCGCTGGTTGTTACGc                                                                                                                                                                                                       <  2:228944/251‑1 (MQ=11)                            agagagCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAAACTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGc                                                                                                                                                                                                       >  1:494529/1‑250 (MQ=255)                                      gTGAAAACGATAAAGCGTGACTACATCAGTCTCATGCCCAAACCAGACGGGATAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGACTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTACTGGGTTAAGTGATAACCGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTg                                                                                                                                                                                            >  1:440495/1‑251 (MQ=255)                                              gATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTc                                                                                                                                                                                    >  1:92325/1‑251 (MQ=255)                                                        gACTACATCAGTATCATGCCCAAACCAGCCGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCa                                                                                                                                                                           >  2:295262/1‑250 (MQ=255)                                                            aCATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACa                                                                                                                                                                      <  2:169417/251‑1 (MQ=255)                                                                cAGTATCATGCCCAAACCAGACGTGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGttt                                                                                                                                                                  <  1:115693/251‑1 (MQ=255)                                                                             aaaCCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAACTATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGc                                                                                                                                                      >  1:24909/1‑250 (MQ=255)                                                                                ccAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATc                                                                                                                                                  <  1:396316/251‑1 (MQ=255)                                                                                                 agcaAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTg                                                                                                                                   >  2:398277/1‑249 (MQ=255)                                                                                                   caAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTGCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGc                                                                                                                               <  2:15369/251‑1 (MQ=255)                                                                                                              gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGTGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGAtt                                                                                                                    <  1:377193/251‑1 (MQ=255)                                                                                                                      gTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATa                                                                                                            <  1:241946/251‑1 (MQ=255)                                                                                                                                       aatggaaccagCATACTGAGCTGGGATATCGCACGCCACAGGAATATATGCGGGAACGCGCTAGTAATGGGTTCAATGATACCAGATGTCTGGACCTATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTATGCTACGAAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGAAATCTGGAAAAATTTGGTGAGGCCCGGCAACATTTAACCATAAAAGAACTGGCTGaaa                                                                                            >  2:457712/11‑250 (MQ=255)                                                                                                                                                   aTAGTGCGCTGGGATATAGCTCGAAACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAATCCAGTTTCCCTCCTCACGCAATCTGGTAAAATTTGGTAAAGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAACAAATTCTTGCACt                                                                               >  2:365250/1‑251 (MQ=37)                                                                                                                                                      gTGCGCTGGATTATCGCTCGCACCGGGACTATCTGCGGCAGCGGGCTTGTAGTGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑CGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGACAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTga                                                                             >  2:350480/1‑250 (MQ=38)                                                                                                                                                            tGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGAGGGCTTGTAATGGGTTACGTGATAACAGATGTCTGAAAATATAGGGGCAAATCAAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAACAGAACTGGCTGAAAAATTTATGGCACTGAAAGCAAc                                                                      >  1:70233/1‑251 (MQ=38)                                                                                                                                                               gTTATCGCTAGCCACGGGAATAGCTGCGGCAGCGGGCTTGTAATGGGTTAAGAGAGAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGTTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGAAAAACAGTTTCCCTCCTCAAGCAATCTGGAAAAATTTGGCGAGGCCTGACAAGAATTAACAATAAAAGAAATGGCTGAAAAATTTCTgttactgatcaaaactgt                                                                   >  2:286063/1‑234 (MQ=25)                                                                                                                                                                     gctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGc                                                             >  2:284505/1‑251 (MQ=255)                                                                                                                                                                                  atatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTc                                                >  1:28415/1‑251 (MQ=255)                                                                                                                                                                                            ccgCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGTCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATAcc                                      <  2:264318/249‑1 (MQ=255)                                                                                                                                                                                                      gTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCGAC‑AGCTGAACTACGCTCCTATGTTTGTTACGAAATAAACACTGGTGTTTTCGACTATGCAAAACAGTTTCACTCATAACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGtcat                            >  1:323416/1‑251 (MQ=255)                                                                                                                                                                                                                                cTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCACTAAAAACTGGCGGTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAACTTTGGTGAGGCCCGACACGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGACGCCAAAACATCACTCAACACATACCGTGCCGACATCAAAAATATCCTGAGCATAATCGGTg  >  2:396874/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      GCGGAGTCCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAC‑AGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTG  >  NC_000913/1651303‑1651776

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.