Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,293,268	IS5 (–) +4 bp	intergenic (‑346/‑256)	hns ← / → tdk	DNA‑binding transcriptional dual regulator H‑NS/thymidine/deoxyuridine kinase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1293268 =	0 (0.000)	20 (0.660)	19/488	0.5	100%	intergenic (‑346/‑259)	hns/tdk	DNA‑binding transcriptional dual regulator H‑NS/thymidine/deoxyuridine kinase
	?	NC_000913	= 2290113	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region
*	?	NC_000913	= 1293271	0 (0.000)	14 (0.460)	13/490	0.9	100%	intergenic (‑349/‑256)	hns/tdk	DNA‑binding transcriptional dual regulator H‑NS/thymidine/deoxyuridine kinase
	?	NC_000913	2101749 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region

GGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCC  >  NC_000913/2289908‑2290346                                                                                                                                                                                                              |                                                                                                                                                                                                                                          ggggTAAAACGGCTCGATGACTTCCACCATGTTTGGCCATGGCAGAATCTGCTCCATGCTGGACAATAAAATCTCTTGTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACAGGATCTCAGATCAGGGACTTGTTCGCACCTTTAATAACGCTGTAGCCACCAGTTCAGATATTGCTGAACGTCAAaga                                                                                                                                                                                               <  2:253044/250‑1 (MQ=14)                                      cATGGCAGAATCTGCTCCATGCGTGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGATCTTTCATCCTa                                                                                                                                                         <  1:222001/251‑1 (MQ=255)                                                      ccATGCGGGCCAAGCAAATCTCTTTTCTGGTCTGACGGCGCTTAATGCTGAATTCACTGTCGGCGAAGGAAAGTTGATGACTCATGATGAACCCTGTTCTATGGATCCAGATGACAAACATGATCTAATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCAACAGa                                                                                                                                                                                                                       >  2:290089/1‑173 (MQ=14)                                                      ccATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGTCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCCCCTTCCATAACGCTGTAGCCACCAGa                                                                                                                                                                                                                       <  1:290089/173‑1 (MQ=255)                                                      ccATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCt                                                                                                                                           <  2:72141/249‑1 (MQ=255)                                                      ccATGCGGGACAAGAAAATCTCTTTGCTGGTCTGACTGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATTAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCCCCTTCCATAACTCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTg                                                                                                                                          <  2:255347/250‑1 (MQ=255)                                                       cATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGtt                                                                                                                                        <  2:46283/251‑1 (MQ=255)                                                             ggACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGc                                                                                                                                  <  1:113689/251‑1 (MQ=255)                                                               aCAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCAt                                                                                                                                <  1:5029/251‑1 (MQ=255)                                                                                  ggTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCAGAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCAGCCTATGAAAGTAATTGCTTTTAAAATCATTGGGTACAGAAAATAccc                                                                                                              <  2:306117/250‑1 (MQ=255)                                                                                                    gctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGt                                                                                                                                                                                                                                 >  1:251053/1‑117 (MQ=40)                                                                                                    gctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGt                                                                                                                                                                                                                                 <  2:251053/117‑1 (MQ=40)                                                                                                                gTCGGCGAAGGTAATGTGATGACACATGATGAACCCTGTTCTATTGCTCCAGATGACAAACATGATCTCATATCAGGGACGTGTGCGCACCTTCCATAACGCTGGAGCCACCAGAACAGATATTTCGGAGCGACAAAGAGAAGCAGAACCAGATTGATGCATGGAGCTTCTATCCTATTAAATTAATTTCTGTTAAAAGCTTTTGATAGAGAAGTTACCCATAGCTCCACACCCGGGGTCatctgctaaaa                                                                               <  2:137488/251‑11 (MQ=16)                                                                                                                            aaGTTGATAACTAATGATGTACACTGTTCTATGGCTCCAGATGCAAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTCTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTAATCCTATGAAATTAATTGCTGTAAAAAGCATTGGGTACAGAAAATACACATATCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTaa                                                                      >  2:270828/1‑248 (MQ=38)                                                                                                                                                                                  ctcATATCAGGGCCTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTg              >  1:251652/1‑250 (MQ=255)                                                                                                                                                                                         cAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACCAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCaaa       >  1:50079/1‑250 (MQ=255)                                                                                                                                                                                             gACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTcc  <  1:219249/251‑1 (MQ=255)                                                                                                                                                                                                              |                                                                                                                                                                                                                                          GGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCC  >  NC_000913/2289908‑2290346

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