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)	15 (1.150)	12/290	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)	15 (1.150)	11/290	0.1	100%	coding (23/2850 nt)	rcsC	histidine kinase RcsC
	?	NC_000913	= 2997691	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

AAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACC  >  NC_000913/2068734‑2069189                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         aaaaTCAGCAACGAAGAGACGCTCGTCGTGACCACGACTCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTCTCCAGACTATTGAGCAACAAACAATGCAGAACTCCACCAATGCATCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACACCTCCTTGGCCGGCCGTACACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCAAGACATCTGTTATCACTTAAccc                                                                                                                                                                                                               >  1:287875/1‑251 (MQ=255) aaaaTCAGCAACGAAGAGACGCTCGGCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAAccn                                                                                                                                                                                                               <  1:28349/251‑2 (MQ=255)                            gtGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAACCCGCTGCCTCAGATATTc                                                                                                                                                                                    <  1:92772/251‑1 (MQ=255)                                    gAATCGTTCGGCAAGCAACCTTTGGGAAAGCACGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGCCATCTGTTATCACTTAACCCATTACACGCCCGCTGCCGCAGATATTCCCGTGgcg                                                                                                                                                                            >  1:292824/1‑251 (MQ=255)                                        cGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgcgagcn                                                                                                                                                                        <  1:203518/251‑2 (MQ=255)                                                             gAAAGCCCGGCGACTATACAGGTTATTGCCCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                                                                                                                                                   >  1:105849/1‑203 (MQ=255)                                                                                                              ccaccaATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTGGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGTGATAGCCCGGCGCACTATGCGCATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTtgctgc                                                                                                   <  1:310243/250‑1 (MQ=255)                                                                                                                                                     ccGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGAAAAAAAATCCGCATTCGTGGCTGGATTTGCCCCTGTATTTCCAGACATCTGTTATCACTTAACCCATTCCAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAAc                                                                                            >  1:293367/1‑218 (MQ=255)                                                                                                                                                                                      gCCGTCAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTCCAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCCAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTAtt                          >  1:249769/1‑251 (MQ=255)                                                                                                                                                                                                          cttGGCTGGATTTGCCCCTATATTTCCAGACANCTGTTATCACNTAANNCATANCNANCCCGCTGCCGCAGATATTCCCGTGGCGAGNGATAACCCAGCGCACTATGCGGANNCCNNNCGNNATAAt                                                                                                                                  <  2:240865/125‑1 (MQ=11)                                                                                                                                                                                                          cGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCAATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATCACCAAGCGCACTATGCGGATGCCATTCGTTATAAt                                                                                                                                  >  1:240865/1‑127 (MQ=255)                                                                                                                                                                                                              gCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTCAAAGCCCGCTGCCGCAGATATTCGCGTGGCGCGCGATACTCCCGCGCACTATGCGGATGCCATTCGTTATAATGATCGCACGCATCTGCATGTTTCTTTGCTGCCGTTAACACGTCTGGTTTGGGCATGATACTGTTGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGGCTCCGCAcc  >  1:217723/1‑251 (MQ=11)                                                                                                                                                                                                              gCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCAcc  >  1:136451/1‑251 (MQ=11)                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         AAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCGCACC  >  NC_000913/2068734‑2069189

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