Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,979,518	IS5 (+) +4 bp	intergenic (‑303/‑232)	insA5 ← / → uspC	IS1 family protein InsA/universal stress protein C

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1397238	NA (NA)	24 (0.950)	23/490	0.2	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	1979518 =	0 (0.000)					intergenic (‑303/‑235)	insA5/uspC	IS1 family protein InsA/universal stress protein C
*	?	NC_000913	= 1979521	0 (0.000)	13 (0.510)	12/492	0.7	100%	intergenic (‑306/‑232)	insA5/uspC	IS1 family protein InsA/universal stress protein C
	?	NC_000913	= 2290113	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region

GGGTTGCTGAAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGTTTTCATTATTGCGGCGTCGGTTGTGACTCGCGCGACCATCGGC  >  NC_000913/1396994‑1397469                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        gggTTGCTGAAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCctaact                                                                                                                                                                                                                                   <  1:10163/251‑1 (MQ=11)               cGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGa                                                                                                                                                                                                                     <  1:132696/251‑1 (MQ=35)                     ccAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGTACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGTGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTGCCCTAACTAATCAATGCGTTGATTGTaa                                                                                                                                                                                                               <  2:14326/251‑1 (MQ=255)                     ccAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTa                                                                                                                                                                                                                >  1:74392/1‑250 (MQ=39)                                       cACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAgaggtgag                                                                                                                                                                                             <  1:174331/251‑1 (MQ=255)                                              gCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGtttgca                                                                                                                                                                                      <  2:247138/251‑4 (MQ=255)                                                 aaCCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGCGGTTTTCag                                                                                                                                                                                     >  2:65304/1‑249 (MQ=255)                                                   ccTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGc                                                                                                                                                                                  >  2:22486/1‑250 (MQ=255)                                                                                              cACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCGAAATAGGCTGATTCAAGGAATTTACGGGAGAAAAAATCGGCTCAAACATTAAGAAATGAAATGACTGAGTCAGCCGAGAAGAAGTTCCCCGCTTATTCGCACCTTCCCGAACTAATCAATGCGTTGATTGTAAATGCAGCTAAGAGGTGATGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGATCCTTTTAGGTAACaa                                                                                                                                      <  2:249375/251‑1 (MQ=18)                                                                                                            gATAACGCCTTAACTGGCGAAGAAACGGTCTAAATCGGCCGATTCAACGCATTTACGGGCGAAACAATCGGCTCACACATGCAGACATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCACTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGACATGACATGGTATTGTTAGGTTTAGGTAACAAGACATTTGTCTGc                                                                                                                         >  2:133846/1‑250 (MQ=21)                                                                                                                                                                        gaAAAAAGCGGCTCAAACATGAAGAAATGAAATGACGGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCTCCTTCCCGAACTAATCAATGCGTCGATTGTAAATCCAGCTAAGAGGTGAGGTTCTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCtt                                                            <  2:1155/251‑1 (MQ=255)                                                                                                                                                                                                aaatgaaatgACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGTTTTCattat                                    <  1:139722/251‑1 (MQ=255)                                                                                                                                                                                                                     cGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGTTTTCATTATTGCGGCGTCGGTTGTGACTcg               <  1:258391/251‑1 (MQ=255)                                                                                                                                                                                                                                  cccGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGTTTTCATTATTGCGGCGTCGTTTGTGACTCGCGCGACCATcggc  <  2:260251/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        GGGTTGCTGAAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATACGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGTTTTCATTATTGCGGCGTCGGTTGTGACTCGCGCGACCATCGGC  >  NC_000913/1396994‑1397469

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