Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	REL606	4,524,522	IS186 (+) +6 bp	coding (494‑499/549 nt)	fimA →	major type 1 subunit fimbrin (pilin)

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	2772854 =	NA (NA)	16 (0.560)	14/144	1.0	100%	noncoding (1/1343 nt)	IS186	repeat region
	?	REL606	= 4524527	0 (0.000)					coding (499/549 nt)	fimA	major type 1 subunit fimbrin (pilin)
*	?	REL606	= 2774196	NA (NA)	22 (0.770)	20/144	0.5	100%	noncoding (1343/1343 nt)	IS186	repeat region
	?	REL606	4524522 =	0 (0.000)					coding (494/549 nt)	fimA	major type 1 subunit fimbrin (pilin)

GCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2774126‑2774196 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAG  >  REL606/4524522‑4524592                                                                                                                                                GCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCA                                                                    <  1:466784/76‑1  CATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGTCCGCAA                                                                   <  1:146479/76‑1  CATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAA                                                                   >  2:765938/1‑76   ATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAAC                                                                  <  1:369459/76‑1             CGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCT                                                        >  2:446756/1‑76             CGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCT                                                        >  1:771604/1‑76                 GCACCGAGGCGTCGCATTCTTCATATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTA                                                    <  1:374446/76‑1                     CGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTAGGGGCCGCAACCCCGGGTGCTGCTAATGC                                                >  2:440316/1‑76                                    TTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGAGTGCTGCTAATGCGGATGCGACCTTCAA                                 >  2:444174/1‑76                                      CAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGG                               >  1:665189/1‑76                                           GGTTCAACCCTTAAGTGAGCGCTTATGGGTCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAG                          <  1:869936/76‑1                                               CAACCCTGAAGTTAGCGCTTATGGGCCGGAACCCCGGGTGCTGCTAATGCGGATGCGTCCTTCAAGGTTCAGTATC                      <  2:876258/76‑1                                                 ACCCTAAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAA                    <  2:755514/76‑1                                                  CCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAT                   >  1:563070/1‑76                                                      TAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACC               >  2:469010/1‑76                                                       AAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCT              <  1:559437/76‑1                                                       AAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCT              >  1:49387/1‑76                                                          TTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACC           >  2:10881/1‑76                                                          TTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACC           >  1:259524/1‑76                                                            AGCGCTTATGGGCCGCAACCCCGGGTGCTGCTACTGCGGATGCGACCTTCAAGGTTCAGTATCAATAGCCTACCCA         <  1:765938/76‑1                                                                CTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTT     >  2:414060/1‑76                                                                   ATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAG  <  2:49387/76‑1                                                                                                                                                GCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2774126‑2774196 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAG  >  REL606/4524522‑4524592

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