Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,633,913	IS5 (–) +4 bp	coding (292‑295/591 nt)	pinQ →	putative recombinase PinQ

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1633913 =	0 (0.000)	17 (0.890)	14/276	0.2	100%	coding (292/591 nt)	pinQ	putative recombinase PinQ
	?	NC_000913	= 2290113	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region
*	?	NC_000913	= 1633916	NA (NA)	20 (1.030)	16/278	0.1	NA	coding (295/591 nt)	pinQ	putative recombinase PinQ
	?	NC_000913	2101749 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region

CGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGAT  >  NC_000913/2289891‑2290325                                                                                                                                                                                                                               |                                                                                                                                                                                                                     cGGCCATTACCAGCCTTGGGGTAAAACGGCTCGCTGAATTACACCATGTTTTCCCATGGCAGGATCCGCTCCATGCGGGACAAGCAGATCTCTTTTCTGGTTTGACGGCGCGTCCTGCTGACTTCACTGTCGCCGGAGGTAAGTTGCTGACTCATGATGAACCCTGGTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAg                                                                                                                                                                                                          >  1:221764/1‑235 (MQ=11)         aCCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGAAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGCTGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATa                                                                                                                                                                                                                   >  1:77972/1‑218 (MQ=14)                                        tCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAAc                                                                                                                                                                     <  1:280561/233‑1 (MQ=255)                                                                                                        gACGGCGCTTACTGCTGAATTCACTGTCTGCGGAGGTAAGTTGATGACTCATGCTGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTTCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGAAAAAGAGAAACAGAACCAGATTGATGCATTg                                                                                                                                                 >  1:480546/1‑189 (MQ=21)                                                                                                            gcgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGNACCCTGGTCNNTGNCTCCANANNNCNNNCANGATCTCATATCAGGGACTTNTTNNCACCTTNCATAACGCTGTAGCCACCAGAACAGa                                                                                                                                                                                             <  2:116885/141‑1 (MQ=17)                                                                                                            gcgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGa                                                                                                                                                                                             >  1:116885/1‑141 (MQ=255)                                                                                                                                                 tgatgaCTCATGATGAACCCTGTTCTATGGCCCCAGATGACAAACATGATCTCATATCAAGGACTTGTTCGCACATTACATagt                                                                                                                                                                                                                 >  1:445128/1‑82 (MQ=255)                                                                                                                                                   atgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATCTCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGtt                                                    >  1:196931/1‑239 (MQ=255)                                                                                                                                                         cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTa                                                             <  1:234543/224‑1 (MQ=255)                                                                                                                                                             atgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTGCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTGCAGCCTATGAAATGAATTGCTGTTAAAAGCATTGTGTACAGAAAATTCCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTGTAAAGAAATGCACAAGTATtgt                              <  1:369510/251‑1 (MQ=39)                                                                                                                                                                                tCCAGNTGACANANNTGATCTCATATCAGGGACTNGTTCGCACCTTCCATNACGCNGTAGcnacna                                                                                                                                                                                                    >  2:191481/1‑66 (MQ=17)                                                                                                                                                                                tCCAGATGACAAACATGATCTCATGTCAGTGACGTGTTCTCACCTTCCATAACGCTTTAGCCACCAGAACAGAGATTGCGGAACGACAAAGAGAAACAGAACCAGATTGAt                                                                                                                                                       <  1:191481/111‑1 (MQ=37)                                                                                                                                                                                 ccAGATAACAAACATGCTCTCATATCAGCGACTTGTTCGCACCTTCCATAACGCTGTAGCAACAAGAACAGATATTGCGGAACGACAAAGAGAACCAGAACCAGATTGAt                                                                                                                                                       >  1:424092/1‑110 (MQ=255)                                                                                                                                                                                         cAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTCGCCACCAGAACCGATATTGCGGAACGACAAAGAGAAACAGCACCCGATTGATGCATTGACCTTTCATCCTATg                                                                                                                                  >  1:376974/1‑123 (MQ=37)                                                                                                                                                                                         cAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTtgat  >  1:165857/1‑251 (MQ=255)                                                                                                                                                                                         cAAACATGATCTCATATCAGGGACTTGTNCGCACCTTCCNTAANNCTGTNGNCNCCAGAACAGATATTGCGGAACGACAAAGANAAACAGAACCAGATTGATGCATTNNTCNNNCTNNCTATg                                                                                                                                  <  2:376974/123‑1 (MQ=37)                                                                                                                                                                                                              gACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGCAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTCATTGCTGTTa                                                                                                                  >  1:191761/1‑118 (MQ=255)                                                                                                                                                                                                                   gTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAg                <  1:290250/211‑1 (MQ=255)                                                                                                                                                                                                                      cGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACACAGAGAAACAGAACCAGATTCATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATtgtg                             >  1:236321/1‑195 (MQ=255)                                                                                                                                                                                                                      cGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAgttgtt           <  1:141261/213‑1 (MQ=255)                                                                                                                                                                                                                        cACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACa                                                                                                   >  1:401998/1‑123 (MQ=255)                                                                                                                                                                                                                          ccTTCCATAACGCTGTAGCAACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTCAAAGCATTGGGTACAGAAAATACCCATa                                                                                      >  1:407393/1‑134 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                     CGGCCATTACCAGCCTTGGGGTAAAACGGCTCGATGACTTCCACCATGTTTTGCCATGGCAGAATCTGCTCCATGCGGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCATAACGCTGTAGCCACCAGAACAGATATTGCGGAACGACAAAGAGAAACAGAACCAGATTGATGCATTGAGCTTTCATCCTATGAAATTAATTGCTGTTAAAAGCATTGGGTACAGAAAATACCCATAGCTCCATACCCGGAGTCAGTTTTTAAAAACTGTTTAAAGAAATGCACAAGTATTGTGATTGATTTTTTAGTTGTTTTTCTTGAT  >  NC_000913/2289891‑2290325

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.