Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913	682,848	22 bp→TTTTTTTTTTAATGATACGGCGA	12.4%	coding (525‑546/1671 nt)	hscC ←	chaperone protein HscC

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 682847	78 (0.830)	8 (0.120) +TTTTTTTTTTAATGATACGGCGA	7/130	NT	12.4%	coding (547/1671 nt)	hscC	chaperone protein HscC
	?	NC_000913	682870 =	75 (0.790)					coding (524/1671 nt)	hscC	chaperone protein HscC

AAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGTACTCAAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/682795‑682847 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CCACCGAGATCAAACACCAGCGAACGGGTATTTTGTTGGGTATGCAGGCCATACGCCATCGCAGCTGCTGTGGGTTCATTAATTAAGCGTAC  >  NC_000913/682870‑682961                                                      |||||||||||||||||||||||                                                                                             AAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGTACTCAAGTTTTTTTTTTAATGATACGGCGCCCACCGAGATC‑TACACCAGCGAA                                                                      <  2:966212/99‑1            AGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGTACTCAAGTTTTTTTTTTAATGATACGGCGACCACCGAGATC‑TACACCAGCGAACGGGTATTTTG                                                           <  2:1895553/99‑1                   GCGTGAACTTCAATCACCGGCGTGGCGTACTCAAGTTTTTTTTTTAATGATACGGCGACCACCGAGATC‑TACACCAGCGAACGGGTATTTTGTTGGGTA                                                    <  2:1541646/99‑1                      TGAACTTCAATCACCGGCGTGGCGTACTCAAGTTTTTTTTTTAATGATACGGCGACCACCGAGATC‑TACACCAGCGAACGGGTATTTTGTTGGGTATGC                                                 <  2:2371660/99‑1                            TCAATCACCGGCGTGGCGTACTCAAGTTTTTTTTTTAATGATACGGCGCCACCCGAGTTC‑TACACCAGCGAACGGGTATTTTGTTGGGTATGCAGGCCA                                           <  2:599893/99‑1                            TCAATCACCGGCGTGGCGTACTCAAGTTTTTTTTTTAATGATACGGCGACCACCGAGATC‑TACACCAGCGAACGGGTATTTTGTTGGGTATGCAGGCCA                                           <  2:2074880/99‑1                             CAATCACCGGCGTGGCGTACTCAAGTTTTTTTTTTAATGATACGGCGACCACCGAGATC‑TACACCAGCGAACGGGTATTTTGTTGGGTATGCAGGCCAT                                          <  2:2206922/99‑1                                    CGGCGTGGCGTACTCAAGTTCTTTTATTAATAATACGGCGACCACCGAGATC‑TACACCAGCGAACGGGTATTTTGTTGGGTATGCAGGCCATACGCCAT                                   <  2:787107/99‑1                                                        TTTTTTTTAATGATACGGCGACCACCGAGATC‑TACACCAGCGAACGGGTATTTTGTTGGGTATGCAGGCCATACGCCATCGCAGCTGCTGTGGGTTCAT               <  2:2534880/99‑1                                                                     TACGGCGACCACCGAGATC‑TACACCAGCGAACGGGTATTTTGTTGGGTATGCAGGCCATACGCCATCGCAGCTGCTGTGGGTTCATTAATTAAGCGTAC  <  2:475387/99‑1                                                      |||||||||||||||||||||||                                                                                             AAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGTACTCAAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/682795‑682847 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CCACCGAGATCAAACACCAGCGAACGGGTATTTTGTTGGGTATGCAGGCCATACGCCATCGCAGCTGCTGTGGGTTCATTAATTAAGCGTAC  >  NC_000913/682870‑682961

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