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)

AATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAC  >  REL606/2772783‑2772919                                                                        |                                                                  aaTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCAt                                                               <  1:421015/76‑1 (MQ=255) aaTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCAt                                                               <  1:857187/76‑1 (MQ=255)     ttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAgc                                                           >  2:909234/1‑76 (MQ=255)        aTTGATCGTGATGNCCGGTTNATCCCGGCTGGCGCGNGGAACACCCTGAACNTTGGGGAAAANCCCCATAAGCGCt                                                        <  2:827452/76‑1 (MQ=255)         ttCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                                                       >  2:89641/1‑76 (MQ=255)              cGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                  <  2:212661/76‑1 (MQ=255)              cGTGATGACCGGTTTATCCCCCCCGGGGCGGGGGACACCCCGAACGGGGGGGGAAAAACCCCCAAACGCTAACTTa                                                  >  2:632058/1‑76 (MQ=255)                tgatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAg                                                <  2:262705/76‑1 (MQ=255)                  atgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAggg                                              >  2:326029/1‑76 (MQ=255)                    gaCCGGTTTATCCCCGCTGGCGCGGGGATCACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                            <  2:440444/76‑1 (MQ=255)                      ccGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAATAACCCCATAAGCGCTAACTTAAGGGTtgt                                          <  1:949755/76‑1 (MQ=255)                        ggTTTATCCCCGCTGGCGCGGGGGACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGGTGTgg                                        >  2:47776/1‑76 (MQ=255)                         gTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt                                       <  2:564870/76‑1 (MQ=255)                         gTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt                                       <  2:407184/76‑1 (MQ=255)                             aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                   <  2:117577/76‑1 (MQ=255)                              tCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                  <  1:616916/76‑1 (MQ=255)                               ccccGCTGGCGCGGGGAAAACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                 >  2:874812/1‑76 (MQ=255)                                     tagcgcGGGGAACACCCTGAACAGTGGGGTAAAACCCAATAAGCGCTAACTTAAGGGTTGTGGTATTACACCTGat                           >  1:830299/3‑76 (MQ=21)                                          cgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                      >  2:604299/1‑76 (MQ=255)                                            gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                    <  1:319053/76‑1 (MQ=255)                                               aaCACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAGCTTAAGGGTTGTGGCATTACGCCTGATATGATTTAAc                 >  1:664379/1‑76 (MQ=255)                                                  acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg              <  1:47776/76‑1 (MQ=255)                                                     cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTNATATGATTTAACGTGCCg           >  1:756512/1‑76 (MQ=255)                                                              ggggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc  >  1:499770/1‑76 (MQ=39)                                                                        |                                                                  AATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAC  >  REL606/2772783‑2772919

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