Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,320,770	IS2 (+) +5 bp	coding (728‑732/1827 nt)	atoS →	sensor histidine kinase AtoS

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2068941 =	NA (NA)	43 (0.690)	41/488	0.4	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	2320770 =	0 (0.000)					coding (728/1827 nt)	atoS	sensor histidine kinase AtoS
*	?	NC_000913	= 2070271	NA (NA)	32 (0.510)	28/490	0.9	100%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 2320774	0 (0.000)					coding (732/1827 nt)	atoS	sensor histidine kinase AtoS

GATAGTTTCTACTGCAACCGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCG  >  NC_000913/2068697‑2069185                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     gATAGTTTCTACTGCAACCGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGAttt                                                                                                                                                                                                                                                >  1:196469/1‑251 (MQ=255)                    gCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGCCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTGATTGACCAACAAACATTGCAGACCTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACa                                                                                                                                                                                                                             <  2:607856/251‑1 (MQ=255)                     cAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACAt                                                                                                                                                                                                                            <  1:449685/251‑1 (MQ=255)                                      aaaaTCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTCAACAAATCCGCATTCGTGGCTGGATTTGACCCTATATTTCCCGACATCTGTTATCACTTAAccc                                                                                                                                                                                                           >  2:630566/1‑251 (MQ=255)                                                            tcgtcgTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGc                                                                                                                                                                                        >  2:334591/1‑248 (MQ=255)                                                             cgtcgtGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGa                                                                                                                                                                                     <  1:334591/250‑1 (MQ=255)                                                                                                       cccGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                                                                                                                                          >  1:548488/1‑251 (MQ=255)                                                                                                          ggCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACACCTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAtt                                                                                                                                       <  2:126702/251‑1 (MQ=255)                                                                                                           gCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGCGCGATAACCCAGCGCACTATGCGGATGCCATTc                                                                                                                                      >  1:312878/1‑251 (MQ=255)                                                                                                             gACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGt                                                                                                                                    <  1:316481/251‑1 (MQ=255)                                                                                                                          aTTGACCAACAAACATTGCAGAAATCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTCCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGa                                                                                                                        >  2:527571/1‑250 (MQ=255)                                                                                                                                                                           aTTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATg                                                                      >  1:394775/1‑251 (MQ=255)                                                                                                                                                                               gCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGGTAACCCGTCTGGTTTGGGCATGATAc                                                                  <  2:594552/251‑1 (MQ=255)                                                                                                                                                                                    aCATCCCCGTAGTAGAGATAACAGACAACGCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCATTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATg                                                             <  1:195197/251‑1 (MQ=255)                                                                                                                                                                                           cGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCCCTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                                                                                                                                          >  2:451077/1‑167 (MQ=255)                                                                                                                                                                                           cGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                                                                                                                                          <  1:451077/167‑1 (MQ=255)                                                                                                                                                                                                  gagaTAACAGACAACTCCTTGCCAGGCCGTAAACAAATCCGAATTCGTGGCTGGATTTGCCCCTCTATTTCCAGACCTCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCAGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTAt                                               >  2:360890/1‑251 (MQ=255)                                                                                                                                                                                                    gaTAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGCACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCg                                             >  2:408973/1‑251 (MQ=255)                                                                                                                                                                                                       aaCAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCAGTCTGGTTTGGGCATGATACTGATGTTGTCACGCTTTATCGttt                                          >  2:127758/1‑251 (MQ=255)                                                                                                                                                                                                           gACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTACCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCAGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCAc                                      >  1:2513/1‑251 (MQ=255)                                                                                                                                                                                                                                           cATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGAc      <  1:402437/251‑1 (MQ=33)                                                                                                                                                                                                                                               cGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTAGGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCg  <  1:630566/251‑1 (MQ=14)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     GATAGTTTCTACTGCAACCGCAACTTCTTCAGAAACAAAAATCAGCAACGAAGAGACGCTCGTCGTGACCACGAATCGTTCGGCAAGCAACCTTTGGGAAAGCCCGGCGACTATACAGGTTATTGACCAACAAACATTGCAGAACTCCACCAATGCCTCCATAGCCGATAATTTGCAGGACATCCCCGGAGTAGAGATAACAGACAACTCCTTGGCAGGCCGTAAACAAATCCGCATTCGTGGCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCCGTTACTCTCCGGACTCCG  >  NC_000913/2068697‑2069185

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