Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,979,486	IS5 (+) +4 bp	intergenic (‑271/‑264)	insA ← / → uspC	IS1 repressor TnpA/universal stress protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1979486 =	0 (0.000)	40 (1.180)	25/260	0.5	100%	intergenic (‑271/‑267)	insA/uspC	IS1 repressor TnpA/universal stress protein
	?	NC_000913	2101749 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region
*	?	NC_000913	= 1979489	0 (0.000)	60 (1.770)	29/260	0.3	100%	intergenic (‑274/‑264)	insA/uspC	IS1 repressor TnpA/universal stress protein
	?	NC_000913	= 2290113	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region

ATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTA  >  NC_000913/2101619‑2101884                                                                                                                                   |                                                                                                                                        aTAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCgg                                                                                                                       >  3:94406/1‑149 (MQ=255) aTAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCgg                                                                                                                       <  4:94406/149‑1 (MQ=255)   aaaGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTaaggaa                                                                                                                                      <  1:23256/132‑1 (MQ=255)   aaaGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTaaggaa                                                                                                                                      >  2:23256/1‑132 (MQ=255)   aaaGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCggg                                                                                                                      >  8:135106/1‑148 (MQ=255)         aCCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAAttc                                                                                                               <  7:135106/149‑1 (MQ=255)                    aGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGa                                                                                                    >  1:5161/1‑149 (MQ=255)                                  acaAAATCTGAAAAGCGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAg                                                                                               >  3:100166/1‑140 (MQ=255)                                  acaAAATCTGAAAAGCGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAg                                                                                               <  4:100166/140‑1 (MQ=255)                                                aCGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGttt                                                                          >  1:183249/1‑147 (MQ=255)                                                aCGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGa                                                                        <  2:5161/149‑1 (MQ=255)                                                                            cAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCtt                                             >  1:12767/1‑148 (MQ=255)                                                                                  acaacCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAg                                      >  5:71144/1‑149 (MQ=255)                                                                                    aacCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGcc                                    >  2:23545/1‑149 (MQ=255)                                                                                                 tAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGtt                       >  8:2690/1‑149 (MQ=255)                                                                                                      aCTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc                  >  1:138937/1‑149 (MQ=255)                                                                                                                 gAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAgg       >  4:99660/1‑149 (MQ=255)                                                                                                                       ttGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTa  <  2:183249/148‑1 (MQ=37)                                                                                                                                   |                                                                                                                                        ATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTA  >  NC_000913/2101619‑2101884

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