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)	34 (0.760)	32/144	0.4	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)	45 (1.010)	36/144	0.2	100%	noncoding (1343/1343 nt)	IS186	repeat region
	?	REL606	4524522 =	0 (0.000)					coding (494/549 nt)	fimA	major type 1 subunit fimbrin (pilin)

GAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAT  >  REL606/2772781‑2772912                                                                          |                                                           gAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCGGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAcccc                                                          <  1:793483/76‑1 (MQ=255)   aaTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCAt                                                        <  1:634595/76‑1 (MQ=255)     ttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGGACACCCCGAACGGGGGGGGAAAAACCCATaa                                                      >  1:21838/1‑76 (MQ=25)           ttCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGGACACCCTGAACGGGGGGGGAAAACCCCATAAGCGCTa                                                >  1:1444404/1‑76 (MQ=38)           ttCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                                                <  2:669008/76‑1 (MQ=255)                 gTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTaa                                          <  2:1448578/76‑1 (MQ=255)                  tgatgaCCGGTTTACCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAg                                         <  1:431953/76‑1 (MQ=255)                       aCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtg                                    <  2:590474/76‑1 (MQ=255)                        ccGGTTTATCCCCCCCTGGGCGGGGGACACCCCGAAAGGGGGGGGAAAAACCCCCAAACGGCAACTTTAGGGGtgg                                   >  2:901231/1‑75 (MQ=255)                         cGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtgtg                                  <  1:60785/76‑1 (MQ=255)                          ggTCTATCCCCCCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTgg                                 <  2:1006838/76‑1 (MQ=255)                             ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                              <  2:1235404/76‑1 (MQ=255)                               ccccccGCGGGGGGGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                            <  1:658464/74‑1 (MQ=25)                                       tGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                    >  1:713125/1‑76 (MQ=255)                                       tGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                    <  1:267460/76‑1 (MQ=255)                                         gcgcGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                  <  1:292182/76‑1 (MQ=255)                                             ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt              <  1:1083389/76‑1 (MQ=255)                                                gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa           <  1:230278/76‑1 (MQ=255)                                                 aaCACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc          <  1:888126/76‑1 (MQ=255)                                                   cacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt        >  1:398782/1‑76 (MQ=255)                                                    acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg       >  2:234171/1‑76 (MQ=255)                                                    acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg       >  2:825619/1‑76 (MQ=255)                                                         gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt  >  1:1053108/1‑76 (MQ=255)                                                         gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt  <  1:564248/76‑1 (MQ=255)                                                                          |                                                           GAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAT  >  REL606/2772781‑2772912

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