Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,792,772	IS2 (–) +5 bp	intergenic (‑41/+33)	nlpC ← / ← btuD	NlpC/P60 family lipoprotein NlpC/vitamin B12 ABC transporter ATP binding subunit

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1469240	NA (NA)	32 (0.790)	30/486	0.4	100%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	1792772 =	0 (0.000)					intergenic (‑41/+37)	nlpC/btuD	NlpC/P60 family lipoprotein NlpC/vitamin B12 ABC transporter ATP binding subunit
*	?	NC_000913	= 1651548	NA (NA)	33 (0.810)	32/490	0.4	100%	noncoding (706/706 nt)	IS2	repeat region
	?	NC_000913	= 1792776	0 (0.000)					intergenic (‑45/+33)	nlpC/btuD	NlpC/P60 family lipoprotein NlpC/vitamin B12 ABC transporter ATP binding subunit

GTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGA  >  NC_000913/1651317‑1651777                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      gTAACGGAATATCAGAGAGCTTCGTGAAAAAGATAAAGCGTGACTAAATCAGTCTCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAGTGGCATCCGCATAGAGAGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAAGGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGctcctc                                                                                                                                                                                                                    >  2:359856/1‑251 (MQ=11)        aaTAGCAGAGAGCTTCGTGAAAACGATAACGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGAGGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCtgtttgt                                                                                                                                                                                                             >  2:123647/1‑251 (MQ=255)                 gagattCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAAGCAGAAGGGGTAACGGGAGCAAAGAACCTGGCAGGGGCGTTGGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATa                                                                                                                                                                                                    <  1:88772/247‑1 (MQ=14)                      tCGTGAAAACGATAAAGCGTGACTACATCAGTATCCTGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCTCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGTGTTAAGTGATATCAGATGTCTGTAAATATAGTGGCAAATCCACAGGTGAACTACTCTCCTCTGTTTGTCACTCACTaaaaa                                                                                                                                                                                                <  2:144295/250‑1 (MQ=14)                               cGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGtttt                                                                                                                                                                                      <  1:368952/251‑1 (MQ=255)                                gATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTc                                                                                                                                                                                     <  1:270176/251‑1 (MQ=255)                                 aTAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGGTTGTTACGCAATAAAAACTGGTGTTTTCg                                                                                                                                                                                    <  1:123647/251‑1 (MQ=255)                                   aaaGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGAc                                                                                                                                                                                  >  1:213140/1‑251 (MQ=255)                                     aGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGATCTACGCTCCTCTGTTTGTTACGCAATAAAAACTGTTGTTGTCGACTa                                                                                                                                                                                <  2:40291/251‑1 (MQ=255)                                          gACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCa                                                                                                                                                                            <  2:76541/250‑1 (MQ=255)                                                     tATCATGCCCAAACAAGGCGGGTTAACGGCGGCAAAGGTCCTTCCGGTGGCGTCCGACCGGTATACCTAGGGGCAACCGCATACTGCGCTGGGTTATGGCTCGACCCAGAGATTTCTGCGGAAGCGGGCTTGTAATGGGTTCAGTGATAACAGATACCTGGAAATATAGGGGCAAATCCACAGGTGAACTCCGATCCTCTGTTTGTTACGCAATAACAACTGGTGTTTTCGGCTATGTACAACAGTTTccc                                                                                                                                                                >  2:268517/1‑251 (MQ=255)                                                           gCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCctcctc                                                                                                                                                           <  1:163859/250‑1 (MQ=255)                                                                           ttAACGGCAGCAAAGAACCTTGCCGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTACTCTGTTTGTTACGCAATAAAAACTGGTGTTTACGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAAt                                                                                                                                          >  2:56943/1‑251 (MQ=255)                                                                               cGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTccc                                                                                                                                                                >  1:81389/1‑225 (MQ=255)                                                                               cGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTccc                                                                                                                                                                <  2:81389/225‑1 (MQ=255)                                                                                      aaaGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGcc                                                                                                                               >  2:3228/1‑251 (MQ=255)                                                                                            ccTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACaa                                                                                                                         <  1:177999/251‑1 (MQ=255)                                                                                                gCAGAGGCGTTCGTGCATTATAACGAATGGCATCCGCATAGGGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCGACGCTGTGTGTTAAGCAATAAAAACTGGGGTTGTCGACTATGCAAAACAGTCTCCCTCCTCACGCATTCTGGATAAATTTGCTGGGGCCCggcacaaat                                                                                                                     >  2:177999/1‑243 (MQ=17)                                                                                                   gagGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTa                                                                                                                   >  1:349789/1‑250 (MQ=255)                                                                                                        gTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATa                                                                                                             <  2:349789/251‑1 (MQ=255)                                                                                                             aGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGa                                                                                                        <  2:213140/251‑1 (MQ=255)                                                                                                               cATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAAc                                                                                                      >  1:82649/1‑251 (MQ=255)                                                                                                                                             cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTgaaagaa                                                                         >  1:185805/1‑250 (MQ=255)                                                                                                                                                                gggAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACACAATAAAAACTGGTTTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGTTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAACCTGAAGTCGGCAAAACATc                                                     <  2:171457/251‑1 (MQ=255)                                                                                                                                                                              cAGCGGGCTTGTAATGGGTTAAATGATAACAGATGTCTGGTAATATAGCAGAAAATCCACAGGTGAACTACTCTCCTCTGTTTGATACGCACTAAACACCGGTGGTTTCGACTATGCAAAACAGTTACCCTCCTCACGCAATCAGGAAAAATTTGGTGAGGACCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTTTGGCACTGAAAGAAACTGAAGTCGCCAAAAAATCACTCAACACATAcc                                       >  1:228943/1‑251 (MQ=255)                                                                                                                                                                                          aaTGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGtcatc                            <  1:353761/250‑1 (MQ=255)                                                                                                                                                                                                               tGTCTCGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCg      >  1:40563/1‑251 (MQ=255)                                                                                                                                                                                                                    ggAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGa  <  1:3228/250‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      GTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACAGGTGAACTACGCTCCTCTGTTTGTTACGCAATAAAAACTGGTGTTTTCGACTATGCAAAACAGTTTCCCTCCTCACGCAATCTGGAAAAATTTGGTGAGGCCCGACAAGATTTAACCATAAAAGAACTGGCTGAAAAATTTCTGGCACTGAAAGAAACTGAAGTCGCCAAAACATCACTCAACACATACCGTGCCGTCATCAAAAATATCCTGAGCATAATCGGTGA  >  NC_000913/1651317‑1651777

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