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

ATTTAACACAGGTGGTCGACCAAAACGTTTCCCTGCGCCGCGGGCTCTTACTATCCCGGAATGAGTGCGTTCAAGTAAAAGGTCTCGTTCAAATTCAGCGACTGCTGAAATTACTTGCATCATCATTTTTCCTGTTGGACTGGTCAGGTCAATGCCCCCCAATGCTAAg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1634080‑1633913 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGC  <  NC_000913/2290113‑2289893                                                                                                                                                                                                                                                                                                                                                                                                       ATTTAACACAGGTGGTCGACCAAAACGTTTCCCTGCGCCGCGGGCTCTTACTATCCCGGAATGAGTGCGTTCAAGTAAAAGGTCTCGTTCAAATTCAGCGACTGCTGAAATTACTTGCATCATCATTTTTCCTGTTGGACTGGTCAGGTCAATGCCCCCCAATGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTT                                                                                                                                            <  1:172142/251‑1      nCACAGGTGGTCGACCAAAACGTTTCCCTGCGCCGCGGGCTCTTACTATCCCGGAATGAGTGCGTTCAAGTAAAAGGTCTCGTTCAAATTCAGCGACTGCTGAAATTACTTGCATCATCATTTTTCCTGTTGGACTGGTCAGGTCAATGCCCCCCAATGCTAAGGCAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTT                                                                                                                                       >  1:58986/2‑251                 CGACCAAAACGTTTCCCTGCGCCGCGGGCTCTTACTATACCGGAATGAGTGCGTTCAAGTAAAAGGTCTCGTTCAAATTCAGCGACTGCTGAAATTACTTGCATCATCATTTTTCCTGTTGGACTGGTCAGGTCAATGCCCCCCAATGCTAAGGAAGGTGCGAACAAGTCCCTGATA                                                                                                                                                                                                      >  1:225628/1‑177                                                       CCCGGAATGAGTGCGTTCAAGTAAAAGGTCTCGTTCAAATTCAGCGACTGCTGACATTACTTGCATCATCATTTTTCCTGTTGGACTGGTCAGGTCAATGCCCCCCACTGCTAAGGAAGGTGCG                                                                                                                                                                                                                     >  1:88628/1‑124                                                         CGGAATGAGTGCGTTCAAGTAAAAGGTCTCGTTCAAATTCAGCGACTGCTGAAATTACTTGCATCATCATTTTTCCTGTTGGACTGGTCAGGTCAATGCCCCCCAATGCTAAGGAAGGTGCGAACAAGTCCCTGAT                                                                                                                                                                                                       >  1:201275/1‑136                                                               GAGTGCGTTCAAGTAAAAGGTCTCGTTCAAATTCAGCGACTGCTGAAATTAATTGAATCATCATTTTTCCTGTTGGCCTGGTCAGGTCAATGCACCCCAATGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATG                                                                                                                                                        >  1:323038/1‑177                                                                                    CTCGTTCAAAGTCAGCGACTGCTGAAATTACTTGCATCATCATTTTTCCGGTTGGACTGGGCAGGTCAATGCCCCCCAATGCTAAGGAAGTTGCGAACAATTCCCTGATATGAGATCATGTTTGGCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTG                                                                                                                            <  1:434994/184‑1                                                                                      CGTTCAAATTCAGCGACTGCTGAAATTACTTGCATCATCATTTTTCCTGTTGGACTGGTCAGGTCAATGCCCCCCAATGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGATCAGATTCTGn                                                          <  1:146832/248‑2                                                                                        TTCAAATTCAGCGACTGCTGAAATTACTTGCATCATCATTTTTCCTGTTGGACTGGTCAGGTCAATGCCCCCCAATGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAG                                                                                                                  <  1:410196/190‑1                                                                                                          TGAAATTACTTGCATCATCATTTTTCCTGTTGGACTGGACAGGTCAATGACCCCCAATGCTAAGGAAGGCGCGAAAAAGTCCa                                                                                                                                                                                                           >  1:205397/1‑82                                                                                                                                   CCTGTTGGACTGGTCAGGTCAATGCCCCCCAATGCTACGGAAGGTGCGAACAAGTCCATGATCTGCGATCATGTTTGTCATCTGGAGCCATCGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGACTTCAGCAGTACGCGCCGTCAGACCAGAACAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCACGGCTG          >  1:219373/1‑251                                                                                                                                           ACTGGTCAGGTCAATGCCCCCCAATGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATGAACTTACCTTCGCCGTCAGTGAATTCAGCATTAAGCGCCGCCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATTTCAAAACATGGTGGAAGTCACCGAGCCGTTTTACTCCAAGGCTGGTAATGGC  <  1:151605/251‑1                                                                                                                                           ACTGGTCAGGTCAATGCCCCCCAATGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGACATCGAGCCGTTTTACCCCAAGGCTGGTAAT     >  1:288211/1‑248                                                                                                                                            CTGGTCAGGTCAATGCCCCCCAATGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGG                                        <  1:324245/212‑1                                                                                                                                                    GTCAATGCCCCCCAATGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAANCGCNNNNNGAnnnnnnnnnn                                                                                          >  2:300401/1‑144                                                                                                                                                    GTCAATGACCCCCACTGCGAAGGAAGGTGCGACCAAGTCCCTGATATGGGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCAGCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCATCAGCCCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAAn     <  1:300401/239‑2                                                                                                                                                           CCCCCCAATGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACC                                                                                                                                         >  1:442481/1‑100                                                                                                                                                                                                                                                                                                                                                                                                       ATTTAACACAGGTGGTCGACCAAAACGTTTCCCTGCGCCGCGGGCTCTTACTATCCCGGAATGAGTGCGTTCAAGTAAAAGGTCTCGTTCAAATTCAGCGACTGCTGAAATTACTTGCATCATCATTTTTCCTGTTGGACTGGTCAGGTCAATGCCCCCCAATGCTAAg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1634080‑1633913 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTTGTCCCGCATGGAGCAGATTCTGCCATGGCAAAACATGGTGGAAGTCATCGAGCCGTTTTACCCCAAGGCTGGTAATGGC  <  NC_000913/2290113‑2289893

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