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

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1979486 =	0 (0.000)	8 (0.450)	7/264	NT	100%	intergenic (‑271/‑267)	insA‑5/uspC	IS1 protein InsA/universal stress protein C
	?	NC_000913	2101749 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region
*	?	NC_000913	= 1979489	0 (0.000)	26 (1.450)	21/264	NT	100%	intergenic (‑274/‑264)	insA‑5/uspC	IS1 protein InsA/universal stress protein C
	?	NC_000913	= 2290113	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region

GCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTT‑ACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCC‑AGTTTTTAGTGAGATCTCTCCCAC  >  NC_000913/2101587‑2101913                                                                                                                                                                    |                                                                                                                                                                      gCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTT‑ACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCCTGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTTAAACAACCTAtttctaccc                                                       >  2:314035/1‑269 (MQ=255)                        tCACAGTACATAAAGTAACCTTGATCGAAGCCATTT‑ACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCCCCCGTAAATGCCTTGAAt                                                                       >  2:333590/1‑236 (MQ=255)                            aGTACATAAAGTAACCTTGATCGAAGCCATTT‑ACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                                                   <  1:333590/236‑1 (MQ=255)                                        caCCATGGTCAAAGCCCTTTAACCCGAACAAAATTTTAAAAACGTACCAGCTTAAATGATCCGGCGAACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCAGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCC‑AGTTTTTAGt                <  2:316847/274‑1 (MQ=18)                                        aaCCTTGATCGAAGCCATTT‑ACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCAGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCC‑AGtttt                    >  1:316847/1‑270 (MQ=255)                                               aTCGAAGCCATTT‑ACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCCCCCGTAAATGCCTTTAATCAGCCTATTTAGACCGTTTCTTCCCCATTTAAGGCGTTATCCCCAAGttttttttaaaaatt         >  2:355543/1‑265 (MQ=255)                                                     gCCATTT‑ACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAAAAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGGAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGACGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTTGCCATTTAAGGCGTTTTTCCC‑AGTTTTTAATTAAATCTCTCCCAc  >  1:179962/1‑275 (MQ=255)                                                      ccATTT‑ACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCC‑AGTTTTTAGTGAGATCTCTCCCAc  <  1:355543/274‑1 (MQ=255)                                                                  cacaAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCg                                                                                     >  1:315783/1‑181 (MQ=255)                                                                      aaaTCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaaa                                                                                 <  2:315783/181‑1 (MQ=255)                                                                                            tAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaaa                                                                                 >  1:351093/1‑159 (MQ=255)                                                                                                tGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGcc                                                                             <  2:351093/159‑1 (MQ=255)                                                                                                    ccTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtca                                                                                                                          >  2:286802/1‑110 (MQ=255)                                                                                                        cgcACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCggg                                                                                                                                                    >  2:173117/1‑80 (MQ=255)                                                                                                        cgcACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttcatt                                                                                                                   >  2:99431/1‑113 (MQ=255)                                                                                                        cgcACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttc                                                                                                                      <  1:286802/110‑1 (MQ=255)                                                                                                            ccAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGaaa                                                                                                                                                <  1:173117/80‑1 (MQ=255)                                                                                                            ccAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttctt                                                                                                               <  1:99431/113‑1 (MQ=255)                                                                                                               aTCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGtt                                                                                                         >  2:246568/1‑116 (MQ=255)                                                                                                                   acaacCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAg                                                                                                     <  1:246568/116‑1 (MQ=255)                                                                                                                          tATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaaa                                                                                 >  1:151369/1‑129 (MQ=255)                                                                                                                          tATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTg                                                                          >  1:267445/1‑136 (MQ=255)                                                                                                                              aGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGcc                                                                             <  2:151369/129‑1 (MQ=255)                                                                                                                              aGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATc                                                                      <  2:267445/136‑1 (MQ=255)                                                                                                                                     aTACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCggg                                                                                                                                                    >  2:5011/1‑51 (MQ=255)                                                                                                                                         tttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGaaa                                                                                                                                                <  1:5011/51‑1 (MQ=255)                                                                                                                                                                    |                                                                                                                                                                      GCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTT‑ACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCC‑AGTTTTTAGTGAGATCTCTCCCAC  >  NC_000913/2101587‑2101913

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