Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC JC	NC_000913	4,542,428	IS1 (+) +9 bp	11.2%	coding (392‑400/597 nt)	fimE →	regulator for fimA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3583428 =	NA (NA)	13 (0.100)	12/176	NT	9.8%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	= 4542436	107 (0.910)					coding (400/597 nt)	fimE	regulator for fimA
*	?	NC_000913	= 3584195	NA (NA)	14 (0.110)	12/176	NT	10.5%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	4542428 =	107 (0.910)					coding (392/597 nt)	fimE	regulator for fimA

ATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3584120‑3584195 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGA  >  NC_000913/4542428‑4542525                                                                                                                                                                                ATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCGCAGACTCATCCTCATATGTTA                                                                             >  1:2719359/1‑99    GGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCGCAGACTCATCCTCATATGTTAAGG                                                                          <  2:2689738/99‑1         AGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCGCAGACTCATCCTCATATGTTAAGGCATGC                                                                     <  2:1933787/99‑1              CATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCG                                                                >  2:1383159/1‑99                      AGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAA                                                        >  1:3433476/1‑99                           TCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGC                                                   >  1:188124/1‑99                                  TTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGT                                            >  1:940169/1‑99                                     TCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGT                                         >  2:562129/1‑99                                         AACATAAAACACTATCAATAAGTTGGAGTCATTACCCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAG                                     >  2:638434/1‑99                                          ACATAAAACACTATCAATAAGTTGGAGTCATTACCCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGA                                    <  2:70560/99‑1                                           CATAAAACACTATCAATAAGTTGGAGTCATTACCCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGAT                                   >  1:1521860/1‑99                                           CATAAAACACTATCAATAAGTTGGAGTCATTACCCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGAT                                   >  2:858368/1‑99                                           CATAAAACACTATCAATAAGTTGGAGTCATTACCCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGAT                                   >  2:865869/1‑99                                                                       ATTACCCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGC       >  1:1428761/1‑99                                                                            CCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGA  <  1:249669/99‑1                                                                                                                                                                                ATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3584120‑3584195 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGA  >  NC_000913/4542428‑4542525

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