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)	30 (1.140)	26/488	0.2	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)	21 (0.800)	21/490	0.3	100%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 2320774	0 (0.000)					coding (732/1827 nt)	atoS	sensor histidine kinase AtoS

TTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATACAACT  >  NC_000913/2070029‑2070511                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 ttACGCCAGAGCAATAACTGGCTGGCTGCTACACCATGTTGCAGGGCAACGAGGGAGACCGTCAGCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATATTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGCCTATCACTTATTTAAGTGATATTGGTAGTCTGGCGATTAAGGGGGACAGTCTAGGGGCGa                                                                                                                                                                                                                                           >  2:200302/1‑250 (MQ=11)            aaaTAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCg                                                                                                                                                                                                                               <  2:211436/251‑1 (MQ=37)              aTAACTGGCTGGATGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCgtg                                                                                                                                                                                                                             >  1:204872/1‑251 (MQ=38)                    ggctggctGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGTGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGtttt                                                                                                                                                                                                                         <  1:12072/249‑1 (MQ=255)                           tgctACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAACGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTACGACTATCACTTATTTAAGTGATATTGTTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTATTTCTCAt                                                                                                                                                                                                                >  2:211043/1‑251 (MQ=255)                                                     gggAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCAGCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGTCTATCACTTATGTAAGTGATATTGTTTTTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCAGCCTCCCGTGTTTTAAGTCTCATTGATGGTCAGGAGGTAACTTATCAgc                                                                                                                                                                                      <  1:160091/251‑1 (MQ=17)                                                                              gCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTggg                                                                                                                                                             >  1:215449/1‑251 (MQ=255)                                                                                   gctgAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACt                                                                                                                                                          <  2:91186/249‑1 (MQ=255)                                                                                        aCAATTGCGATCTTTGCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATa                                                                                                                                                   <  1:57489/251‑1 (MQ=255)                                                                                          aaTTTCGATCTTTTCCTGTGTGGTACGCCGGCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGGGTTTTAATGCTCATTGATGGTCAGGAGGTAACTGATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGa                                                                                                                                                 <  2:23698/251‑1 (MQ=255)                                                                                                                                                                ccAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCaa                                                                             <  2:215449/249‑1 (MQ=255)                                                                                                                                                                   aTGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTTGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGc                                                                        >  1:145270/1‑251 (MQ=255)                                                                                                                                                                                  aCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTtc                                                         <  1:211043/251‑1 (MQ=255)                                                                                                                                                                                      gTATTCAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTGTGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTtcatc                                                      >  1:213624/1‑250 (MQ=255)                                                                                                                                                                                              aCTATCACTTAGTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATGCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAg                                             <  1:185337/251‑1 (MQ=255)                                                                                                                                                                                                                  tatTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCGCCCGTGGTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGTTGAGTCTGCGCTGGAGCGTGGTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATAACCAAAAAAGGTGAATCCCCGGACTCCCt                         >  2:204253/1‑251 (MQ=255)                                                                                                                                                                                                                   atTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGAAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCCTCACCAAAAAGGGTGAATCTCCGGACTCCCt                         >  1:146556/1‑250 (MQ=255)                                                                                                                                                                                                                             tGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGCTGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGGATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTa               >  1:263929/1‑250 (MQ=255)                                                                                                                                                                                                                                agcttCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGGGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGGGAATCTCCGGACTCCCTATATCACTTAAAtt           >  2:145186/4‑251 (MQ=255)                                                                                                                                                                                                                                         gggCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATACacct  >  1:226360/1‑248 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 TTACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATACAACT  >  NC_000913/2070029‑2070511

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