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)	18 (0.520)	17/486	0.7	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)	17 (0.480)	15/488	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

AAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGG  >  NC_000913/2101504‑2101973                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                 aaaaaaaCTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTaaggaagg                                                                                                                                                                                                                             >  1:238352/1‑251 (MQ=255)            gATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGc                                                                                                                                                                                                                  >  1:110055/1‑251 (MQ=255)                            aTAATGTATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCCATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGTGAAATTCTTCTCg                                                                                                                                                                                                  >  2:14669/1‑251 (MQ=255)                                 gTATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTg                                                                                                                                                                                              <  2:265504/250‑1 (MQ=255)                                      gCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAg                                                                                                                                                                                        <  1:18919/251‑1 (MQ=255)                                       cATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt                                                                                                                                                                                       >  1:235169/1‑251 (MQ=255)                                              aGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGACTTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttc                                                                                                                                                                                 >  1:287849/1‑250 (MQ=255)                                                       aTAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCAGTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAAACAACAACCGTATCAGAATAGAGACTGTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCAt                                                                                                                                                                       <  2:115935/251‑1 (MQ=255)                                                                                       tCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTATCCCGTAAATGCCt                                                                                                                                       >  1:164403/1‑251 (MQ=255)                                                                                         aaGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTg                                                                                                                                     <  1:107223/251‑1 (MQ=255)                                                                                            cacaGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAAt                                                                                                                                  >  1:43257/1‑251 (MQ=255)                                                                                                           tCACAGTACATAACGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAcc                                                                                                                   >  2:122062/1‑251 (MQ=255)                                                                                                                                                         aaTCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGAtct                                                                     >  1:296699/1‑251 (MQ=255)                                                                                                                                                                                    aTCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGATTTATTACGGAAGGCGCGAATAAGCGGGGACATTCTTCTCGGCTGACTCAGACCTTTCATTTCTTCCTGTTTGAGCCGATTTTTTCTCCCGTAACTGCCTTGAATCACCCTATTTAGTCCGTTTCTTCGCCATTTAACGCGTTATCCCCAGTTTTTTGTGAGATCTCTTCCACTGACGTATCATTTGTTCCGc                                          >  2:8928/1‑251 (MQ=255)                                                                                                                                                                                           gcACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAAc                                   <  2:110055/251‑1 (MQ=255)                                                                                                                                                                                                                             tatttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACGCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTGCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTgg  <  2:287849/248‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                 AAAAAAACTTCGATTGTCATGATGAGCATAATGTATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGG  >  NC_000913/2101504‑2101973

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