Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,319,850	IS2 (+) +5 bp	coding (23‑27/2850 nt)	rcsC ←	histidine kinase RcsC

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2068941 =	NA (NA)	9 (0.940)	8/314	0.1	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	2319850 =	0 (0.000)					coding (27/2850 nt)	rcsC	histidine kinase RcsC
*	?	NC_000913	= 2319854	0 (0.000)	12 (1.250)	11/314	0.0	100%	coding (23/2850 nt)	rcsC	histidine kinase RcsC
	?	NC_000913	= 2997691	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

GCTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGACAA  >  NC_000913/2997447‑2997933                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   gCTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTGTCCTGTGTGGTACGCCGTCTGCGGTGCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTAGTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTg                                                                                                                                                                                                                                              <  1:120826/251‑1 (MQ=9)       gCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCttt                                                                                                                                                                                                                                        >  1:36238/1‑251 (MQ=33)           gagaAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTttctt                                                                                                                                                                                                                                    >  1:182721/1‑251 (MQ=35)           gagaAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTttctt                                                                                                                                                                                                                                    >  1:182716/1‑251 (MQ=35)                  aaCTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGgtgt                                                                                                                                                                                                                             >  1:105587/1‑251 (MQ=38)                         tcgggagaaggcacgctgaccGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTGTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTAGTTAAGTGATATTTGTTGCCTTGAGATTCATGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATc                                                                                                                                                                                                                        <  1:207134/230‑1 (MQ=14)                                                   aCGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAAGCGTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGTGCCAGTCTAAACTTGCTCTGTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACAt                                                                                                                                                                                            <  1:107699/251‑1 (MQ=255)                                                              cGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGcn                                                                                                                                                                                 <  1:53645/251‑2 (MQ=255)                                                                            cAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCttt                                                                                                                                                                                                                     <  1:91432/201‑1 (MQ=255)                                                                                    tgctgAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGAGCCAGTCTAAACTTGCTCTTTTCGTCGTGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACAt                                                                                                                                                                                            <  1:100046/218‑1 (MQ=255)                                                                                                                                    tcCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTCCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGt                                                                                                            >  1:81244/1‑250 (MQ=255)                                                                                                                                                                              cTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGAtttn                                                                                                                                                                                                           <  1:130672/113‑2 (MQ=255)                                                                                                                                                                              cTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTg                                                                                                                                                                                                           >  2:130672/1‑113 (MQ=255)                                                                                                                                                                                aaaaaCTAGGATTAAGACTATCACTTNTTTAAGNNATNTTGGTTGTCTGGAGATTCANGGNNNCNNNCNAAACTNGCNNTTTTCTTCTNGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCNGCTTCTCTGNNATGTTTATCTTTTTCCAGATANCCTNNGNNTNCTcncc                                                                                                                          >  2:45481/1‑192 (MQ=38)                                                                                                                                                                                                     cACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGGTTATCTGTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATATCCTGGGGAGTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTGGTTGAAAGAGTGTCATnc                                          <  1:169995/251‑1 (MQ=255)                                                                                                                                                                                                                         gTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATc                                                                                                                   >  1:49527/1‑158 (MQ=255)                                                                                                                                                                                                                         gTTGNNTGGAGATTCNGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTNTAGGATGAANNACNTTAAGNANNNCNNNATNGGCTTCTCTGTAATGTTTANCTNNTTCCAGNTAACCTGGGGATTCTCCCCCTAAATc                                                                                                                   <  2:49527/158‑1 (MQ=255)                                                                                                                                                                                                                                       cAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCgta                 >  1:87575/1‑242 (MQ=255)                                                                                                                                                                                                                                           gggCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGAc    >  1:10728/1‑251 (MQ=255)                                                                                                                                                                                                                                            ggCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTATGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGAc    >  1:147610/1‑250 (MQ=255)                                                                                                                                                                                                                                              ncAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGACaa  >  1:62883/2‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   GCTTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACTTGCTCTTTTCTTCTGGTGTTTATCTTTGAAAGATTTGTAGGATAAATAACATTAAGTAAAGCCGCATCGGCTTCTCTGTAATGTTTATCTTTTTCCAGATAACCTGGGGATTCTCCCCCTAAATCAAGCTGTAACGTTATTCTTTCATAACCTTTATCTATGAGGTGGTTCAACCTCGTTGTTGAAAGAGTTTCATCCAGTTCATTGACTTTAGCTTCGCGTAGTAAATCAGAGACAA  >  NC_000913/2997447‑2997933

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