Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_007779	2,865,278	+TTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGC	9.3%	coding (813/876 nt)	rpoS ←	RNA polymerase, sigma S (sigma 38) factor

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_007779	2865225 =	159 (0.840)	9 (0.080)	7/104	NT	9.3%	coding (866/876 nt)	rpoS	RNA polymerase, sigma S (sigma 38) factor
	?	NC_007779	= 2865278	146 (0.770)					coding (813/876 nt)	rpoS	RNA polymerase, sigma S (sigma 38) factor

TCCTGCAAACGCAGGGGCTGAATATCGAAGCGCTGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_007779/2865261‑2865225 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATC  <  NC_007779/2865278‑2865236                                      ||||||||||||||||||||||||||||||||||||||                                            TCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAA                   <  1:2146644/101‑1  CCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAAC                  <  1:9352657/101‑1     GCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCA               <  1:9308917/101‑1       AAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGG             >  1:7411437/1‑101       AAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGG             >  1:9406682/1‑101         ACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGG           <  1:2426324/101‑1         ACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGG           <  1:8396434/101‑1                 GCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATAT   <  1:3957532/101‑1                  CTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATCGAAGCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATC  <  1:8506558/101‑1                                      ||||||||||||||||||||||||||||||||||||||                                            TCCTGCAAACGCAGGGGCTGAATATCGAAGCGCTGTT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_007779/2865261‑2865225 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCGCCGTTTGCGCGAAATCCTGCAAACGCAGGGGCTGAATATC  <  NC_007779/2865278‑2865236

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